





lume XII July 1958 Number 63 


uNIVCEKSITY 
OF MICHIGAN 


APR 29 1959 


Mathematical Tables: 


and other 


Aids to Computation 





Published Quarterly 
by the 
National Academy of Sciences—National Research Council 


Published quarterly in January, Aprii, July, and October by the National Academy of Sciences— 
National Research Council, Prince and Lemon Sts., Lancaster, Pa., and Washington, D. C. 


Second-class postage paid at Lancaster, Pa. 





Editorial Committee 
Division of Mathematics 
National Academy of Sciences—National Research Council 
Washington, D. 0. 


C. B. Tompxtins, Chairman, University of California, Los Angeles, California 
C. C. Craig, University of Michigan, Ann Arbor, Michigan 

ALAN FLETcHER, University of Liverpool, Liverpool 3, England 

Eugene Isaacson, New York University, New York 3, New York 

A. H. Tavs, University of Illinois, Urbana, Illinois 

C. V. L. Smrru, Ballistic Research Labs., Aberdeen Proving Ground, Maryland 


Information to Subscribers 


The journal is published quarterly in one volume per year with issues numbered 
serially since Volume I, Number 1. Subscriptions are $5.00 per year, single 
copies $1.50. Back issues are available as follows: 


Volume I (1943-1945), Nos. 10 and 12 only are available; $1.00 per issue. 

Volume II (1946-1947), Nos. 13, 14, 17, 18, 19, and 20 only available ; $1.00 
per issue 

Volume III (1948-1949), Nos. 21-28 available. $4.00 per year (four is- 
sues), $1.25 per issue ; 

Volume IV (1950 and following), all issues available; $5.00 per year, $1.50 
per issue 


All payments are to be made to the National Academy of Sciences and forwarded 
to the Publications Office, 2101 Constitution Avenue, Washington, D. C. 


Agents for Great Britain and Ireland: Scientific Computing Service, Ltd., 7 
23 Bedford Square, London W.C.1 


Microcard Edition 


Volumes I-X (1943-1956), Nos. 1-56 are now available on Microcards and may 
be purchased from The Microcard Foundation, Box 2145, Madison 5, Wisconsin, 
at a cost of $20.00 for the complete set. Future volumes will be available on 7 
Microcards in the year following original publication. F 


Information to Contributors 


All contributions intended for publication in Mathematical Tables and Other Aids — 
to Computation and all books for review should be addressed to C. B. Tompkins, © 
Department of Mathematics, University of California, Los Angeles 24, California. © 
The author should mention the name of an appropriate editor for his paper if 

this is convenient. 











An Inverse Method for The Generation of 
Random Normal Deviates on 
Large-Scale Computers 


By Mervin E. Muller 


1. Introduction. Many applications of electronic computers require efficient 
generation of large numbers of random normal deviates; see, e.g. [7], [12], 
[14], [17]. Tables of random normal deviates are of course available, for example 
[18], [23], but they are not sufficiently extensive for many purposes and an out- 
side source of this kind cannot usually be used effectively by the computer. What 
is required is some method of generation which can be rapidly carried out by the 
machine itself. 

Methods are available by which pseudo random numbers may be produced, 
| see for example [2], [5], [6], [10], [11], [13], [20], [22]. Judging by the results 
given, for example, by [5], [6], [20], [21], the most satisfactory procedure now 
in use for generating random numbers is that based on “residue class’’ techniques, 
see for example [13], [20]. We shall not consider here the validity of these 
methods but shall assume that some satisfactory method of producing uniform 
random deviates is available from which random normal deviates are to be 
produced. Several ways of generating random normal deviates are known; see, 
for example, [1], [6], [21], but as is shown in [15], these methods are not as 
| fast nor as reliable, especially for extreme tail values, as that developed here. 

In this paper the inverse of the Normal Distribution Function is approximated 
stepwise by rational functions. Hence a uniform random deviate may be trans- 
formed into a normal deviate. The details of the approximations have been 
obtained for use on a large scale binary machine with index registers, where it is 
possible to utilize memory space in order to obtain greater speed. 

The present approach can be used, if desired, in the course of generating 
arbitrary m-variate normal deviates or Chi-squared variates. 

2. Inverse method. In principle, the inverse method of generating a normal 
deviate x from a uniform deviate u is well known. The problem is to find the 
inverse relationship x = x(u) given that 


1 z 
z= 
V2ed-~ 


| The actual determination of x(u) offers certain numerical difficulties when it is 
| desired to generate reliable normal deviates, especially for large values of x. We 
have considered this problem for two aspects. First, speed and secondly and 
| equally important, reliability. We have insisted that the maximum absolute error 
in x should be less than 4 X 10~ in the range — 5 < x < 5, where we assume 
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u is correct, and where Prob {— 5 < X < 5} > 1 — 6 X 10-’ so that x is correct 
to within 4 X 10~ except for an event of probability less than 6 X 10~’. 

The relation x = x(u) is approached in stages. The interval [0,1] for u is 
subdivided so that over each sub-interval it is possible to obtain a reliable and 
fast procedure for computing x. Over most of (0, 1) x = x(u) is approximated by 
Chebyshev polynomials. As x becomes large in absolute value it is necessary to 
increase the degree of the approximating polynomial. However, even though the 
degree of the polynomial increases, the frequency with which these approxima- 
tions are needed decreases ; hence this method will use, on the average, a low order 
approximating polynomial for x = x(u). Due to symmetry, it is actually only 
necessary to study x = x(u) for 1/2 < u < 1. 

For speed in computing it is best to have the Chebyshev polynomials of as 
low degree as possible, subject to the specified level of accuracy. Following this 
approach one could find sub-intervals of greatest or optimum width, see [3]. 
Two drawbacks then become apparent. If on a binary computer the sub-interval 
widths are not some negative power of two, considerable computing time is spent 
in mapping a given sub-interval onto the range (—1, 1) which is needed when 
using a Chebyshev approximation. More important, if the intervals are not of 
uniform width and a negative power of two, considerable machine time, and 
memory space, is required to select the appropriate approximation. For, it will 
not then be possible to use an index register to select the appropriate approxima- 
tion. Taking these facts into consideration, and taking into account memory space 
requirements, it was decided to have the widths of the sub-intervals equal 1/128. 
With this choice of interval width the largest part of (0,1), namely 1/8 < u < 7/8, 
could be approximated by linear functions. Quadratic and quartic approximations 
are used for 7/8 < u < 127/128. For 127/128 <u< 1, and by symmetry 
0 <u < 1/128, x = x(u) has a singularity of logarithmic type; consequently, 
for u in this sub-interval an approximation of more subtle type than Chebyshev 
polynomials is needed. 

3. Use of polynomials. Since the techniques used in approximating x = x(u) 
have possible application when considering the generation of other possible 
pseudo random variables we shall review the techniques in detail. 

To approximate x = x(x) in the jth sub-interval, 


64+j-1 “+i) : 
——_—_————_ < < = 
( 128 oem oe yf 1(1)63, 





by Chebyshev polynomials we shall make use of the technique of Interpolating 
Chebyshev Polynomials as developed by C. Lanczos, see for example [8], [9]. 
To begin, the jth sub-interval is transformed onto (—1, 1) by using 


r = 256u — 127 — 2}. 


Let r = cos 0,0 < @ < x. Let T;(r) = cos [i(cos™ r)] denote the ith Chebyshev 
polynomial. Thus for the jth sub-interval, x = x,(r) will be approximate by a 
polynomial of mth degree, namely, 


Xaj(r) = 1/2 do; + e dijT;(r) ; 
t=] 
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where the d;;’s are obtained from 


2 n 
di; = ei | X #(re) Til), 





and where the r;’s are the zeros of the (m + 1)th Chebyshev polynomial. Thus 
te= cos (x(2k + 1)/2(n + 1)). This approach approximately minimizes the 
maximum error throughout the sub-interval. 

For computing purpose it is more convenient to express the 7;(r)’s as poly- 
nomials in r, see Lanczos [8], [9], for example 7:(r) = r, T:(r) = 2r? — 1. One 
then obtains 


(1) Xnj(r) = E asp 


The coefficients a;; are given in section 4. 

In calculating the d;;’s, and hence the a;;’s, it is necessary to evaluate x = x(u) 
for u;’s where 

r, + 127 + 27 
uj = 256 ’ 

and this was done for the seven intervals 7 = 57(1)63, i-e., 7/8 < u < 127/128. 
The values of x(uz;) were obtained by cubic inverse interpolation in the National 
Bureau of Standards Tables of the Normal Probability Function, [19]. 

However, for the intervals, 7 = 1(1)56, i.e., 1/2 < u < 7/8, where the linear 
approximations are adequate the essential computations for the a,;;'s are concerned 
with finding the inverse values, namely x = x(u,;). However, for the linear cases 
the a;;’s given in section 4 were not obtained in the above manner since we had 
previously needed to have values of x corresponding to u, = (128 + s)/256, 

= 0, 1, 2, ---, 128. Consequently we avoided the additional labor of finding the 

necessary x’s by utilizing the available x(u,)'s. These values are included as Table 
5 in the Appendix. Thus, while not fitting the straight lines by first degree 
Chebyshev polynomials we were able to obtain the desired level of precision as 
follows: In the jth interval, where it is appropriate to fit a straight line, consider 
Uj, Uja1, ANd Uj, 12 = (Uz + 4541)/2. Let xj, Xj41, Xj, 1/2 denote the corresponding x 
values. The jth line is fitted such that the deviations at the ends of the class 
interval are equal and such that the deviation at the mid-point, namely 1;, 1/2, is 
equal in magnitude to the deviations at the end points, but is of opposite sign. 
This approach gives essentially a Chebyshev-type approximation in that we are 
striving to minimize the maximum error in a given sub-interval. By straight- 
forward computations one then obtains that 


2 (xs44 + X5,1/2) — (Xi41 — x5) (255 + 47) 
ay = 4 ‘ 
ay = 128 (x54: = %3). 





In the selection of the width of a sub-interval, w;, or for determining the appro- 
priate degree of the Chebyshev Polynomial for a given sub-interval, the following 





result is used, see [3] for more general details, 


1 
w; = 4 E= €|x*) (c;) a D 
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where « is the largest absolute error that is tolerable, in our case « = 4 X 10-4, 
and where x‘"*» (u) is the (7 + 1)st derivative of x = x(u) with respect to wu. 

4. Table of coefficients. The following tables provide the necessary coeffi- 
cients for the approximations given by equation (1) of section 3. 


j 40; 
1 — 1.25339 
2 — 1.25388 
3 — 1.25487 
4 — 1.25639 
5 — 1.25845 
6 — 1.26107 
7 — 1.26428 
8 — 1.26809 
9 — 1.27255 

10 — 1.27767 

11 — 1.28348 

12 — 1.29003 

13 — 1.29735 

14 — 1.30547 

15 — 1.31445 

16 — 1.32432 

17 — 1.33514 

18 — 1.34698 

19 — 1.35986 

20 — 1.37390 

21 — 1.38914 

22 — 1.40567 

23 — 1.42357 

24 — 1.44295 

25 — 1.46390 

26 — 1.48655 

27 — 1.51103 

28 — 1.53749 

J 

57 
58 
59 
60 
61 
62 

j Qo; 


449 
344 
723 
365 
203 
329 
016 
727 
126 
105 
794 
594 
187 
572 
096 
485 
882 
984 
540 
816 
726 
396 
635 
134 
578 
769 
773 
093 


a0; 


1.56668 
1.63732 
1.71722 
1.80989 
1.92135 
2.06352 


TABLE 1. Linear Cases 


a1; j Qo; 
2.50678 851 29 — 1.56607 
2.50775 044 30 — 1.59698 
2.50967 688 31 — 1.63039 
2.51257 300 32 — 1.66656 
2.51644 669 33 — 1.70571 
2.52130 843 34 — 1.74815 
Z.a00k4. 192 35 — 1.79421 
2.53405 216 36 — 1.84424 
2.54196 946 37 — 1.89869 
2.55094 571 38 — 1.95803 
2.56100 644 39 — 2.02283 
2.57218 075 40 — 2.09375 
2.58450 135 41 — 2.17155 
2.59800 495 42 — 2.25713 
2.61273 256 43 — 2.35154 
2.62872 974 44 — 2.45607 
2.64604 704 45 — 2.57223 
2.66475 763 46 — 2.70077 
2.68485 495 47 — 2.84833 
2.70651 014 48 — 3.01223 
2.72973 042 49 — 3.19590 
2.75461 642 50 — 3.40829 
2.78126 053 51 — 3.65403 
2.80976 489 52 — 3.93953 
2.84024 264 53 —4.27256 
2.87281 914 54 — 4.67044 
2.90763 359 55 — 5.15058 
2.94484 078 56 — 5.74560 
TABLE 2. Quadratic Cases 
aij 
859 —0.03343 48405 
538 —0.03745 15701 
812 —0.04284 26652 
233 —0.05049 00254 
077 —0.06226 30013 
790 — 0.08299 31005 
TABLE 3. Quartic Case 
a2; a3; 


aij 
63 +2.26622 681 +0.12757 931 +0.01844 432 +0.00424 42872 +0.00104 06032 


a1; 
868 2.98461 
101 3.02714 
945 3.07264 
010 3.12136 
765 3.17357 
989 3.22957 
327 3.28972 
956 3.35441 
401 3.42410 
516 3.49930 
720 3.58062 
508 3.66874 
347 3.76449 
085 3.86881 
997 3.98282 
734 4.10785 
456 4.24551 
415 4.39645 
659 4.56815 
276 4.75713 
941 4.96703 
673 5.20760 
652 5.48350 
404 5.80124 
248 6.16869 
630 6.60394 
375 7.12472 
472 7.76467 

a2; 
.00087 5575 
.00114 804 
.00157 546 
.00230 549 
.00372 027 
.00708 977 
a4; 


325 
356 
746 
699 
468 
838 
698 
736 
298 
401 
028 
687 
389 
137 
083 
604 
603 
556 
521 
372 
562 
458 
516 
828 
364 
516 
205 
397 














mG 


25 
56 
46 
09 
68 
38 
98 
36 
8 
)1 
28 
37 
39 
37 
33 
)4 
)3 
6 
1 

2 

2 

8 

6 
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5 
7 
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5. The interval (127/128 < u < 1). In view of the fact that x(u) may be 
looked upon as a function having a singularity of logarithmic type in this interval, 
it is necessary to abandon the use of a polynomial type approximation here. 
A satisfactory rational approximation is obtained by using a truncated continued 
fraction expansion, see for example [4], [16]. For any given value of u the 
corresponding value of x is approximated by the following recurrence relation. 
The kth convergent to x() is given by 


Mi (u) 
Ni(u) i 





Xk (u) a 


where M;,(u) and N;(u) are determined as follows: 
(2) Mis1(u) = ad, M;(u) a (u = Ux—1)My_-1(u); M,(u) = i; M,(u) = dg; 
(3) Nesi(u) = deNi(u) + (u — ue-s)Ne-i(u); No(u) = 90; Ni(u) = 1; 


and where d; denotes a selected value of the kth inverted divided difference of x(x). 

In order to insure sufficient precision it is suggested that the machine program 
test to see if wu > .99999. When x is greater than .99999 it is suggested that the 
additional significant figures for u,, k = 10(1)14, be utilized. Double precision 
operations are not necessary if these u, are stored with the first three nines sup- 
pressed, « must then be shifted the appropriate number of places before per- 
forming u — t,. 

Following [4], page 406, x.(u) can be computed from 





: KF (—1)*#*(u — uo) (u — m1) --- (uw — Un—1) | 
Pi, es eae No(1s) Noxa(t) 


To insure that x(u) can be approximated to within 4 X 10~ it is necessary to 
stop at an appropriate value of k = k(u) since the approximation being employed 
is essentially a divergent expansion. Numerical evaluations verified that it is 
sufficient to have the machine program terminate as soon as one of the following 
three conditions is satisfied, namely: (i) select x,(u) as the approximation to x(u) 


TABLE 4. Constants for Rational Approximation 


k Us d, 

0 .99223 97464 2.42000 0000 

1 .99461 38540 0.01826 2366 

2 .99653 30262 —0.65518 1080 

3 .99813 41867 0.02399 7757 

4 .99903 23968 —0.26737 1460 

5 .99931 28620 0.01654 1263 

6 .99966 30707 —0.14194 9840 

7 .99984 08914 0.00997 32778 
8 .99992 76519 —0.06004 9158 

9 -99996 83287 0.00511 81541 
10 -99998 66542 Si —0.02329 9296 

il .99999 45874 56 0.00241 07770 
12 .99999 78875 45 —0.00863 34192 
13 .99999 92066 72 0.00100 76316 


i4 .99999 97133 484 —0.00308 28145 
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if r~ = |xe(u) — xe-1(u)| < 4 X 10+, (ii) select x,(u) as the approximation to 
x(u) if re — Tex. < 0, and (iii) if k = 14 terminate the process and select x14(u) 
as the approximation to x(x). 

The value of k increases as u increases in the interval (127/128 <u <1 
— 3 X 10’), and consequently so does the computing time. However, this sub- 
interval will, on the average, increase the necessary computing time very little. 
It should be kept in mind that while this method is very appropriate and con- 
venient for a machine which performs “floating point’? multiplication, this ap- 
proach would create serious scaling difficulties for a ‘fixed’ point’? mode of 
operation. 

6. Concluding comments. A detailed comparison of this method with other 
known methods is given in [15]. The technique proposed in this paper yields 
accuracy comparable with, or better than, most previous proposals using about 


one-quarter the computing time while requiring about twice as many memory 
locations. 

The author wishes to express his appreciation to Dr. G. E. P. Box, Dr. Collin 
Mallows, and Dr. John W. Tukey for their interest and generous comments 
concerning this work. The author is also indebted to Mrs. A. Schay for her ex- 
cellent help with the numerical computations. 


Princeton University 
Princeton, New Jersey 
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Appendix 
This table gives values of x(u) corresponding to the Normal Distribution for 
F(x;) = 1/2 + 7/256, 7 = 1(1)127. The values were obtained from the National 


Bureau of Standard Tables of the Normal Probability Function, [19], by using 
the inversion formula: 


x(u) = we + +4(2) A a pA 


where up is the nearest tabulated entry to u, and u — uy = a, and where 


Qo = 1 e—=(uo)]*/2, 
2x 
TABLE 5. Inverse Values for the Normal Distribution 
j F(x;)=1/2+ 7/256 Xj j F(x) =1/24+-7/256 Xj 
0.50390 625 0.00979 167 | 33 0.62890 625 0.32895 791 


0.50781 250 0.01958 429 | 34 0.63281 250 0.33931 161 
0.51171 875 0.02937 878 | 35 0.63671 875 0.34970 180 
0.51562 500 0.03917 609 | 36 0.64062 500 0.36013 003 
0.51953 125 0.04897 716 | 37 0.64453 125 0.37059 729 
0.52343 750 0.05878 294 | 38 0.64843 750 0.38110 545 
0.52734 375 0.06859 437 | 39 0.65234 375 0.39165 587 
0.53125 000 0.07841 241 | 40 0.65625 000 0.40225 007 
0.53515 625 0.08823 802 | 41 0.66015 625 0.41288 960 
0.53906 250 0.09807 215 | 42 0.66406 250 0.42357 608 
0.54296 875 0.10791 578 | 43 0.66796 875 0.43431 116 
0.54687 500 0.11776 987 | 44 0.67187 500 0.44509 652 
0.55078 125 0.12763 542 | 45 0.67578 125 0.45593 392 
0.55468 750 0.13751 340 | 46 0.67968 750 0.46682 512 
0.55859 375 0.14740 482 | 47 0.68359 375 0.47777 199 
0.56250 000 0.15731 068 | 48 0.68750 000 0.48877 641 
0.56640 625 0.16723 201 | 49 0.69140 625 0.49984 034 
0.57031 250 0.17716 982 | 50 0.69531 250 0.51096 581 
0.57421 875 0.18712 516} 51 0.69921 875 0.52215 488 
0.57812 500 0.19709 908 | 52 0.70312 500 0.53340 971 
0.58203 125 0.20709 265 | 53 0.70703 125 0.54473 251 
0.58593 750 0.21710 695 | 54 0.71093 750 0.55612 559 
0.58984 375 0.22714 306 | 55 0.71484 375 0.56759 132 
0.59375 000 0.23720 211 | 56 0.71875 000 0.57913 216 
0.59765 625 0.24728 522 | 57 0.72265 625 0.59075 066 
0.60156 250 0.25739 353 | 58 0.72656 250 0.60244 945 
0.60546 875 0.26752 821 | 59 0.73046 875 0.61423 129 





bo DO DO DO DO DO DO DO DO FF FO tt tt 





0.60937 500 0.27769 044 | 60 0.73437 500 0.62609 901 
0.61328 125 0.28788 143 | 61 0.73828 125 0.63805 558 
0.61718 750 0.29810 241 | 62 0.74218 750 0.65010 407 
0.62109 375 0.30835 463 | 63 0.74609 375 0.66224 768 
0.62500 000 0.31863 936 | 64 0.75000 000 0.67448 975 








j F(x) =1/245/256 


65 0.75390 
66 0.75781 
67 0.76171 
68 0.76562 
69 0.76953 
70 0.77343 
71 0.77734 
72 0.78125 
73 0.78515 
74 0.78906 
75 0.79296 
76 0.79687 
77 0.80078 
78 0.80468 
79 0.80859 
80 0.81250 
81 0.81640 
82 0.82031 
83 0.82421 
84 0.82812 
85 0.83203 
86 0.83593 
87 0.83984 
88 0.84375 
89 0.84765 
90 0.85156 
91 0.85546 
92 0.85937 
93 0.86328 
94 0.86718 
95 0.87109 
96 0.87500 


625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 
000 


TABLE 5—Continued 


xj 


0.68683 
0.69928 
0.71184 
0.72451 
0.73730 
0.75021 
0.76325 
0.77642 
0.78972 
0.80317 
0.81676 
0.83051 
0.84441 
0.85848 
0.87272 
0.88714 
0.90175 
0.91655 
0.93156 
0.94678 
0.96222 
0.97789 
0.99381 
1.00999 
1.02643 
1.04315 
1.06018 
1.07750 
1.09518 
1.11319 
1.13157 
1.15035 


375 
330 
220 
438 
400 
538 
304 
176 
652 
257 
542 
088 
508 
447 
589 
656 
411 
667 
283 
176 
320 
754 
591 
017 
306 
826 
048 
557 
065 
428 
656 
938 





J 
97 
98 
99 
100 
101 
102 
103 
104 
105 
106 
107 
108 
109 
110 
111 
112 
113 
114 
115 
116 
117 
118 
119 
120 
121 
122 
123 
124 
125 
126 
127 


F(x;) =1/2+ 3/256 


0.87890 
0.88281 
0.88671 
0.89062 
0.89453 
0.89843 
0.90234 
0.90625 
0.91015 
0.91406 
0.91796 
0.92187 
0.92578 
0.92968 
0.93359 
0.93750 
0.94140 
0.94531 
0.94921 
0.95312 
0.95703 
0.96093 
0.96484 
0.96875 
0.97265 
0.97656 
0.98046 
0.98437 
0.98828 
0.99218 
0.99609 
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625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 
000 
625 
250 
875 
500 
125 
750 
375 





xj 


1.16953 
1.18916 
1.20926 
1.22984 
1.25099 
1.27268 
1.29502 
1.31801 
1.34171 
1.36620 
1.39153 
1.41779 
1.44507 
1.47345 
1.50310 
1.53412 
1.56668 
1.60100 
1.63732 
1.67594 
1.71722 
1.76167 
1.80989 
1.86273 
1.92135 
1.98742 
2.06352 
2.15387 
2.26622 
2.41755 
2.66006 


661 
435 
123 
876 
172 
865 
241 
090 
784 
382 
749 
714 
258 
903 
294 
054 
859 
866 
538 
192 
812 
041 
233 
187 
077 
789 
790 
469 
681 
902 
747 


Algebraic Approximations for Laplace’s Equation 


in the Neighborhood of Interfaces 


By J. W. Sheldon 


1. Introduction. Let it be required to solve the following problem: 


Problem A: 


Let C;,7 = 1, 2, be two simple, closed plane curves with continuous curvature. 
Let C; enclose all the points of C;. Let G: be the region interior to C;. Let Gz be 
the region interior to C; and exterior to C,. Let W be a continuous, bounded func- 
tion of position on C2. Let it be required to find harmonic functions V“ such that 


(a) V® = W on Cz. 
(b) V® is regular and bounded in G,. 
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(c) V® = V® on Cy. 
avm ave 


an on 
along the normal to C, the positive normal direction being chosen as the direction 
from G; into Gz. 








te) 
(d) Di on Ci, where D; > 0 and aa denotes differentiation 


Problems with boundary conditions (c) and (d) connecting distinct harmonic 
functions on a common boundary between adjacent regions arise in electrostatics, 
and specific cases are treated in books dealing with electrostatic theory. (See, for 
example, Jeans [1] and Smythe [2].) For electrostatics, D; and Dz are the 
dielectric constants in G; and Gz. Boundary condition (c) in electrostatic theory 
is derived from the requirement that no work be done taking a charge around a 
closed path which crosses C,. The required condition is found to be V = V® 
+‘‘constant” and we set “‘constant” to zero. Boundary condition (d) is derived 
from Gauss’s electric flux theorem. Boundary condition (d) may also be derived 
by a variational method. If we require that D, times the Dirichlet integral for G, 
plus D times the Dirichlet integral for G; be an extremal, and that V = V® on 
C;, then (d) is obtained as the natural boundary condition on C;. 

Problem A may also arise for steady-state heat conduction, steady-state 
diffusion, steady-state current flow, and for problems of incompressible, immiscible 
fluid flow in porous media. In the latter problem C, may be a moving curve [3]. 
Boundary conditions (c) and (d) may also arise in transient problems of dif- 
fusion type. 

We suppose that we wish to solve Problem A by finite difference methods. We 
introduce a rectangular network over a region covering G; and G2. We introduce 
a function g, defined at each meshpoint interior to C2, which is to approximate in 
some way the solution of Problem A. In order to determine ¢ we require a system 
of equations equal in number to the number of meshpoints interior to C2. We call 
each equation in the system an “algebraic approximation.’ In the way that we 
will treat this problem, there will be an algebraic approximation associated with 
each meshpoint inside C, in the sense that each algebraic approximation approxi- 
mates to the conditions of Problem A best in the neighborhood of a unique mesh- 
point and no two algebraic approximations approximate Problem A best in the 
neighborhood of the same meshpoint. The meshpoints interior to C; may be 
classified in the following way: 


(i) Meshpoints adjacent to C2. A meshpoint is adjacent to C; if a straight-line 
segment from the meshpoint to one of the four nearest neighbor meshpoints 
intersects or extends to C>. 

(ii) Meshpoints adjacent to C;. A meshpoint is adjacent to C; if the meshpoint 
is situated on C; or if a straight-line segment from the meshpoint to one of the 
four nearest neighbor meshpoints intersects C, (but not if the segment only 
extends to C). 

(iii) Normal meshpoints. All other meshpoints interior to C; are called normal 
meshpoints. 


The algebraic approximation to be associated with the normal meshpoints is 
the well-known five-point formula for approximating the Laplace equation. 
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Algebraic approximations at the meshpoints adjacent to C, have been obtained. 
See, for example, Shortley and Weller [4] or Viswanathan [5] (Algebraic ap- 
proximations have been obtained for Dirichlet or Neuman or mixed conditions 
on C:.) The present paper is concerned with obtaining algebraic approximations 
for meshpoints adjacent to Ci. 

Problems with boundary conditions (c) and (d) of Problem A have been 
treated by finite difference methods in connection with the ‘‘group-diffusion”’ 
problems of reactor design [6, 7]. To the author’s knowledge, in this field it has 
always been assumed that the interfaces between materials have contours such 
that analogues of equations (21) or (22) apply. If the interface does not meet 
the conditions for equations (21) or (22), it has been customary to perturb the 
interface so that it does. In some problem this may be an unsatisfactory procedure. 
In the case of ‘“‘moving boundary”’ problems, for example, the essence of the 
problem consists in getting the velocity of the moving boundary accurately. 

2. Existence of solution to Problem A. Books on potential theory which give 
existence proofs for the Neuman and Dirichlet problems do not (to the author’s 
knowledge) treat Problem A. An existence proof can be constructed following the 
integral equation technique described by Kellogg [8]. This method consists in 
postulating charge and dipole distributions on the boundary curves, and then 
obtaining and proving the existence of a solution of the equations which must be 
satisfied if the charge and dipole distributions are to produce potential functions 
which satisfy the given boundary conditions. The equations obtained are linear 
integral equations. The method can be extended to show that Problem A possesses 
a solution. It can be shown that V® is the potential of a continuous dipole dis- 
tribution on C, plus the potential of a continuous charge distribution on Cy, 
V® is the potential of a continuous charge distribution on C;. (The charge dis- 
tribution on C; to produce V™ is different from that on C; which contributes to 
produce V®.) We shall not give the detailed proof. In view of the fact that 
Problem A is the solution of a problem in the calculus of variations, and that we 
have required C; and C: to be smooth and W to be continuous, it would be very 
surprising if the solution did mot exist. Hence we would not really be justified in 
presenting a rather long existence proof here. 

3. The behavior of V“ and V® on C;. Let us suppose that C; is an analytic 
arc or that C; is composed of segments of analytic arcs. We have the following: 


THEOREM 1. A segment T of C; which is an analytic arc is not a natural boundary 
for V® or V®, 

This means that every point of I has a neighborhood such that there exists 
a harmonic function which is regular in the region consisting of G, and this 
neighborhood which coincides in G; with V™. And every point of I has a neighbor- 
hood such that there exists a harmonic function which is regular in the region G2 
and this neighborhood and which coincides in G; with V®. 


In order to prove Theorem 1 we first solve the following problem : 


PROBLEM B. Let K be a circle of radius R and center at the origin of the x, y 
coordinate system. Let G, be the region interior to K lying below the x-axis and 
let Gz be the region interior to K lying above the x-axis. Call the part of the x-axis 
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interior to K, ['. Let r, 6 be polar coordinates in the x, y plane. Let W™ (6), 
—2r<26< 0, W(6),0 < 0 < x be given bounded continuous functions of @ on K 
5 tle aw” o Iwo 
satisfying W(0) = W®(0), W®(—2) = W(x), pa = pr”, 
i (2) 
p, 2” (—7) a p,2@ (x) 
00 00 
(a) V®(R,0)=W%), —x<0<0, 
V@(R, 0) = W (6), 0<0<-r. 
(b) V® is regular and bounded in G®. 


. Let it be required to find functions V such that 


(c) vo (r,0) = V® (r, 0) on f. 
(d) p, 2729 Le p, 2229 on #. 
06 06 


Problem B can be solved by the method of circular harmonics. We obtain: 


(1) V@(r,0) = 440+ > a.(Z) cos 26 


n=l 


+£3,(5) sinn#, —x<0<0, 


n=1 


(2) V(r,6)=4404+ 5. 4.(Z) cos at 


n=l 


+23 3(5) 8 no, O<0< 
D.& R sinné, O<@0<z, 


(3) An = Died. W® (6) cos méd@ 
1 <6 


oe f (2) 
+ (Di + DarJo W® (6) cos médé@, 


= ol ©) (9) si fe —e | 
sy (D: + D:)xLJ-« W® (6) sin m6dé@ + § W® (0) sin méd@ | , 
m=Q0,1,2--- 


(4) Bn 


With (3), (4) substituted into (1) andr = R, it is apparent that (1) represents 
the sum of two Fourier cosine series plus two Fourier sine series. The formal sums 
of the two series are easily identified. The two cosine series give us 


D; 


—xr<60<0. 
D. + D. WO +5 W?(—8), x<0<0 
The two sine series give us 
D: wn 
—7r<6<0. 
bap” = D+ D, wW?(—8), x™<e<0 


Adding the above terms, the sum of the two cosine series plus the two sine series 
gives us W1(6), — x < 0 < 0. W'(6), W*(6) satisfy Dirichlet conditions and are 
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furthermore continuous within their range of definition. This means that the 
cosine series converge to their formal sums for — + < @ < 0, and the sine series 
converge to their formal sums for — x < @ < 0. The sine series converge to the 
value zero, at 6 = 0, 0 = x. The sum of the four series gives W'(0) for —r<0<0. 
A similar argument shows that if r = R the series on the right-hand side of (2) 
converges to W*(6) for 0 < 6 < rz. It is readily verified that V™(r, 0), V(r, 6) 
as defined by the series in (1), (2) satisfy the remaining conditions of Problem B 
so that (1), (2) is the solution of Problem B. Then we have: 

Lemma 1. Lf is not a natural boundary for V(r, 6) or V(r, @) of Problem B. 

The series in (1) and (2) are each uniformly convergent for r < p < R for all 
@ and hence each of (1), (2) defines a harmonic function throughout the interior 
of K. Q.E.D. 

It is, of course, a well known theorem from function theory that under a 
conformal transformation a harmonic function maps into a harmonic function. 
We have: 

Lemma 2. Let a neighborhood of C; be mapped into a neighborhood of C., where 
C, is the mapping of C., by a conformal transformation. Let V map into V. 

av av 
Then D,— = D,— on €.. 
on on 

The Lemma follows from the definition of a conformal transformation. 

With the help of the solution to Problem B and Lemmas 1, 2 we can prove 
Theorem 1. The method of proof follows closely a proof by Bieberbach [9] of 
a lemma from potential theory to the effect that if V were regular in a region G, 
continuous in the extended domain obtained by adding to G the points of an 
analytic arc that belongs to the boundary of G and if V were zero on this arc, 
then this arc is not a natural boundary for V. We may consider V“ to be the 
real or imaginary part of a function of a complex variable. Let the analytic arc T 
be represented in the complex plane by z(a), a < a < b, where a is the real part 
of a complex variable y = a + 78. a is also chosen as the arc length along I. 

Then the function z(7) maps a neighborhood K in the y-plane of every point 
B = 0, a < ao < 5 onto a neighborhood U of the point z = 2(ao); this mapping 
is conformal. Since a represents arc length along I x’ (ao)? + y’(ao)*? = 1. Hence 
2’ (ao) # 0. Hence the mapping is also simple. We choose for K a circle with center 
on the real axis at ao. This determines a region U in the z-plane. Let G; and V™, 
4 = 1, 2, refer to the regions and functions of Problem A. From the intersection 
of G; and the boundary of U we obtain the values of V“ on the boundary of U. 
By the inverse of the mapping mentioned the boundary values on U are trans- 
formed to boundary values on K. These satisfy the conditions for Problem B. 
From the results of Problem B each of the transforms of the functions V™ to K 
is analytic in K, and hence by an inverse transformation each of the V“ is analytic 
throughout U. Q.E.D. 

The domain of analyticity of the functions continued across T is difficult to 
determine. It will depend on the location of the zeros and poles of z(7), as well 
as on the location of non-analytic points of C:, and the distance of C; to C2, and 
also the location of the mapping of Co. 

4. The connection between power series for V and V. Let O be the center 
of the circle K. Let ¢ and m be the abscissa and ordinate measured from O. We 
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label the value of V“ at O, V§? and we let Vi, V&, VP, etc., be partial deriva- 
tives of V“ with respect to ¢ and at O. We may expand V“ about O in power 
series convergent in K: 


(5S) VO = HP + Vin + Wit + Viet + Vint + AVP! + ---, 
(6) V® = 09 + 19 n+ V+ Vn + Vat + 40! + ---, 


THEOREM 2. If the power series for V™ is given the coefficients in the power series 
for V™ may be determined. 


Since V® = V® on PL it follows that 


(7) VP = Vf, VF = Vf. 
Since D.V® = D,V%, it follows that 
(8) DV? = DV? 


Since 0}? = Vi, WY + V2 = 0, i = 1, 2, it follows that V2 = V°2. We 
have expressed V and all its derivatives through second order in terms of 1 
and its derivatives through second order. Let § = D:/D2. We have 


(9) V® = VP + Vn + Wt + AV? + EVP nt + 40 ue’. 


We shall be concerned with terms through second order only. However, one 
may show that (7), (8) and the derivatives of the equations Vj) + Vi2=0 
always give enough equations so that the derivatives of V can be expressed in 
terms of those of V. 

5. Five point algebraic approximations. We consider first the case in which 
T is, locally, a straight line. We suppose that I passes through a meshpoint making 
an angle a with the x-axis. Let the meshpoint be labeled “a” and let the four 
nearest neighbors be labeled “‘b,” ‘‘c,” “‘d,”’ “e’’ as shown in figure 1. From 
Theorem 1 it follows that V“ should satisfy Laplace’s equation at “‘a.’"’ Let the 





C=-hsin« 
N=-hcos~ 
+ 
Cc “& 
a 
Ps L=-hcos« 
beh ox 7 N=hsinw 
= hcosx ™ X 
N=-hsinex 4 a b 
t=hsinc 2 
N= hcosx ? G 
Gi \x 
Ficure 1. 
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network values of V“* at a, b, c, d, e be Va, Vs, Ve, Va, Ve. It is not necessary to 
superscript these quantities since there is only one function value associated with 
each-meshpoint. Let 4 be the mesh spacing. We have the following theorem: 


THEOREM 3. When £ + 1, except for a = 7 ,n=0,1, ---,7, there exists no 


algebraic approximation for Laplace’s equation in terms of Va, Vs, Ve, Va, Ve which 
ts such that the local truncation error tends to zero as h > 0. 


For definiteness, assume that 7/2 < a < x. Then meshpoints 6} and ¢ are on 
the boundary of or inside Gz and d and e are on the boundary of or inside G. 

We choose an origin at ‘‘a’’ and introduce coordinates ¢, m where n is the normal 
at “a” from G; into G2 and ¢ is the distance along I from a. The direction of 
positive ¢ is chosen so that the #, m coordinate system is right-handed. It is as- 
sumed that a circle about a which just encloses }, c, d, e has radius less than the 
circle of convergence of the power series (5), (6) and that I is straight inside this 
circle. Since IT is straight it is not necessary to distinguish between T, f; V“, 
V, etc. The values assumed by ¢ and at b, c, d, e are given in figure 1. Partial 
derivatives V,, V:, etc., refer to V,™, V.™, etc., in what follows. Making use of 
(5), (9) we have 


(10) Via-— Va = thsinaV, — hcosaV; + $h? sin? aVan 
— th’ sin a cos aVn; + $h? cos? aVee + O(h'), 


(11) V.-— Va = — theosaV, — hsinaV; + $h? cos? aVan 


+ th? sin a cos aVa: + $h* sin? aVig + O(h*), 


(12) Va— Va = —hsinaV, + hcosaV; + $h’ sin? eVny 
— sin acos aVnz + $h* cos? aVix + O(h*), 


ll 


hcosaV, + hsinaV; + $h? cos? aVan 
+ h* sin a cos aVae + 3h? sin? aVy, + O(h'), 


(13) Ve-— Vu 


(14) O Van + Vie 


ll 


Suppose we could solve (10)—(13) for V,, and Vi. By substituting the solution 
into (14) we would have an algebraic approximation to Laplace’s equation. The 
O(h*) terms would be O(h) after carrying out this process and would tend to zero 
as h +0. A necessary condition for being able to carry out this process is that 
the determinant of coefficients of (10)—(14) vanish. In order to evaluate the de- 
terminant, we take the following combinations of equations (10)—(14) : 


Eq. 15 = Eqs. [(10) + (11) + (12) + (13) — n2(14)], 


i+é 
2 





Eq. 16 = Eqs. [ (10 + (11) + &(12) + &(13) —- rca) | , 


Eq. 17 = Eggs. [an — 4 sin? ah?(14) + €(12) — * cos" ait(14) | ; 


Eq. 18 = Eqs. [(11) + (12) — 3#°(14)], 
Eq. 19 = Eq. (14). 








to 
th 


10 


ch 


on 
he 


at 
le- 








ALGEBRAIC APPROXIMATIONS FOR LAPLACE’S EQUATION 181 


The transformation of (11)—(14) to (15)—(19) is non-singular for — ~ 1. We 
obtain 


(15) (& — 1)(sina — cosa)hV, = Vit+ Ve + Va + Ve — 4Va + O(h'), 


(16) (€ — 1)(sina + cosa)hV, = Vi + V. — 2Vz 
+ &(Va + Ve — 2V.) + O(h'), 


(17) —é&(sina + cos a)hV, + (Ecosa — sin a)hV, 


1-¢ 
+ 2 





(cos? a — sin* a)h?V,, = Ve — Va + &(Va — Va) + O(h*), 


(18) —(sina + écosa)hV, + (cosa — sina)hV, 
+ (& — 1) sina cos ah? V,, = V. + Va — 2V, + O(h*), 
(19) Fim a Vu = 0. 


The equation system (15)—(19) is triangular. Its determinant D is 


(20) D = 3(& — 1)4(sin? a — cos? a)? sin a cos ah’. 
3 
D is zero for § # 1, 5 Sa < x only whena =5,a=—,a=-. By sym- 
metry, D would be zero also only for a = r¥ n=0,1,---,7 forO <a < 2z. 


. — . 3 
This proves Theorem 3. The algebraic approximations which we get, for a = = 


and a = = are, respectively : 


(21) ve + vo —-—22 [ = 2% Ye 

















Di+D: h* 
vig =H 00 
(22) VP + vg=4= a TN a, " r “*) 
+ af ms = ve) + O(h) = 0 


(21) and (22) may also be obtained by an application of Gauss’s theorem 
[6, 7]. 


If T is straight but does not pass through meshpoints, we still obtain an alge- 
" oe i : : nx 
braic approximation in terms of four neighbors if a = 3° For example, suppose 


a = w and I passes a distance 9 above “a” leaving “a” in G, “c’’ in Gs. We 
choose an origin O on [ where I intersects ac. We let V;, V,, etc., represent 
derivatives of V at O. The continuation equations analogous to (10)—(13) are 


(23) Vo a Vo = —1Va + dn? Van + O(h'), 


(24) Vi — Vo = —nV. +hVi+ 31? Van 7 hn Vas + 5h V4 + O(h'), 
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(25) Ve- Vo = &(h — )V, + (hk — 2)?Van + O(h), 
(26). Va — Vo = —9Vn — hVi + 49°Vanw + hn Vat + $2°Vie + O(h!), 


(27) Ve — Vo = —(n + A)Vn + (A + 20)?Van + O(h'). 
Subtracting (23) from (25) we obtain 


(28) Ve — Va = [&(A — 9) + 9) Vn + 3 (0? — 2h) Van + O(h*). 
Subtracting (27) from (23) 
(29) Va — Ve = hVn — $(2nh + h?) Van + O(h*). 


Dividing (28) by (4 — ») + 7 and subtracting (29) divided by h, we obtain 








Ve — Va Va — V. | h? — 2hn | 
30 ~ = =) ——— +h + 2] Van + OCF). 
me a-d+s 8 "Hib o+s . ” 
Taking (24) minus twice (23) plus (26) and then dividing by h?, we have 
Vi — 2Va+ V 
(31) — SSF 6 Ve + OU: 


Dividing (30) by the factor multiplying V,. in (30), then adding equation (31) 
and making use of V,, + Vi = 0, we cbtain 


2 ie Ve — Ve . e 
(32) Van + Vu = 2[ P= 4 5 + 20] [ _ 4 








t(n—n) + th-n) +n hh 
Vs — 2V. + V, 
+— at = + O(h) = 0. 


Equation (32) is our algebraic approximation in this case. It reduces to (22) when 
n = 0. If we accept the idea that V“ be continuable to neighboring meshpoints 
with terms through second order in a Taylor Series, as is implied by equations 
(23)-(27), then equation (32) is the unique algebraic representation of our 
problem for this case. Terms through second order are necessary in order that the 
local truncation error vanish as 4 — 0. Some schemes which appear on the surface 
to represent our problem with local truncation error tending to zero when h — 0 
do not do so. Douglas, Garder, and Peaceman treated a moving boundary prob- 
lem. According to Douglas, Garder, and Peaceman [10] to treat the case just 
considered we introduce an additional meshpoint at 0. Then we can easily develop 
a divided difference approximation to Laplace’s equation at a in terms of Vz, Vs, 
Vo, Va, Ve. We obtain 





GS) VotVae [ Zea vs _ Ya— ve) 


h+ n h 
Ve — 2Va + Va 


+o + of) =0. 





However, we now have an extra unknown V» so that we must have an additional 
equation. As our additional equation we choose an algebraic approximation to 
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condition (d) of Problem A, namely 


on yi- V fa D 
(34) ti NE = D,———* — [2+ 2 a-a| Van + O(h*). 
-—417 n 2 2 


The truncation error term which is O(h) is explicitly given in (34). We can compare 
(33), (34) with (32) by using (34) to eliminate V» from (33). Solving (34) for 


Vi =, , , 
— and then subtracting — from both sides, we obtain 
7 u] 


(35) 





Vo-—Ve Vi-Vo s| ot 
fh-—a)+n 2LEh—2) +1 
Substituting (35) into (33), we obtain 
2 [ Ve — Va x V.- ve) 4 Vs — 2V,+ Va 
h+nli(h—a)+ h h? 


+[(Sobteios 
n+i(h—nh+y 


| (h _ n) Ven + O(h*). 


(36) Veo + Viz = 











| Van + O(h) = 0. 





h 
Equation (36) does not agree with (32) unless 7 = i 


or 7 = h. Furthermore, 


in general, the local truncation error does not vanish as 4 — 0. When D,; = Dg, 
it is easily shown that V® and V® are parts of the same harmonic function. 
(This is a theorem of potential theory, cf. Kellogg [11 ]). Since this is so, it would 
be expected that (32) and (36) should reduce to the ordinary five point formula 
for approximating the Laplace equation. (32) does this but (36) does not. The 
local truncation error in (36) remains finite when 4 — 0 even for ¢ = 1. The fact 
that the Jocal truncation error in (36) remains finite as h — 0 does not mean that 
the error in the solution remains finite as 4 — 0. The local truncation error may 
be regarded as a ficticious charge density distribution and the error in the solution 
as the potential of this distribution. As 4 — 0 the local truncation error in (36) 
remains finite. Hence the error in the solution approaches the potential of a 
charged layer whose thickness approaches zero while the charge density in the 
layer remains finite. The potential of such a distribution will become everywhere 
arbitrarily small as the layer thickness goes to zero. 

Equation (34) causes difficulties when the system of algebraic approximations 
is to be solved by an iteration method. The rate of convergence of most methods 
will tend to zero as 7 — 0. Of course, the local truncation error in (36) may be 
reduced by approximating V,,, by a finite difference approximation. 

6. Algebraic approximations with more than five meshpoints. Except for 
special cases as we have discussed, no algebraic approximation with local trunca- 
tion error which vanishes as h — 0 is possible for the points adjacent to I’ in terms 
of five meshpoints. If we make use of one additional neighbor and annex one 
additional equation to (10)—(14), then there is no difficulty in obtaining an alge- 
braic approximation. We solve (10)—(13) together with our extra equation to 
obtain V,, and Vz, in terms of values of V at meshpoints. Then we substitute 
these expressions for V,, and V,, into equation (14). 
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But which extra neighbor do we choose, and why? The problem becomes 
very unsymmetrical. If we make use of several additional neighbors, say the four 
corner points in figure 2, then a symmetrical formulation is possible, but now the 
solution is no longer unique, there being more equations than necessary. 


6 2 5 
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; 
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Ficure 2. 


In figure 2 we show the four neighbors of figure 1 relabeled ‘‘1,’’ ‘‘2,”’ “‘3,”" “4.” 
We show the central meshpoint ‘‘a’’ relabeled ‘‘0,”” and we show the corner points 
5,” “6,” “7,” “8.” Let us suppose that the discontinuity intersects one of the 
line segments 0—1, 0-2, 0-3, 0-4, so that a special representation is necessary. We 
choose an arbitrary origin on the discontinuity and let the coordinates of the 
meshpoints relative to this origin be (¢;, 2;),71 = 0,1, ---, 8. Then the continua- 
tion equations analogous to (10)—(14) are 


(37) Vi-— Vo = (i; — om) Vn + (ts — to) Vi + 3 (22 — 10?) Vann 
+ (Bind; — Pomote) Vaz + $062 — to?)Vu, ¢ = 1,2, 3,4. 


where 5; = 0 when meshpoint 7 is in G1, 6; = 1 when meshpoint 7 is in Gz. 
In addition we have the continuation equations for the corner points, 


(38) Vi-— Vo = (ein; — mo) Vn + (ts — to) Vi + 3 (m2 — 107) Van 
+ (nds — Emoto) Vaz + $(t2 — te?)Viur, * = 5,6, 7, 8. 


The values of V,, and V;; obtained from any solution of (37), (38) are independent 
of where the origin is chosen on the discontinuity. To see this, suppose that the 
origin is displaced on the discontinuity so that t; is replaced by t; + At. Then the 
coefficients of V,, V; and V,, in (37), (38) are unchanged, but the coefficients of 
Vat, Vee are increased by (f'n; — £9) At and (t; — to) At respectively. 

These latter quantities are proportional to the coefficients of V, and V;. When 
we solve equations (37), (38) for Van and Vi, (gn! — Fono)V, and (t; — to) V: 
must be eliminated, and when this is done the extra terms (f'n; — £m) AtV a: 
and (t; — to) AtV«: will also drop out of the equations. 

In order to obtain a non-singular set of 5 equations from (37), (38), we may 
choose the set of four equations in (37) and as a fifth equation choose 


(39) Vs — Vet Vi — Va = aiVn + a2Vi + O3Van + O4Vne + O5Vier, 
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where @, 2, @3, @4, @5 are the coefficients of V,, Vi, Van, Vat, Vig obtained when we 
form equation (39) from (38), by taking the combination of equations (38) 
indicated by the left side of (39). The system of equations (37), (39) is a plausible 
set to use for obtaining an algebraic representation. Loosely speaking, (37) does 


not give a “‘well-set’’ determination of V,, while (39) does. | V,. drops out of 


(37) when a = =.) 


From the system (37), (39) we obtain (21), (22), and (32) in the respective limit- 
ing cases. Arms, Gates, and Zondek [12] have shown that extrapolated line 
relaxaticn may be used to solve difference equations involving nine-point formulas. 
Of course, it is not practical to solve (37), (39) to obtain a explicit formula analo- 
gous to (21), (22), (32) for the general case. However, it is not difficult to use’a 
computer to calculate the coefficients in the algebraic representation for any given 
case. First we compute the matrix of coefficients of equations (37), (39) and then 
invert this matrix, and finally we combine elements of the inverse to obtain the 
values of the coefficients which multiply each V; in the algebraic representation 
of Von _ Viz z= 0. 

Other methods for determining an algebraic approximation may be based on 
an integral formulation of our problem. Referring to figure 2, we let 


& ‘ k 
(40) Ex = f DnV,™dx, En = f DnVe™dy, 
i 3 


i.e., Ey is the line integral of D,,V,™ along the line /k shown in figure 2, Ey the 
line integral from j to k. Superscript m will be one over a segment of the line of 
integration lying in region one, and will be two over a segment in region two. 
Then as may be shown from the conditions of Problem A, we have 


(41) Eu + Ex, ~~ Ea _ Ej; = 0. 


The integrals Ex, Ej, etc. may be approximated by algebraic approximations. 
Again values of V at some of the corner points 5, 6, 7, 8 will be required in making 
the approximations. Algebraic approximations based on (40), (41), have the 
advantage that the resulting equations obey a conservation law which is an 
algebraic analogue of Gauss’s theorem £D,,.VV™dS = 0. In the case of moving 
boundary problems this method also has the advantage of giving a unique 
velocity for the moving boundary. 

7. Curved Interfaces. If T is an arc with continuous curvature, it may be 
approximated in neighborhoods by simple analytic arcs, for example, by osculating 
circles. The mesh spacing must be small enough so that the radial distances from 
the osculating circle to the nine meshpoints to be used in the algebraic approxima- 
tion are less than and are’preferably small compared to the radius of curvature. 
There are two ways in which we may obtain the continuation equations. One 
method consists in mapping the osculating circle onto a segment of the real axis in 
the y-plane. This mapping is easily obtained. The mapping Z = Z» + Re®/Rei(a/®) 
maps the real axis of the y = a + 18 plane for 0 < a < 27 intoa circle with center 
at Zp and radius R in the Z-plane. The inverse of this mapping will map a circle 
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onto the segment (0, 27) of the real axis. Under the mapping the nine meshpoints 
to be used in the algebraic approximation are mapped into nine meshpoints in 
the 7 plane. The coordinates of these meshpoints in the y plane, relative to an 
origin on the real axis, are substituted into the continuation equations (37), (38). 

A second method consists in expanding V“ in a power series about a point O 
on the osculating circle, and choosing ¢ to be arc length from O on the osculating 
circle. We have as before V2? = ¢V{, Vi? = V{?, VIP? = VP, VP = ev? but 
V2) ~ V%). Expressing Laplace’s equation in polar ple whe with origin at 
the center of curvature, we have 


r 1 4 : 
(42) VQ + — Vn + ViP = 0 
whence it follows that 
(43) v2 = yo 4 $F yw, 
rT 


Then the continuation equations (37), (38) become 


wy r= ten [(em eat) (em entt)) 


+ (te = to) Ve + 3 (m2 — 10?) Van + (Einsts — Eorote) Vie 
+ 3(t? — to?) Vier, 1 sad 1, 2, 3, .; 7 8. 





In (42), (43), (44) r is positive if the center of curvature is in G;, negative if 
in G2. This follows because m was the normal from G; into Ge. 

8. Other equations. It can be shown that the power series method here de- 
scribed may be used to obtain algebraic approximations to analogues of Problem 
A for certain other elliptic or parabolic equations, for example, for the Helmholtz 
equation and for the diffusion equation. 
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Simplified Calculation of the Exponential Integral 
By James Miller and R. P. Hurst 


Calculations of the negative exponential integral function [1] 
oo t 
(1) -Fi(-) =f — at 


are usually made using one method for x < 3, another for 3 < x < 20, and a 
third for x > 20. The intermediate region is of considerable interest for many 
applications, and it is here that the function is most difficult to compute. The 
Ruedenberg [2} and Kotani [3] methods are valid for any argument, but diffi- 
culties in application limit their usefulness to the intermediate region where less 
complicated methods fail. 

The authors wish to present a simple method which is useful over the entire 
range of arguments. The method is based on a Taylor’s expansion of the slowly 
varying exponential product function, — e* Ei (— x). An interpolation formula is 
obtained by which this function may be computed from a table value, 
— e* Ei (— xo); then — Ei (— x) is obtained upon dividing by the exponential. 
The accuracy of the result is limited only by the accuracy of the table value and 
that of the computed exponential. The positive exponential integral 


zat 
Ei (x) = at 


may be computed by analogous means. 
The formulae for the negative exponential product function are developed as 
follows : 


Let 
(2) g(x) = — e* Ei (— x), 
then 

g’ (x) = g(x) — x 

g(x) = g(x) +x 

g(x) = g(x) — 2b 

giv(x) = g(x) + 3x. 
One may show by induction that 
(3) g(x) = g@-Y (x) + (— 1)*(m — 1) ke. 
The Taylor’s expansion, then, is 
(4) — e* Ei (— x) = g(xo) + g’ (xo) (x — x0) + (1/2!)g” (xo) (x — x0)* 


+ (1/3 !)g’” (x0) (% — x0)® + --- + (1/m!)g™ (xo) (x — x0)" + Rosi 


Received 20, September 1957. This research was “in¥™ by the Robert A. Welch Founda- 
tion of Houston and the Office of Ordnance Research, U. S. Army. 
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where 


are D gt) (x1) (x ‘ae %o)**1 


~ (n+ 1)! 
with x < x1 < Xo. 
The preceding equations may be combined to produce the simple formula 


(5) ~eEi(-2) =M+ OTM 
where 
(6) To = g (Xo) = — e~ Fi (- Xo) 


and the 7’s are related by the recursion relation 


(7) T, = == #7, , pn &- 2") (n > 1). 


(— xe)* 


For the positive exponential product function we have 


(8) e~* Ei (x) = Uo+ > Un 

where 

(9) Us = e-* Ei (xo) 

and 

(10) B's =—*ly., + ga | (n > 1). 
n (— xo)” 


Equations (5) through (10) are very convenient for computing the negative 
or positive exponential product functions where a table value is known. For ex- 
ample, with a < = convergence to 16 figures is obtained with 12 or 

0 
fewer terms; convergence to 8 figures requires 6 or fewer terms. The spacing of 
the arguments in the accompanying table is such that this condition is fulfilled 
over most of the range 2.5 < x < 80. 

To prevent loss of figures, x — x9 should be negative for computing 
— e* Ei (— x) and positive for computing e~* Ei (x). The accompanying table 
of — e* Ei (— x) was computed starting with the value for x = 80, where 
— e* Ei (— x) was computed using the asymptotic expansion; each computed 
value became the T> for the next lower argument. The table of e~* Ei (x) was 
computed going up from the value of the function for x = 0.2, where Ei (x) was 
computed using the usual MacLaurins series, then multiplied by e~* to obtain 
the product function. All calculations were made on an IBM 650 computer using 
double precision, floating point arithmetic. 
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EXPONENTIAL INTEGRALS AND EXPONENTIAL PRODUCT FUNCTIONS 
x eo” 71 (x) Ei (x) ~e*Ei(-x) ~Ei(-x) 
0.20 -0.6728006649831373( 0) =-0.821760587902400% 0) 061493348746932240(1 0,1222650544183893( 1 
0.25 -0,.4225331193688149( 0) =-0.5425432666619137 0) 061340865644831393(1 0.1044282634643738( 4 
0.30 -0.2242223687152438( 0) -0.302668539265825% 0) 0.3222535605080586(1 0.9056766516758467( 9 
0.35 -0 6630230 7596780976(=1) -0.89434001918464643(-1) 061127045349862318(1)  0,7942154346208358( 0 
0,h0 0.7022622616783961(-1) 041067652186193248( 3 01047828008456006 (1 0.7023801186656624( 0) 
0 0.1816318255906013( 0) 0.26468554053465651( 0) 0.9807147233583903(0 0 .6253313163232692(4 0 
0.50 0.2754982985512703( 0) 0.45462199048631736 0) 0.9229106324837305(0 0.55977359477616084 © 
0.55 0.3569918279209135( 0) %6152906570621864 0) 0.87246573131596157 (0 0.5033660813923938( 0 
0,60 0.4225198103287018( 0} %.7698612899373594 0) 0.8279334352735088 (0 0 .4543795031894021( 9 
0.65 0.4799618229587596( 0, 49193866682454431( 0) 0.7862775692635576 (0 0.4115169759494795( 0 
0.70 0.5268172627521610( 0} %1066907194624291( 1 0.7526780200295871 (0 0.3737688432335091( 0 
0.95 005703036403055990( 0) 01207332816001222( 1 0.7205015065870451(0 0.34034081291123004 9 
0. 0.6054242952633644( 0} 61367396548212326( 1 0.6912453978028315(0 0.3105965785455430( 0 
0.85 0.6350166235265754( 0) 61485714176225254( 1) 06665047874381068(0 0.2840192685924615( 0) 
0.90 0.6597860062607939( 0) 61622811713696867% 1) 0.6399492266392997(0 0.2601839393259996( 0 
0.95 0.6803337743392053( 0) 61759145611867690H 1 0.6173059209526070 (0 0.2387375236537346( O 
1,00 0.6971748832350661( 0) %1895117816355937% 1 0.596347 3623231941 (0 0.2193839343955203( 0 
1,05 0.7107540774523513( 0) 062031087184099664 1) %.5768820704052510(9 0.2018728132201966( 0 
1,10 0.7214575567019491( 0) %.216737827956340% 1) 65587475561 702364(0 0.1859909045360402( 0 
1.15 0.7296223877856724( 0) 0.230628825185562% 1 0.5418049014535823(0 0.1715553536299986( O 
1,20 0.7355440612554694( 0) 0.244209228519265a1 1 0.5259345318947846(0 0.1584084368514626( 0) 
1.25 0.7394826255791207( 9) %2581047974355476{ 1 0.5110328836740476(9 0.1464133725259102( 0 
1,30 007616677163903875( 0) %+2721398880232024 1) 224970097482616443(0 0.1354509578491291¢ 0 
1.35 0.7423027172881568( 0) 0.2863377453173074 1) 0.46837861375242837(0 0.1254168443886644( 0 
1 0.7415682306020866( 0) .3007207464150640 1 0.4712925524860814(9 0.1162193125713579( 0 
- 0073962499463951793( 0) %+3153106049367744 1) 04594675683131685(0 0.10777743986897664 0 
1. 0.7366163509600129( 0) %+330128544912979% 1 0.4482566692915830 (0 0.1000195824066327( 0 
1.55 0.7326703489608585( 0) 0.3451954502797438( 1 0.4376112831118401(0 0.9288210816420927(-1 
1,60 0.72790 15439636994( 0) 0.360531994901946H 1 0.4274879752994281(0 0,8630833369753978(-1 
1.65 0.7224125688366733( 0} 0.3761587568694135{ 1) 0.4176677732671577(0 0.8024762667334322(-1 
1.70 0.7162953753046061( 0) 0.3920963201354 1 0.4086555960005607 (0 0 .7465464440125306 (+1 
1.6 0.7096326000236425{ 0) 0.4083653659076957( 1 0.3998797703779726(0 0,6948868460463879 (-1 
1 0.7024983823417449( 0) 0.24249867557487934( 1) 0.3914916189703562(0 0.6471312936366886(-1 
1,85 0 6949593560157989( 0) 0.441981607980407S 1) 0.3834651071497194(0 0.6029496682837346{=1 
1,90 0.6870754132617678( 0) 0.459371368695358 1) 0.3757765396688868 (2 0.5620437617453485 (-1) 
1.95 0.6789003963937472( 0) 0.477177876458987% 1) 0.3684062987138854 (0 0.5241437956799866 (=1 
} 2.00 0.6704827097900733( 0) 046954234356001890( 1) %.3613286166882226 (0 0,4890051070806112 (+1 
2.05 0.6618658628802201( 0) 0.5141308568693968( 1) 0.3545313797517170(0 0.4564056299118849 (=1 
2.10 0.6530889531759691( 0) 0.5333235359655902( 1) 0.3479959534707185(0 0,4261434150851507 (+1 
2.15 0 .6441870969978426( 0) 2.553025500901557% 1) 0.3417070338905328 (0 0,3980345664801850 (-1 
2.20 ©.6351918146137759( 0) 0.5732614699814381( 1) 0.3356505134529960(0 0.3719113705193257 (=1 
ve 2025 026261313739651271( 0) 0.594056907503733a( 1 0.3296133666814699 (0 0,3476207311944460 (-1 
2.30 06170311019689313( 0) 0.6154380791333508( 1) 0.32418635424648885 (0 0.3250226717162159(=1 
X- 2.35 0.6079136605259716( 0) 0.6374321074059444( 1) 0.3187498781096143(0 0.3039890883162248 (=1 
2 0.5987992978087617( 0) 0.660067027635053% 1) 0.3135020126067469(0 0.2844026093579567 (=1 
or 2 0.5897060737337514( 0) 0.6833718444858624( 1) 0,.3086303187941678 (0 0,2661555839389467 (=1 
2. 0.5806500637222667( 0) 0.7073765894578601( 1) 0.3035258364859861(0)  0,2491691787026973(-1 
2.55 0.5716455429153259( 0) 0.7321123795013333( 1) 0.2987602152032416(0) 0,2332925697125217 (+1 
of 2. 0.5627051529145893( 0) 0.7576114769787598( 1) 0.2941856623768467 (0 0.2185022180408219 (=1 
2.65 0.5538600528724073( 0) 0.7839073511713031( 1) 0.2897348968068591 (0 0.2047012192534298 (-1 
od 2.70 0.5450600565375012( 0) 0.8110347615222037( 1) 0.2854211067605018 (0 0,1918187182658636 (-1 
2B 0.5363737566762730( 0) 0.8390297228026057( 1) 0.2812379121729163(0 0.17978938215454641(-1 
2 0.5277886381277471( 0) 0.8679297723609165( 1) 0,.2771793306851186(0 0, 1685529244521605 (+1 
1 2.85 0.5193111806092019( 0) 0,.8977738397739503( 1) 0,.27323974657133851(0 0,15805367526804 31 (=1. 
. 2.90 0.5109469522665601( 0) 0.9286024186566349( 1) 0,2694138803956038(0 0.1482401922726118 (-1 
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0.5027006948559761 ( 
0.4945764013486612 
0, 48657 73866677791 
0,4787063521935835 
0,.4709654446071541 
0, 463356309587 1836 
0, 455880 1408222987 
0, 4485377247567375 
0.4413294814467530 
0, 43425550 18691712 
0. 427315581991 3683 
0, 4205092538831298 
0. 4138358141250072 
0.4072943497445684 
0, 40088637618910396 
0, 3946027874400052 
0. 3884500187028375 
0. 3824239214001773 
0.3765228510448905 
0. 37074506 78673489 
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A Note on Summation Formulas of Powers 
of Roots 


By Daniel C. Fielder 


Introduction. Through application of a method due to Newton [1], it is 
possible to find solutions for the sum S; of the bth powers of the m roots of f(x), 
where f(x) is an integral rational function of mth degree in x. As anyone who has 
used this method for large m and 6 can attest, the procedure involves tedious and 
repetitive manipulations which increase rapidly in complexity as b is increased. 

It is shown herein that it is possible to reduce the labor of finding such algebraic 
sums by rejecting redundant information in a general formula for the sums and 
using only pertinent terms to predict the numerical and literal coefficients. 
Inductive reasoning is employed to arrive at a simplified method. 

Conventional means for determining sums. Given the expression 


(1) F(x) = aoe® + ayx* + ---Gn_it + an, 


where f(x) is an integral rational function of mth degree in x, and do, ai, etc., are 
the coefficients of the powers of x. The sum of the bth powers of the roots of 
f(x) are 


(2) Sp = c + cP + ---C,", 


where ¢1, C2, etc., are the roots of f(x). Newton's method makes use of the follow- 
ing tabulation to find S;. ' 


la; ae Sido = 0 
2a2 + Sidi + S2ao = 0 
(3) 3a3 + Side + Soai + Sao = 


| 


For the case of m = 6, to get S, it is necessary to go through » operations, each 
successive one becoming increasingly involved. The only respite occurs if one is 
able to find a table of S,’s of desired scope. Unfortunately, available tables seem 
to end with the fourth or fifth summation. A typical listing of sums and a summary 
of Newton’s method is given by Adams [1]. 

Development of a summation formula. A general summation formula is now 
presented. While in the interests of simplicity it is not desirable to use this formula 
directly to calculate S,, it is found that a limited number of terms of this formula 
for a particular S, can be exploited advantageously to yield complete information 
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about the S,. By induction, it can be shown 


n—1 


S.=-— ~ (na) +— +: x (arae-0)| 


2 A=1 


: -Al: 3 = Kg (an00s-v-.)| 


B=1 A=1l 


1 [2 n—3 n—2—C n—1—C-—B 


* ) pm = (aca x0 10»-c-2-4)| 


4 C=1 B= A=1 


re 





(n — 1)>’s n literal terms 


One needs only to consider S, for m equal to 11 or 12, for example, to note that 
the use of equation (4) results in summations of many repetitive combinations of 
a’s. However, it can be seen that S, is an mth order integral symmetric function 
[2] in ai, a2, etc., with the weight of the symmetric function necessarily equal 
to n. For example, if m were nine, one of the literal terms of the coefficients of 
— 1/a° in equation (4) would be a,’a2*a;'. The weight is the sum of the products 
of corresponding subscripts and exponents. In this case, it is 1 K 2+2X 2 
+ 3X 1 = 9. A second necessary condition is that the degree of the product of 
the a’s (exclusive of ao) must equal the power of the corresponding a». For example, 
a;*a;' would not be found in the coefficient of any 1/a*. Thus, the coefficient of 
(— 1)*/ac* contains all the combinations of a’s (exclusive of a9) which simul- 
taneously have weight of m and degree of k. 

Development of a practical summation form. To specify completely the coeffi- 
cient for each (— 1)*/ao* of a particular S,, one must know (a) what the com- 
binations of the a’s are and (b) what the individual numerical coefficients of the 
combinations are. Obviously, equation (4), after required manipulations, yields 
the answer. However, if the tedious approach of equation (4) is to be avoided, the 
information must be obtained in a somewhat different manner from direct ap- 
plication of equation (4). 

The answer to (a) above can be found from a slight rearrangement of a table 
which Euler first developed in connection with enumeration of partitions. The 
table is presented in extended form in Euler’s “Introductio,”’ Lib. 1, Chap. XVI, 
and is reproduced in Chrystal [3]. The number of literal coefficients for each 
(— 1)*/ao* is found directly from the rearranged table which is shown herein as 
Table 1. A brief explanation of the mechanics of constructing the table is given 
at the end of this paper. The literal coefficients are found by systematically setting 
down combinations of the a’s with weight and degree consistent with m and k 
until the number specified by Table 1 is formed. For example, if » = 11 and 
k = 3, Table 1 indicates that there are ten literal coefficients. 
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The answer to (b) above can be found from a consideration of the exponents 
of the a’s. In the coefficient of each (— 1)*/ao', there always appears the combi- 
nation @;*—a,_141 which has the numerical coefficient . The product m X (k — 1)! 
X 1! divided by the product of the factorials of the exponents of a particular 


5S; = 
S: = 
S3 = 
Sy = 
Ss = 
Ss = 


S; = 


Ss = 


Ss = 


Siro _ 


Su = 


TABLE 2. Values of S,, m = 1(1)11 
— a;/do, 
a;*/ae — 2a2/do, 
— a;°/an? + 3a;a2/ae — 3a3/ao, 
a;*/aot — 4a;°a2/ae® + (40:03 + 2a2*)/ae? — 4a4/do, 
> a1°/ao® + 5a1°a2/act saad (Sa;*a3 + 5a,a2") /ao® + (5a1d4 + 5a203)/ac* - 5as/do, 
a3°/ao° — 6a14a2/ao°+ (6a;%a3+ 9a;2a2*) /ap*— (6a :2a4+12a:0203+ -2a2*) /ag* 
+ (60:05 + 60204 + 3a;")/ae — 6a¢/ao, 
— a;"/ao’ + 7a;'a2/ac® — (7ax‘a3 + 14a;°a2*)/ao® 
+ (7a8a44+21aa203+ 70102") /ao! — (70*,5+ 140,0204+ 7072034 70,0";) /ag? 
+ (7a1d5 + 7a2a5 + 7a304)/a0e = 7a7/do, 
ai°/a® — 8a;8a2/ao’ + (8a1°a3 + 20a;‘a2*)/ac® 
— (8a;'a4 + 32a:5a2a3 + 16a;’a*2)/a05 
+ (8a;*a5 + 24a:°a2a, + 12a;7a;? + 240,020; + 2a2*) /ac! 
— (8a;%a5 + 16a;a205 + 160,030, + 8a2*as + 8203") /ac® 
+ (8a;a7 + 8206 + 80305 + 40,7)/ac? — 8as/a0, 
— a;9/ae? + 9ay7a2/ac® — (9ai°as + 27a;5a2*)/a0’ 
+ (9ay5a4 + 45a;*a2a3 + 30a;*a2*)/a0® 
— (9a;‘as5 + 36a;2a2d, + 18a;°a;7 + 54a;7a-7a3 + 9a;a2") /ao® 
+ (9a;%a5 + 27a;a205 + 274;*a304 + 27a;0203" + 9as'a3)/a0 
— (9aya7 + 18a:a206 + 1800305 + 9asa¢ + 9a2*as + 1820904 + 3a;3*)/a0? 
+ (9aiag + 9427 + 9asd5 + Yadds)/ae — Yae/do, 
a}"°/ao'° sae 10a;°a2/a0? + (10a,’a;3 + 35a °a2") /ac® 
— (10a;°a, + 60a;'a2a3 + 50a;‘a2*)/ao’ + (10a;°as + 50a;*a2a, 
+ 25a;*a3" + 100a;%a-7a; +4 25a17a2*)/ao® = (10a;4a¢, a 40a;*a2ds 
+ 40a;°a3a4 + 60a;*a2*a4 + 600;7a3"a2 + 40a2*a3a1 + 2a2°)/ao® 
+ (10a;*a7 + 30a;2a2a¢ + 30a;°a305 + 150:%a, + 300,020; 
+ 60a;420304 + 10as’a, + 15as%a3* + 100103") /ac* 
— (10a;2ag + 20a,a207 + 20a;a06 + 20a;asa5 + 20020305 
+ 10a2a. + 10a2%a¢ + 10a3*a4)/ac® + (10a:a9 + 10a2a3 
+ 10a3a7 + 10asa4 + 5Sa5”)/ac? — 10a10/ao, 
— ay" /ag4 + 11a;°a2/ao" = (11a;°a3 + 44a;a2")/a0? 
+ (11@;"a,4 + 77a;°a2a3 + 77a1°a2*)/ao® = (11a;°as + 66a:°a20,4 
+ 33a;5a?, = 165a;‘a-a3 + 55a;°a2*)/ao? + (11a,5a¢, a+ 55a;‘aeds 
+ 55asa304 + 110a;7a.%a; + 110a;°a2"a,4 + 110a;*a2a;* + 11a,a2°) /ag® 
= (11a;‘a7 + 44a:%a006 + 440:°a305 + 22a;°a, + 660;7a2"a5 
+ 132a;2a2a304 + 44a,a2%a,4 + 66a;a2"a;” + 11la2‘a; a 22a;*a3*)/ao® 
+ (1lasas + 33a;2a2a7 + 33a;%as06 + 33a;%asas5 + 33010275 
+ 664102030; + 33a;a2a2 + 33a103%a4 + 11a2*°as + 33a2*asa4 + 1103°a2)/aot 
— (Alaya9 + 22a a2a3 + 22010307 + 22a :0406 + 11.0; + 11a:*a; 
+ 22acaga¢ + 22a2asas + 1lasae? + 110;7a5)/ac® + (1laiai0 + 1lasds 
+ 11asag + 1lasaz + 11asa6)/ac? — 11a1:/a0. 
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combination of a’s yields the numerical coefficient for that combination. For 
example, in the summation for Si; given in Table 2, 1/ao* has as the a;*—'ay_241 
coefficient a;°a4g which, in turn, has the numerical coefficient 11. The a;*aga¢5 
combination has 2! X 1! X 1! as the product of the factorials of the exponents. 
If 11 X 3! X 1! isdivided by 2! K 1! X 1!, the result is 33 which is the numerical 
coefficient of the a:°a3a. combination. 

Construction of Table 1. Table 1 shows the number of coefficients of — 1/a» 
as the extreme right-hand entry corresponding to a given m. The next entry to 
the left is the number of coefficients for +-1/a,?, etc. The pattern for the number 
of — 1/ao and + 1/a¢? coefficients is obvious. It should also be noted that entries 
to the left of the stepped division remain fixed for 2 and k. The entries in a column 
for k = (n — h), where h = 0, 1, 2, etc., are formed from cumulative sums, 
starting from the right, of the entries in the row for m = h. For example, 
for n = 6, h = 5, and k = 1, the entry is 1. This entry is equal to the first 1 on 
the right in the row for m = 5. For m = 7,h = 5, and k = 2, the entry is 3. This 
entry is equal to (1 + 2) which is the sum of the first two right-hand entries of 
the row for » = 5. A column is formed in this fashion until the sum of the entries 
in the row (in this case, the fifth) is arrived at, after which the entries in the 
column remain fixed at this value. The fixed values are those to the left of the 
stepped division, as explained above. 

Summary. The method of finding S, may be summarized briefly as follows. 
First, obtain the number of coefficients of each (— 1)*/ao* from Table 1 or an 
extension thereof. Second, by some systematic procedure obtain as many different 
combinations of the a’s as there are coefficients specified by Table 1, making 
certain that the weight of each combination equals m and the degree equals k. 
Third, find m X (k — 1) ! divided by the product of the factorials of the exponents 
to obtain the numerical value of a combination. A listing of S, from m = 1 to 
m = 11 is given in Table 2. 

Georgia Institute of Technology 
Atlanta, Georgia 
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An Alternating Direction Method for Solving 
The Biharmonic Equation 


By S. D. Conte and R. T. Dames 


1. Introduction. The biharmonic operator, defined by 


= ot 
AA = rary 
at? axtay? + oy 


plays a role in the classical theory of plates similar to that of the Laplace operator 
in potential theory. In particular, the biharmonic equation, AAW = 0, together 
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with appropriate boundary conditions, describes the small deflections in a thin 
homogeneous plate under various forces exerted on the plate boundary. From a 
computational point of view, direct numerical solutions of such boundary value 
problems are of considerable interest even in those cases where analytic solutions 
in series form are available. 

The use of relaxation methods for the solution of the biharmonic equation has 
been considered by Frankel [3] and more recently by Held, Heller and Lubell [2]. 
It is the purpose of this paper to present an implicit alternating direction iterative 
scheme which appears to be more efficient than any of the known relaxation 
methods. The procedure is an extension of a method developed by Douglas and 
Rachford [1] for Poisson’s equation. 

2. The differential and difference equations. For definiteness we shall consider 
the problem of determining the deflections W in a square homogeneous plate 
supported along its boundaries and bent by moments distributed along two 
opposite edges. The boundary value problem associated with W may be written: 


ow ow ow 














_ — =0 ’ ’ 
@ aaw a +2 tS (x,y) eR 
7) 

(b) W(e,9) = =0 atx =0,1 forO0 <y <1, 

(2) (c) W(x,y) =0 aty=0,1 for0 <x <1, 
2 

(d) a ap aty =0 for 0 <x <1, 
oy? 
ew 

(e) = g(x) aty =1 for 0 <x <1, 
oy" 


where R is the open region (0 < x < 1, 0 < y < 1) and f(x) and g(x) describe 
the bending moments. 
, , ‘ind 1 
We superimpose a square mesh over the region R with grid size h = 7” where 
M is a positive integer, and employ the following notation: 
w(th, jh) = w,,;, 

Views, 5 = Wi41,5 — 204.5 + Wi-r;, 
and 

VWs, 5 = Wiz2,5 — 4Wigs,5 + 6s, 5 — 44-15 + Wia;. 
The usual finite difference approximations to (2) are: 


(a) Vw. j + 2V,7V wi, 5 + V,'w, j = 0 (th, jh) € R’, 


(bt) wm; =0 ' 

(3) big engs- Gamay any nh De 
(c) Wi0 = Wim = 0 (¢ = 1,2,--- M—1), 
(4) wea = —wW-1 + Hf (th) (s = 1,2,--- M—1), 
(e) Wim. = —Wi, mi + A’g(th) (¢@ = 1,2,--- M—1), 


where R’ is the set of lattice points (ih, jh) in R. 
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Since formulation (3) represents a system of (M — 1)* equations in (M — 1)? 
unknowns, the amount of computing time involved in solving (3) will be imprac- 
tical for large M unless one can employ a numerical method which avoids the use 
of (M — 1)*th order systems. Even when the usual point relaxation methods are 
used, they converge very slowly. Thus for the solution of (3) we consider the 
following iteration scheme which involves line rather than point relaxation: 


(a) aft? = uf) — rags Viul? — 27, 4,V20208) — 144,V5us"), 
(b) ft? = uftt) — 7, vtultt + 7, 4,VSul") (ih, jh)  R’, 
(c) us") =w,;=0 
(4) u $s, oo ai, j (¢ = 0, M) 
(d) uff=w;=0 G=0,M) 
ae hf (ih 
) aa to is) (¢ = 1,2,--- M —1). 
Ui M+. = —Uz,u—1 + A’g(th) 


j =1,2,--- M—1), 


1,2, --- M ~ 1), 


Here rat is an iteration parameter which will be chosen to accelerate convergence 
and u? }, for (th, jh) in R’, is an initial approximation to w,, ;. 

Equation 4(a) is implicit in x alone and represents a system of M — 1 equa- 
tions in M — 1 unknowns along a single line in the x-direction. A similar state- 
ment can be made for equation 4(b) and the variable y. Hence, the numerical 
procedure for executing one iteration consists of first solving 4(a) for u{*j*» at all 
points in R’ and subsequently solving 4(b) for u{";*”. At every stage of this itera- 
tion scheme, the matrices involved are all of order M — 1 with at most 5 non-zero 
elements in any row (either on the main diagonal or on two adjacent diagonals). 
Such “‘quidiagonal” systems can be solved efficiently by an extension of an al- 
gorithm for solving tridiagonal matrices due to L. H. Thomas and will be discussed 
later in Section 5. 

3. Convergence of the iterative procedure. To study the convergence of u{") to 
W;, ;, it is convenient to eliminate uth from 4(a) and 4(b). This yields the single 
equation 


ust) = uh) — yas (Veale? + vials? + av2v2ui)) — 2a vive (ust? — uf"). 
Now let the error e{") be defined by 


ef) = us") — W%, j- 


Then e") satisfies the following equation. 
(5) eft) - ef") = Pua (Vieltt + Vie{n*)) oe 2V2V2e5")) 

— ras Viv5 (eft? — ef)), (ih, jh) € R’ 
together with homogeneous boundary conditions corresponding to relations 4(c), 


4(d) and 4(e). 
The eigenfunctions of (5) are of the form 


eS") = A sin (pmih) sin (qnjh). 
Now suppose 


et) =n, 6, ~,¢=1,2,---M—1. 











 =——.|CUo 
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By substituting this expression into (5) we find that 


(1 — an415,?S,7)? 


(6) Aea(anta) = Ons1(Sp* + Se!) + ng 1Sy'S,*’ 





where we have set 


Qa+1 = 167241, 


. prxh 
S, = — 
» = sin 
and 
= sin f2- 
S, = sin 3 


Since the initial error e{?} can be expanded in the form 


M-—1 
ef} = Co, sin (prih) sin (qxjh) 


P.q=l 
we have 
M-1 
ei} = XL Cir, sin (pxih) sin (qxjh) 
P,q=l 
where 


CH, = IT Ap. (ax) CY. 


k=1 


For all positive ax, it is clear from (6) that 0 < Ap, ¢(ax) < 1 for all and g. 
Thus, after iterations each component of the error decreases by a factor of 


Il Apa (ax) 


k=1 


as long as a (and hence 7x) is positive. The minimum value of A,,,(a) is zero and 
occurs when 





2 
«SPS? 
, M(M —1), " : 
This implies, at least theoretically, that after aa” iterations, with appro- 


priate choices for the a, one could make all components of the error vanish and 
thus obtain the exact solution. In practice, however, it is more efficient to reduce 
all components of the error by a fraction of their original value. With this in mind, 
we seek a method for choosing the a; which optimizes the convergence rate of the 
iterative procedure. 

Let 


1 — aS,?S,? y 
1 + aS,*S,* 


Sn o(a)= ( 
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Since S,‘ + S,* > 25S,?S, it follows from (6) that 
Ap. a(a)< Xp, a (a). 


In view of this inequality, any convergence results based on X,,,(a) will also 
apply when 4i,,,(«) is replaced by A,,,(@). In what follows we shall base our con- 
siderations on the more convenient quantity Kg, ¢. 

4. Determination of the iteration parameters. We next consider the problem 
of finding a set of iteration parameters a, (k = 1, 2, --- ¢) such that the maximum 
of the function 

t t (am — S2S2\2 
Gi(SyS,?, ax) = It Xp, e(ax) = tl (=) 


is a minimum for p and g in the range 1, 2, --- M — 1. This problem may be 
resolved efficiently by treating it as a Chebyshev minimax problem and is 
discussed in [4]. Applying the results obtained there to the function G;(S,?S,’, az) 
we find that a recommended set of iteration parameters is 


(7) am =x", k=1,2,---¢# 
for 0<x<1. This gives an upper bound on G,(.S,?S,?, ax) defined by 








ime 2 \ 
(8) GulSPSe, 08) < Pie) = (LZ eT) 
provided that 
4 
- 1 
(9) x1 < S;* and Sun, “1. 


The iteration parameters (7) are not asserted to be the best possible set. A measure 
of the loss of convergence due to not using the set of optimum parameters 
& (k = 1,2, ---,é) is given by 

GUSpPSé ax) 

Gi(S,*S,?, ax) 





(10) 


The formal process of choosing the iteration parameters a; proceeds as follows: 
First one decides on a value of P;(x). When the initial error is to be reduced by a 
given factor, the choice of P,;(x) determines the number of cycles necessary for 
this reduction. Next, x is determined from (8) while the number ¢ of iterations 
per cycle is computed by solving (9) for #; i.e., 

(11) t#>1+ 4 log Si 

log x 
Hence, after one cycle, consisting of ¢ double sweeps, every component of the 
initial error will be reduced by a factor of P;(x). 

The significance of relation (10) is important for the initial selection of P;,(x) 
and hence x. In order to carry out the iteration procedure in the most efficient 
manner, the value of x should be small enough so that this ratio is close to unity. 
Although a decrease in x means an increase in P,(x), it still may be more efficient 
in terms of total sweeps, to increase the number of cycles while at the same time 
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decreasing the number of iterations within a cycle. At any rate, for this scheme to 
give the best results, P,(x) should be such that x < 0.2. For this range 


GulSs'Seh Be) . 
G(S,?S/?, ax) 
As an example of the method, consider the solution of system (4) using a 
20 X 20 grid. Suppose we desire to reduce the initial error by a factor of 10-*. 


We choose x = 0.2 so that P,;(x) ~ 0.01. This means that 3 cycles are necessary 
for the reduction. The number of double sweeps per cycle is given by 


4 log sin — 


pug penta eg 
+n 35 


Hence, it takes three cycles of eight iterations each, using 


_ (0.2) 
+s nog 





k=1,2,---8 


Tk 


to reduce the error by approximately 10-*. This involves 48 single sweeps over 
the grid. 


5. The solution of the quidiagonal systems. The quidiagonal systems arising 
from the iteration procedure (4) may be written in the general form: 
f G1 D Ey )f uy ) { F; ) 
B. C: Dz E: Ue F; 
A; Bs C; D; E; . 
(12) A, B, Cy D, Ey - ° = 


Am-2 Bu-2 Cu-2 Du-» 
L Au-1 Bu. Cui) | Um-1) | Fue J 




















Such systems may be solved by the following elimination scheme: Let 


a, = Ci, 

Ai = D;/w1, Bo = 0, Bu. = 0, 

7 E,/w, yo=0, Yu-1 = Yu-2 = 0, 
and define recursively 


6; = B; — AjBj-2, 


a; = C; — A,yj-2 — 68j-1, for2<j< M—1, 
Bs = (Dj — 8773-1)/0; 
15 = E;/w;. 
Next form 
ho = 0, 
hy = F;/w1, 


hs = (Fj — Ajhj-2 — Sjhj-1)/oj, 2S GS M—1. 
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The values of u are then obtained successively from the formulas 


Umi = hu-1, 
uz = hy — Bilj41 — Vite, J = M— 2, M —3,---1. 


For the matrices arising from equations 4(a) and 4(b) we have 


B; = D; = —4 for all j, 


(13) A; = E;=1 for all j, 


while the diagonal elements depend on the iteration parameters r;. The com- 
ponents F; on the right hand side of the general system (12) will, in our case, 
involve values of u obtained on the previous iteration and may be evaluated in a 
straightforward manner. When relaxing along lines in the y-direction or along 
either of the lines y = 4, y = (M — 1)h, these components will also require values 
of f(x) and g(x). For one complete iteration (i.e., a double sweep over the grid) 
one must solve 2(M — 1) systems of the form (12) with each system of order 
M — 1. Using the elimination scheme described above, this takes approximately 
30(M — 1)? additions, 13(M — 1)? multiplications and 6(M — 1)? divisions. In 
arriving at this estimate, we have made use of relations (13). Regarding both a 
multiplication and a division as equivalent timewise to 5 additions, we may 
estimate the total time T for one iteration as 25(M — 1)? where r is the average 
time in seconds for one multiplication. For a computer capable of executing 4000 
multiplications per second we have approximately 


(14) T = 6.25 X 10-* (M — 1)* seconds. 


Thus, in the example at the end of Section 4, it would take approximately 54 
seconds to reduce the initial error by a factor of 10-*. . 

6. Remarks. The particular fourth order problem considered here can, of 
course, be solved by more direct methods. The second-order Richardson method 
seems to have been successfully used although it has recently been shown by 
E. Windsor [5] that the unextrapolated Richardson Method does not converge. 
Another successful approach has been to reformulate the fourth order equation 
as two coupled second order elliptic equations which can then be solved by 
standard methods. 

On the basis of the theory presented here, the alternating direction method 
seems to be decidedly superior to other proposed relaxation methods, the factor 
of superiority generally increasing as the number of mesh points increases. 

The method proposed here for solving the biharmonic equation can readily 
be adapted to the case of rectangular regions with other boundary conditions. 
However, for arbitrary regions and for general boundary conditions, the proof 
given here does not hold and hence convergence is not assured. 
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TECHNICAL NOTES AND SHORT PAPERS 
On The Numerical Evaluation of Cauchy Principal Values of Integrals 


By I. M. Longman 


Introduction. This note demonstrates a simple method for the numerical 
evaluation of the Cauchy principal value of an integral 


Pf s(e)ax (x real), 


when the integrand f(x) has singularities at one or more points in the finite 
interval [a, 5] of integration. The necessity for such numerical evaluation arises 
sometimes in the solution of problems in applied mathematics. Such integrals can 
occur, for example, in the solution of wave propagation problems. For instances 
of this see Lamb [1], and Pekeris [2]. Sometimes a singularity can be avoided by 
a suitable transformation in the complex plane. The present note presents a 
simple direct method of dealing with such a singularity. 

Formulation of the method. The method is quite general and can be applied 
separately to each singularity of f(x), and so, for simplicity, we will suppose that 
f(x) has only one singularity which, without loss of generality, we may suppose 
to be atx = 0. Suppose then we require to evaluate by some numerical quadrature 
method the integral 


(1) r= P [* f(e)dx 
where the range of integration is made symmetrical about the origin by splitting 


up the integral into two ranges if necessary. We split up the function f(x) into 
its odd and even components 


(2) f(x) = 40f(«) — f(—*x)] + 4%) + f(—x)] = g(x) + A(e), say. 
Then 


(3) I=P [ewe + Pf" i(eyax, 
Now since g(x) is an odd function the first integral vanishes, and so 
(4) I=2 f ” h(x)dx, 

0 


where either there is now no singularity at x = 0, or, if there is, it is not such 
as to make the integral (4) diverge, since we suppose the Cauchy principal value 
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(1) to exist. If (4) still contains a singularity (e.g. of the form log x) the problem 
is no longer that of finding a Cauchy principal value, but the entirely different 
one of the numerical evaluation of an infinite integral. In such a case the singu- 
larity in (4) can usually be removed by a substitution of the form x = y*, for a 
sufficiently large value of m. But, in general, such a singularity arises from the 
presence of a pole on the path of integration and this must be of odd order if the 
Cauchy principal value is to exist. In this case the numerical evaluation has been 
reduced to the determination of (4) which no longer contains a singularity. 
However, in (4), at least in the first instance, the integrand may reduce to the 
indeterminate form © — © or equivalently 0/0, at the lower limit of integration, 
but this is not a serious difficulty. In the first place it may be that lim h(x) is 


known (as in examples 1 and 2), or it may be possible to eliminate the inde- 
terminate form by algebraic manipulation as in example 3. In any case if h(x) 
or at least the part of it which gives the indeterminateness is an analytic function, 
values at and near x = 0 can be obtained by expansion. Another way of over- 
coming the difficulty is the use an integration formula, e.g. Gaussian [3], [4], 
which does not involve the end points of the range of integration. 

Examples. The method is illustrated by three examples. 


(1) To evaluate 


1 
[= pf (e7/x)dx. 
—1 
We have 


ta Pf’ (1/x) cosh x dx +P (1/x) sinh x das. = 2 f (1/x) sinh x dx 
aml -1 0 


which can now be evaluated numerically with no disturbing singularity. 
(2) To evaluate 


«/2 ez 
I=P ———————_ &. 
0 sin xX — cosSx 


Here the singularity is at x = 7/4, and making the substitution x = 7/4 + y 
a little manipulation shows that the integral is reduced to 


a/4 
[= dior f (sinh y)/(sin y)dy. 
0 
(3) To evaluate 





I=P f on dy 
+ —2V1+siny —1 
The singularity is at y = 0, and proceeding on the same lines as above we obtain 
the form 
«/2 d 
r=2f cd 
° Vi —siny + Vi +siny 


which does not contain any singularities. 
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Mersenne Numbers 
By Hans Riesel 


During 1957 the author of this note had the opportunity of running the 
Swedish electronic digital computer BESK in order to examine Mersenne num- 
bers. The intention of the author’s investigation on the BESK was to check some 
known results, and to examine some Mersenne numbers not previously examined. 

Mersenne numbers are numbers M, = 2? — 1, where # is a prime. See [1], 
which contains a more complete list of references. The Mersenne numbers have 
attained interest in connection with digital computers because there is a simple 
test to decide whether they are prime or composite. This is Lucas’s test [1 ]. 
Furthermore the number 2?-'M, is a perfect number, if M, is a prime, and all 
known perfect numbers are of this form. 

In the beginning of 1957 a program for testing the primeness of the Mersenne 
numbers on the BESK was worked out by the author. This program, using 
Lucas’s test, works for all » < 10000. As a test, this program was run for the 
following values of p: 2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 
2203, and 2281. The result was that the numbers M, are prime for these values 
of p, thus confirming known results. 

After these tests, values of > 2300 were to be tested. These values had not 
been tested before. But since the testing of M, for one value of p of this order 
takes several hours on the BESK, a special program for calculating the smallest 
factor of M,, if this factor is <10-2% = 104 85760, was worked out. This special 
program is based on the following well-known theorem : 

All prime factors g of M, (p > 2) are of the form g = 2kp + 1, and of one 
of the two forms g = 8/ + 1. 

The proof of the theorem is quite simple. If ¢g is a factor of M,, 2”? = 1 (mod q). 
Since ~ is a prime, since all numbers for which 2" = 1 (mod g) form a module, 
and since 2? = 1 (mod g), this module consists of all integral multiples of the 
prime p. Now 2! = 1 (mod q) if g is a prime, and hence g — 1 is a multiple of p, 
and in fact an even multiple (since g must be an odd number). This is the first 
part of the theorem: g — 1 = 2kp (k = 1, 2,3, ---). The second part follows 
immediately from the theory of quadratic residues. Since x* = 2 (mod g) has 
the solution x = 2+*»/? (mod g), we see that 2 is a quadratic residue mod gq, 
hence g = + 1 (mod 8).’ 

By the above mentioned special program for small factors of M, the values 
of M, (mod q) for all primes g of the theorem were now calculated. When this 
residue was =0, the factor g was printed out. When no factor g < 10-2” was 
found, the BESK turned to the next value of ». This program was run for all 
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values of » < 10000. The running time on the BESK for one value of p lay 
between 0 and 2 minutes, depending on the size of the smallest factor. These 
smallest factors are given below in a table. The known Mersenne primes are also 
included in the table. Finally, Cole’s factor of M¢z, and Robinson’s factors of 
Mioo and M;,;, though greater than 10-2”, have been printed in the table. On 
checking the values against other sources, a misprint in Archibald’s note [2] 
was detected. This misprint, which concerns the value of the smallest factor of 
Mi3, seems to have come from Kraitchik [3], p. 24 and p. 92. As a further 
check, all factors in our table have been looked up in Lehmer’s Prime Tables, 
except a few, which are too large, and a separate calculation was made, to check 
that they are of the form 2kp + 1. Since, however, there are disturbances in 
digital computers, it is not absolutely sure that all these numbers really are 
factors of the corresponding Mersenne numbers M,. Those primes ~, for which 
no factor <10-2” of M, was found, are, except the known Mersenne primes, 
omitted in the table. 

When this table had been calculated, the BESK examined the omitted values 
of » with Lucas’s test. This was done for all values of » between 2300 and 3300. 
Since the available running time for testing Mersenne numbers on the BESK 
was limited, every value of p was tested only once. The final remainder, see [1 ], 
was printed out in hexadecimal form. On September 8th, 1957, a run indicated 
that the number M3217 is prime. This result was repeated on September 12th. 
All other numbers tested turned out to be composite ; however, this result could 
be false, since the running time is very long. The testing of M3217 took about 
5* 306™ on the BESK, for one run. 

A table of the smallest factor of 2? — 1, p prime follows. Primes 2? — 1 are 
indicated by a dash. 
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Table of Factors of 2? — 1: fin ot 
2 3 5 7 11 13 
jaa sites - ae 23 wt 
17 19 23 29 31 37 
— — 47 233 — 223 
41 43 47 53 59 61 
13367 431 2351 6361 1 79951 —_ 
67 71 73 79 83 89 
1937 07721 2 28479 439 2687 167 —_ 
97 107 109 113 127 131 
11447 — 7459 88807 3391 _ 263 
151 157 163 167 173 179 
18121 8521 33201 1 50287 23 49023 7 30753 359 
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4 Table of Factors of 2” — 1: ow. NE. Continued 
o 181 191 197 211 223 229 
f 43441 383 7487 15193 18287 15 04073 
n 233 239 251 263 277 281 
7 1399 479 503 23671 11 21297 80929 
f 283 311 317 337 353 359 
9623 53 44847 9511 18199 9 31921 719 
. 367 383 397 419 431 443 
“ 12479 14 40847 2383 839 863 887 
“ 449 461 463 487 491 499 
h 12 56303 2767 11113 4871 983 20959 
s, 521 547 557 571 577 587 
ase 5471 3343 5711 3463 5 54129 
es 593 601 607 617 619 641 
0. 1 04369 3607 at 59233 1 10183 35897 
K 643 659 683 701 719 743 
}, 31 89281 1319 1367. 7 96337 1439 1487 
od 757 761 173 811 827 829 
98 15263 4567 68 64241 3 26023 66161 72953 
oe 839 857 859 863 877 881 
26849 6857 72 15601 82 58911 35081 26431 
a 883 907 911 929 937 941 
8831 11 70031 1823 13007 28111 7529 
953 967 977 1009 1013 1019 
3 43081 23209 8 67577 34 54817 6079 2039 
7 1021 1031 1033 1039 1049 1051 
40841 2063 1 96271 50 80711 33569 35 75503 
1091 1093 1097 1103 1117 1129 
87281 43721 9 80719 2207 53617 33871 
1153 1181 1187 1193 1201 1213 
7 2 67497 47 42897 2 56393. 1 21687 57649 3 27511 
1223 1229 1231 1249 1279 1303 
: 2447 36871 5 31793 97423 — 104 44849 
1321 1327 1361 1367 1381 1423 
7927 27 30967 8167 10937 8287 16 99063 
} 1433 1439 1447 1451 1453 1459 
20063 2879 57881 2903 8719 14591 
) 1481 1489 1499 1511 1523 1531 
. 71089 71473 2999 3023 25 22089 88799 
L 1543 1559 1583 1607 1609 1663 
3 1 01839 3119 3167 28927 23 94193 16631 
9 1667 1693 1697 1709 1721 1723 
9 13337 10159 12 35417 379399 1 75543 17231 














210 


1741 
10 02817 


1847 
33247 


1931 
3863 


2003 
4007 


2081 
2 66369 


2179 
2 48407 


2297 
57 65471 


2399 
4799 


2539 
25391 


2617 
78511 


2689 
71 58199 


2767 
6 25343 


2903 
5807 


3001 
32 17073 


3121 
11 23561 


3257 
97711 


3359 
6719 


3499 
35 82977 


3539 
7079 


3623 
7247 


3779 
7559 
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1777 
10663 


1861 
10 23551 


1973 
1 22327 


2011 
21 71881 


2083 
72 69671 


2203 


2333 
3 12623 


2411 
19289 


2543 
5087 


2621 
15727 


2699 
5399 


2819 
5639 


2909 
1 10543 


3011 
20 23393 


3137 
20 01407 


3299 
6599 


3361 
5 57927 


3511 
35111 


3541 
7 64857 


3659 
41 05399 


3793 
60689 


93 


3 


43 


1789 
39359 


1871 
14969 


1987 
67559 


2017 
38711 


2111 
41983 


2207 
23593 


2339 
4679 


2417 
14503 


2549 
63889 


2657 
37199 


2707 
73249 


2837 
22697 


2939 
5879 


3023 
6047 


3163 
55473 


3319 
33191 


3391 
98409 


3517 
62721 


3557 
26841 


3701 
11031 


3797 
91129 


ctor 
1801 
28817 


1877 
15017 


1993 
11959 


2039 
4079 


2113 
22 31329 


2213 
53113 


2351 
4703 


2447 
81 63193 


2591 
7 66937 


2663 
63913 


2711 
5 85577 


2843 
1 42151 


2953 
88591 


3037 
18223 


3181 
19087 


3323 
23 12809 


3433 
44 42303 


3527 
63487 


3571 
64279 


3719 
43 51231 


3803 
7607 


1811 
3623 


1879 
6 05039 


1997 
3 95407 


2063 
4127 


2131 
32 73217 


2251 
4 00679 


2389 
71671 


2459 
4919 


2593 
15559 


2677 
3 64073 


2741 
82231 


2857 
36 79817 


2963 
5927 


3041 
24329 


3217 


3329 
26633 


3457 
20743 


3529 
1 05871 


3581 
21487 


3733 
13 14017 


3823 
89 68759 


1 


63 


18 


47 





1823 
20319 


1913 
01423 


1999 
07097 


2069 
26903 


2141 
89663 


2281 


2393 
33503 


2531 
89673 


2609 
36527 


2687 
98839 


2753 
95727 


2897 
17383 


2999 
71777 


3119 
24953 


3221 
44201 


3347 
26777 


3491 
6983 


3533 
90783 


3593 
21559 


3767 
30137 


3851 
7703 








_—- OS aS Se 








33 


26 


16 


29 


65 


25 


3853 
90641 


3931 
18047 


4057 
97369 


4211 
8423 


4289 
34313 


4457 
31263 


4597 
27583 


4723 
09369 


4861 
29167 


4957 
48711 


5021 
40169 


5167 
31089 


5231 
10463 


5303 
10607 


5431 
54311 


5519 
62281 


5639 
11279 


5791 
26561 


5903 
11807 


6113 
30783 


6199 
61991 
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3863 
7727 


3967 
63473 


4073 
56 69617 
4217 
21 84407 


4297 
8 93777 


4483 
58 54799 


4603 
6 26009 


4759 
1 14217 


4871 
9743 


4967 
3 27823 


5039 
10079 


5171 
10343 


5233 
9 94271 


5347 
3 10127 


5437 
20 87809 


5521 
88337 


5683 
8 75183 


5801 
34807 


5939 
6 17657 


6131 
12263 


6203 
1 48873 


5 


—s 


36 


3907 
31353 


3989 
91473 


4099 
73783 


4229 
29863 


4373 
61223 


4513 
35391 


4639 
94839 


4787 
14889 


4903 
49031 


4993 
79889 


5051 
10103 


5179 
48593 


5237 
22017 


5393 
32359 


5441 
56353 


5531 
40271 


5711 
11423 


5807 
39369 


6011 
88529 


6143 
84593 


6263 
12527 


tor Continued 
3911 3917 
7823 4 07369 
4001 4013 
24007 1 20391 
4127 4129 
74287 15 85537 
4243 4271 
1 01833 8543 
4391 4441 
8783 26647 
4517 4549 
27103 1 36471 
4649 4657 
5 95073 11 45623 
4793 4801 
67103 28807 
4919 4933 
9839 29599 
4999 5003 
2 09959 10007 
5087 5101 
40697 81 10591 
5189 5197 
1 55671 31183 
5261 5279 
42089 10559 
5399 5413 
10799 53 26393 
5449 5483 
25 17439 8 88247 
5557 5563 
33343 5 34049 
5717 5743 
34303 § 43217 
5843 5849 
7 12847 41 05999 
6029 6067 
84407 2 06279 
6163 6173 
5 91649 37039 
6277 6287 
37663 50 29601 


211 


3923 
2 19689 


4019 
8039 


4153 
91367 


4273 
25639 


4447 
71153 


4561 
72977 


4691 
5 72303 


4813 
28879 


4943 
9887 


5011 
80177 


5107 
51071 


5227 
11 29033 


5297 
44 49481 


5417 
4 33361 


5501 
6 93127 


5623 
11 35847 


5779 
9 24641 


5861 
46889 


6101 
79 06897 


6197 
2 97457 


6311 
36 47759 





212 


6317 
6 94871 


6389 
1 91671 


6491 
12983 


6569 
66 21553 


6719 
12 90049 


6883 
18 85943 


6967 
10 72919 


7057 
6 21017 


7207 
1 72969 


7487 
1 79689 


7591 
18 97751 


7691 
15383 


7853 
47119 


; 8017 
1 28273 


8171 
8 00759 


8243 
16487 


8377 
50263 


8461 
14 04527 


8629 
2 58871 


8719 
28 77271 


8803 
28 34567 


8951 
17903 
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19 


48 


80 


~ 


89 


11 


70 


43 


92 


dS 


25 


NS 


6323 
12647 


6397 
44689 


6521 
90751 


6607 
34113 


6737 
41071 


6899 
13799 


6977 
41863 


7079 
14159 


7211 
14423 


7529 
05407 


7639 
50009 


7741 
64079 


7873 
80951 


8059 
43567 


8179 
18513 


8269 
27673 


8387 
76023 


8467 
03209 


8663 
17327 


8731 
14529 


8807 
28983 


8963 
63263 


6343 
90 45119 


6421 
56 89007 


6529 
7 31249 


6637 
52 03409 


6781 
1 08497 


6917 
55337 


6983 
13967 


7103 
14207 


7297 
18 68033 


7537 
23553 


7643 
15287 


7759 
1 39663 


7883 
15767 


8101 
7 12889 


8219 
19 88999 


8273 
10 42399 


8419 
4 88303 


8501 
6 80081 


8669 
26 18039 


8737 
6 81487 


8821 
4 05767 


9001 
15 12169 


~ 


67 


21 


18 


10 


65 


6353 
38119 


6449 
51593 


6551 
13103 


6673 
66423 


6803 
17697 


6947 
80623 


6997 
69071 


7129 
67799 


7411 
22719 


7561 
89759 


7673 
84153 


7789 
72369 


7901 
47407 


8111 
16223 


8221 
20753 


8293 
99033 


8423 
53767 


8513 
72417 


8693 
60791 


8741 
69929 


8839 
23183 


9007 
90071 


6367 
63671 


6473 
27 96337 


6553 
14 54767 


6691 
3 21169 


6841 
41047 


6949 
68 79511 


7039 
12 52943 


7151 
14303 


7417 
1 18673 


7573 
45439 


7681 
39 94121 


7823 


- 15647 


7919 
63353 


8161 
67 08343 


8231 
2 14007 


8317 
3 82583 


8429 
67433 


8597 
4 81433 


8699 
4 17553 


8761 
1 92743 


8867 
70937 


9029 
7 22321 





6373 
38239 


6481 
6 22177 


6563 
13127 


6703 
4 42399 


6871 
2 88583 


6959 
55673 


7043 
14087 


7193 
54 09137 


7459 
1 34263 


7589 
2 88383 


7687 
4 45847 


7841 
23 67983 


8011 
80111 


8167 
31 36129 


8237 
10 54337 


8353 
50119 


8431 
8 09377 


8599 
85991 


8713 
16 03193 


8783 
2 28359 


8929 
1 96439 


9059 
18119 
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Table of Factors of 2” —1:, © —Continued 


actor 

9109 9127 9137 9157 9173 9181 
49 55297 1 46033 10 05071 30 76753 4 95343 13 77151 
9221 9283 9311 9323 9337 9341 
7 19239 59 41121 14 15273 895009 26 14361 74729 
9343 9371 9391 9397 9403 9419 
1 49489 18743 93911 2 25529 25 76423 18839 
9421 9431 9461 9479 9491 9497 
85 73111 6 79033 75689 18959 §=5 31497 5 31833 
9511 9521 9539 9601 9613 9619 
95111 3 61799 19079 22 85039 57679 6 15617 
9791 9811 9829 9833 9851 9859 
19583 77 70313 7 07689 45 62513 78809 11 04209 

9883 9949 9973 

1 58129 80 98487 2 99191 


A Computation of Some Bi-Quadratic Class Numbers 
By Harvey Cohn 


A fascinating chapter in computational number theory began when Lagrange 
showed that every positive integer is representable as the sum of at most 
four perfect squares [1]. Clearly three would not suffice in every case, as 
7 = 2? + 1 + 1? + 1? would be an exception; nevertheless, the problem of ex- 
pressing some positive integer m as the sum of at most three squares soon achieved 
a very special role. For, Gauss showed that r(m), the number of such representa- 
tions, is connected in a very simple way with the much studied (intrinsically 
positive) class number, h, of the field generated by ~—n. Specifically for n 
square-free (and m ~¥ 1, 3 where hk = 1), 


(1) r(n) = gh 


where g = 12 for m = 1, 2, 5, 6 and g = 24 for m = 3 (mod 8). Thus we could 
conclude the existence of at least one such representation for the indicated n. 
Gauss and later, Kronecker, reversed the direction of these equations by making 
large scale tabulations of 4 from r(m), although, unfortunately, no location for 
Kronecker’s alleged tabulation (for odd up to 10,000) seems to exist in the 
literature. In tallying the representation m = x; + x,* + x; it might be noted 
that one must count each ordered triple (x1, x2, x3) of positive, negative, or zero 
integers as a separate unit, so that as much as 2*-3! = 48 could be contributed 
to r(m) when such a decomposition into squares is expressed as triples. 

In more recent times, the representation theory was extended to integers in 
the field k generated by ¥5, i.e., to the quantities 1 = (a + bY5)/2 where a and } 
are of the same parity. Here we seek to represent, necessarily, only those integers 
u which are positive together with their conjugate (i.e., totally positive). Thus, 
e.g., @ >|bV5| > 0. The special surd ¥5 must be used because then, as Gétzky 

Received 8 January 1958. 











214 TECHNICAL NOTES AND SHORT PAPERS 


showed [2], each totally positive integer in k could be expressed as the sum of the 
squares of at most four integers in k. Later, Maass [3] made the more remarkable 
discovery that at most three squares would always suffice ; in fact he arrived at a 
formula analogous to that of Gauss. Since Maass’ formula is the basis of a machine 
calculation, we avoid irrelevant complexities by making certain further assump- 
tions. First of all uz is to be free from square integral divisors in k except for powers 
of (V5 + 3)/2 = [(W5 + 1)/2} = «. Secondly a > 5b > 0,since if 55 > a > 5b, 
we can continually replace u by w/e. Then R(u), the number of representations 
of « as the sum of three squares, is linked to an intrinsically positive quantity H, 
namely the class number of the bi-quadratic field generated by ¥5 and ¥—xz, by 
means of the following formula (which excludes 4 = 1, (5 + ¥5)/2, and 3 
where H = 1): 


(2) H = R(p)/G. 


Here G = 12 when (a, db) $ (1,3), (1,5), (2,4), (5,1), (5, 7), or (6,0) mod 8 
(actually, when »? + uw = 0 mod 4 is unsolvable for 7 in k). Otherwise, G = 120 
except when (a, b) = (1,5), (1,11), (5, 7), (5,9), (6,0), (9, 3), (9, 13), (13, 1), 
(13, 15), or (14, 0) mod 16 (actually, when 7? + » = 0 mod 8 is solvable for 7 in 
k); in these cases, G = 96. 

A tabulation of R and H for 446 selected values of » = [a, b] = (a + bvV5)/2 
was made on the stored program electronic computer, the IBM 650. The values 
of a, 6 were selected with the restrictions 


100 >a > 5d > 0, 


and that » have (except trivially for powers of €) no square divisors in ; or, 
in terms of ordinary integral arithmetic, the condition is that d, the g.c.d. 
of (3[a + 5], 3[a — 5]), be relatively prime to 5(=(¥5)?), and both d and 
(a? — 56*)/4d? be square free. The machine assembled 446 such couples auto- 
matically into the highest four decimal positions of 446 individual ten digit storage 
locations, in lexicographic order. 

The machine next tallied the decompositions » = £;? + &? + £? where 
&; = [a,, b;], a; = 5; mod 2. Here the three couples ; were scanned in lexico- 
graphic order, with the restriction 0 < 5; < 7, while a; > 0 and a’ < 100, (see 
below). Thus, taking all sign possibilities, with 


* = [a? + a? + a? + 5(b2 + 2? + b;*) ]/2, 
b! = + ayb; + Gob2 + yds, 


the couples a’, b’ were constructed and compared with the 446 cases stored in the 
memory. Whenever a matching entry was located the count was augmented and 
accumulated in the last six decimal positions of the memory word. It might be 
appropriate to mention that the IBM 650 has a special “table look-up” operation 
that searches the memory at high speed for the appropriate entry. Without such 
an instruction the search would have had to be programmed with a considerable 
loss of running time. 

In the final phase the ‘‘words”’ a, b, R(u) were unpacked and the congruences 
were examined automatically to calculate H and to produce the output consisting 
of one IBM card per value of u. (See attached table). 
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| Tabulation of R and H for » = (a + 6V5)/2 
a b R Hia 065 R Hia 05 R Hia 0b R H 
2 0 6 11/133 3 240 2/147 9 144 12/59 1 144 12 
4 0 12 11/133 5 288 3/48 2 216 18/59 3 168 14 
5 1 24% 1/34 0 % 8/48 4 168 141/59 5 % 8 
6 0 32 1134 2 9% 8/48 6 144 12|}59 7 % 8 
7 1 #M 2'|34 44 120 1/149 #1 144 12|)59 9 144 12 
9 1 22% 2134 6 72 6|49 #3 20 2/59 11 144 12 
$ 10 2 2 2/135 1 2% 8/|49 5 288 3/60 2 144 12 
11 1 2 2135 3 144 12|49 7 % 8/160 6 192 16 
, 12 0 48 4/135 7 48 4/49 9 120 10/60 8 192 16 
; 12 2 48 4136 2 % 8|50 2 144 12/61 1 288 3 
13 1 9% 1/136 6 % 8|50 6 144 12|61 3 144 12 
14 0 % 1/137 #1 20 2|50 8 2% 8/61 5S 144 12 
142 2% 2137 #3 168 14|51 1 144 12/61 7 20 2 
Ss. 2 2 @is $6 @6 £2184 8 @6 Bia oo OB 
. 15 3 48 4137 7 288 3/51 5 168 14/61 11 % 8 
9 16 2 2 2138 O 38 4|51 7 144 12}62 O 576 6 
| 17 1 48 4138 2 9% 8/51 9 9% 8/|62 2 240 20 
, 17 3 1200 1/138 4 % 8{|52 O 144 12/62 4 216 18 
s 18 2 72 6|38 6 144 12|52 2 144 12|62 6 264 22 
19 1 48 4139 1 72 6|52 6 192 16|62 8 480 5 
19 3 48 4139 3 9% 8&|52 8 120 10| 62 10 192 16 
20 2 48 4139 5 72 6/|52 10 192 16| 62 12 240 20 
211 #1 120 1/139 7 206 8/53 1 360 3/63 1 288 24 
21 3 48 #4140 2 9% 8/53 3 168 14|63 3 192 16 
22 0 192 2140 4 9% 8/153 5 144 12|63 7 240 20 
22 2 #72 #6140 6 144 12153 7 480 5|64 2 120 10 
r, 22 4 72 #6141 #1 #72 «6/153 #9 480 5/64 4 168 14 
1. 23 1 #=72 #6141 #3 288 3/]54 2 120 10| 64 6 144 12 
d 4 2 #72 #6141 #5 120 1/54 4 168 14|64 10 % 8 
" m4 4 48 #41414 7 9 8/154 6 9% 8| 64 12 168 14 
4 25 1 48 4142 © 192 16/54 8 384 4/65 1 144 12 
re 25 3 192 2142 2 168 14/|54 10 144 12|65. 3 480 4 
2% O % 8142 #4 360 3/55 1 144 12|}65 7 192 16 
re 2% 2 48 #44142 6 144 12/]55 3 192 16|65 9 144 12 
nm 2% 4 120 1/43 1 #120 10|)55 7 144 12|65 11 384 4 
27 +1 «#120 10143 3 #192 16|55 9 144 12|65 13 240 2 
~ 7 3 9% $143 #S 144 12/55 11 9 8/|66 O 192 16 
27 5 % 8143 %7 #144 12|56 2 144 12|66 2 144 12 
2 0 9% 8144 #O 144 12|56 4 % 8|66 4 3600 3 
28 2 % 8144 2 2% 8156 6 120 10| 66 6 144 12 
29 1 192 2144 6 % 81|56 10 120 10| 66 8 144 12 
29 3 48 4144 8 72 6/157 1 «216 18) 66 10 192 16 
he 29 5 %72 #6145 1 384 4/57 3 576 6|66 12 240 2 
nd 30 2 2% 8/145 3 2% 8157 5 480 4/67 1 240 20 
be 30 4 9 8145 7 240 2(|57 #%7 192 16|67 3 192 16 
a 30 6 48 4146 O 288 31|57 9 192 16|67 5 216 18 
ch 31 1 48 4146. 2 120 10|57 11 360 3/67 7 168 14 
31 3 #72 6146 #4 72 #6158 O 288 24|67 9 312 26 
ole 31 5 48 4/146 8 288 3|58 2 144 12/67 11 216 18 
32 2 72 6147 1 2% 8(|S58 4 360 3/67 13 240 20 
“es 32 4 120 10|47 3 216 18|58 6 288 24/68 O 288 24 
ng 32 6 120 10147 5 192 16|58 8 192 16|68 2 192 16 
33 1 #144 12/147 7 #120 10|58 10 120 10| 68 6 288 24 
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Tabulation of R and H for u = (a + bvV5)/2—Continued 


R 


264 
192 
240 


192 
192 
288 
192 
240 
192 
144 
144 
192 
192 
144 
144 


H 


22 
16 

2 
12 
20 

5 
16 
12 
12 
16 
16 

6 
16 
20 
18 
16 
10 
12 
12 
14 
12 
30 
28 
16 
18 
16 
18 

7 

5 
16 
24 

3 

6 
16 
12 

2 
14 
12 
12 

4 
12 
16 
16 
24 
16 
20 
16 
12 
12 
16 
16 
12 
12 





J 
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14 


R 


672 
288 
168 
480 
600 
288 
672 
768 
360 
240 
288 
768 
312 
240 
240 
192 
264 
144 
168 
168 
144 
144 
144 
192 
192 
336 
144 
240 
144 
480 
480 
216 
672 
360 
192 
384 
216 
480 
336 
192 
312 
720 
240 
240 
240 
288 
240 
288 
360 
168 
288 
336 
192 


H 
7 





a 


b 


2 
6 
8 
10 
14 
3 
7 
9 
11 
13 
17 


R 


240 
288 
168 
240 
192 
288 
576 
768 


576 


168 
240 


288 
288 


336 
288 
168 
240 
192 


H 


20 
24 
14 
20 
16 
24 

6 

8 
12 
16 

2 

7 
14 
16 
16 

4 
12 
22 
12 

a 
32 
24 
30 
36 
24 
22 
28 
22 
22 
16 
28 
22 
28 
24 
14 

6 

3 
16 
20 

4 

6 
14 
14 
20 

+t 
24 
24 

4 


28 | 


24 
14 
20 
16 





a 


91 
91 
91 
91 
91 
91 
92 
92 
92 
92 
92 
92 
92 
92 
93 
93 
93 
93 
93 
93 
93 
93 
93 
94 
94 
94 
94 
94 
94 
94 
94 
94 
94 
95 
95 
95 
95 
95 
95 
95 
95 
96 
96 
96 
96 
96 
96 
97 
97 
97 
97 
97 
97 


b 


7 
9 
11 
13 
15 
17 
0 
2 
6 
8 
10 
14 
16 
18 
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Tabulation of R and H for « = (a + bV5)/2—Continued 











er 2 wie « a mre *) ora | & 
97 13 600 5/98 6 432 36) 98 18 360 30| 99 13 288 24 
97 15 432 36|98 8 336 28|99 1 192 16|99 15 192 16 
97 17 288 24); 98 10 336 28|99 3 192 16|99 17 240 20 
97 19 360 3/|98 12 80 7|99 5S 288 24;);99 19 288 24 
98 2 288 24|98 14 288 24|99 7 264 22 

98 4 840 7/|98 16 240 26|99 11 192 16 


The computation was monitored for about the first 200 tally operations to 
make sure the score-keeping was correct in all possible cases. The tallying was, 
as before, basically a question of seeing that every permutation and change in 
sign in the triple (£,, £2, £3) counted as a unit. The total running time was roughly 
2.5 hours. One might remark that the human time involved in computing these 
class numbers H from basic algebraic concepts would have to be measured in 
“‘life-times,”” not ‘‘man-hours.”’ 

The computation was completed 18 April 1958 and was sponsored in part 
by the National Science Foundation Grant G-4222. 


Department of Mathematics 
University of Arizona 
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drei quadraten,” Abhandlungen aus dem Mathematischen Seminar Hansischen Universitat, 
v. 14, 1941, p. 185-192. 


Multiplication Time on The IBM 709 
By D. D. Wall 


Average multiply time is useful for roughly estimating problem running time 
for various problems, as well as for roughly comparing different computing 
machines. Determining average multiply time for the 709 is complicated, however, 
due to its zero-skipping feature, and requires an investigation of runs of zeros 
in binary sequences. The particular problem we solve is that of evaluating 
R(n, 1) = total number of runs of length / in all the 2* words of m bits each, and 
S(n,l) = S3-; R(n, x) = number of runs of length > / in the 2" words of n 
bits each. The resulting 709 average multiply time is 193 microseconds fixed 
point, or 170 microseconds normalized floating point, and the purpose of this note 
is to derive these two numbers. 

We make use of a device which we call ‘‘differencing modulo 2,” which obtains 
an ” — 1 bit number from a given m bit number by writing 1 or 0 according as 
the successive bits in the given number exhibit a change or no change. For 
example, each of the (complementary) 8 bit numbers 11010001 and 00101110 
gives the same result 0111001 as its 7 bit difference modulo 2. 


Received 30 April 1958. 
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By differencing modulo 2, each run of length / in a given word of n bits gives 
rise to a run of 0’s of length / — 1 in the m — 1 bit difference, and, conversely, 
each n — 1 bit number with a run of / — 1 0’s may be obtained by differencing 
either of two unique (complementary) m bit numbers with a run of length /. 
Finally, among all the 2* numbers of bits each, the number of runs of 0’s of 
length / is the same as the number of runs of 1’s of length /, and from these 
considerations we conclude that 


R(n, 1) = R(n —1,i-1) = --- = R(n —1+1,1), 


where, again, R(m,/) = number of runs of length / in the 2* words of m bits each. 

We evaluate R(n, 1) by counting the occurrences of an isolated 0 in each of 
the possible positions of an m bit number: 2*-* of these numbers begin 01, 2*-* 
of these numbers begin 101, etc., and 2*~* of these numbers end 10. Since the total 
number of runs is twice the number of runs of 0’s, we have 


R(mn,1) = 2[2*-? + (m — 2)2"-* + 2"-?] = (m + 2)2"-*, n > 1, and R(1, 1) = 2. 


Therefore, R(n, 1) = (n — 1 + 3)2"-*",n > 1, and R(n, n) = 2. 

To find the number of runs of length > / we use the formula for summation by 
parts, + u(x) Av(x) = u(b + 1)0(b + 1) — u(a)v(a) — Sl v(x + 1)Au(x), and 
obtain 


S(n,l) = > R(n, x) = (n — 1+ 2)2"-". 
z=l 


Now, 709 multiplication requires 64 for each 1 bit in the multiplier plus 6x 
for each run of 0’s, plus 6u additional for each of these runs of 0’s of length 7-12, 
plus 64 more for each of these runs of length 13-18, etc., in addition to certain 
time charges independent of the bit structure. The average number of 1’s in an 
n bit number, all m bit numbers equally likely, is 2/2, and the average number of 
runs of 0’s is (” + 1)/4. This last quantity may be found by again differencing 
modulo 2, since the number of runs in an m bit number is 1 plus the number of 1’s 
in the m — 1 bit difference. 

For fixed point multiplication, 124 are required for instruction time, 12 for 
data transfer, 64 are reserved for other purposes, and 34 must be added as an 
average waiting time at the end of the operation in order to total a multiple of 
12 basic cycles. Therefore, with » = 35 and all m bit numbers assumed equally 
likely, average fixed point multiply time is given by T = 33 + 6(35/2) + 6(36/4) 
+d = 192 +d, where d = 6-}-2-*8[.S(35, 7) + S(35, 13) + ---] is the delay 
occasioned by runs of 0’s of length >6. On evaluating this expression we find that 
d < 1 and so obtain the 193, result. Although the 709 disposes of a zero multi- 
plier in a different and quicker manner, this exception does not modify the 193y 
answer because of the assumption that all 2** numbers are equally likely. 

For normalized floating point multiplication, the mantissa has 27 bits of 
which the first is a 1, and 12y are required to add the exponents. The result- 
ing average time is T = 45 + 6(1 + 26/2) + 6(27/4) + d = 169.5 +d, and 
d = 6-4-2-°*[S(26, 7) + S(26, 13) + ---]is <4}, thereby giving the 170y result. 


International Business Machines Corporation 
112 East Post Road 
White Plains, New York 
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Diagonalization of Hermitian Matrices 
By Carl-Erik Fréberg 
As is well known, real symmetric matrices can be diagonalized very efficiently 
by the method of Jacobi, consisting of repeated two-dimensional rotations. A 
similar method can be constructed for Hermitian matrices (cf. ref. 1, p. 160); 


in this case the corresponding transformation matrices are unitary. 
Apart from trivial factors the general two-dimensional unitary matrix has 


the form 
cos ¢ —sin g-e~# 
“4 ee ( i é, ) 
sin g-e” cos ¢ 


A two-dimensional Hermitian matrix H = ( ¢ ' i = is transformed to 
b+ 1% d 
diagonal form by U-'HU = D, where @ and ¢ are defined by: 
tg@ = c/b 
tg 2¢ = 2Vb? + 2/(a — d) 

The matrix elements of D are: 

dy, = 3(a+d+ R) 

du=}a+d-R) R= (@-d'*+4h +e) 

die _ do = 0 


This procedure is now applied to Hermitian matrices of arbitrary order. Let 


a1 bie — iCi2 Dis — tei 
H = | by + ten G22 bes — ices --- 


and define a real matrix A = (ax) in the following way: 


an = bu ti<k 
ai = Gi 


Qk = Ces t>k 


After a two-dimensional rotation (i, 8) we arrive at a new Hermitian matrix H’ 
and a corresponding real matrix A’ = (a’,,). For computation of the elements 
a’, we introduce the following notations (i < k): 


x1 =a,cosg | Vi =axsing | 2 = y3cosé + y,sin@ 
X2=aircosg | y2=arsing | 2 = —yssin@ + y,cosé 
X3 = Gr COS — Y3 1 cos @ — ye sin 0 
Or sing | 2 = y, sin @ + yecosé@ 


ll 
II 


Gk SIN — Z3 


X4 = Qr COS & V4 
Received 10 December 1957. 
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where 
{ tg 0 = axi/an 
tg 2p = 2Varn + ai/ (aii — one) 
Then we have for r <zorr>k: fori <r<k: 


, , r-# 

Ors = X1 + 21 Ar = X%1 — 22 
, , af 

Q' ip = Xo + Be OQ ir = Xo + 2 
, , 

Ark = X3 — 23 Ark = X3 — 


, , 
Qkr = Xa — % Oke = X4 + 23 


and further a’;; = 3 (ai + an +R), a’ = $e +o — R), o'n = on; = 0. 
Here R? = (ai; — ax)? + 4(a?% + a%;). When both r and s differ from i and k 
we have a’-, = ars. The indices 7 and k can be chosen conveniently, e.g. to fulfill 
the condition |a«|+|a.;| = max. 

The whole procedure is now repeated until all eigenvalues have been estab- 
lished with sufficient accuracy in the main diagonal. 

This method has been successfully tested on SMIL, the electronic computer 
of Lund University, up to » = 15. The matrices had all elements in the main 
diagonal = 1, all upper-diagonal elements = 1 — 7 and all lower-diagonal ele- 
ments = 1 + 7. The eigenvalues of these matrices can also be computed directly : 
ue = cot (x(4k + 1)/4n); k = 0,1, --- (2 — 1). An accuracy of about 8 
decimal digits was easily obtained ; further it turned out that the time consumed 
was about 4 times longer than for the corresponding diagonalization of real 
symmetric matrices. 


Department of Numerical Analysis 
Sélvegatan 14 
Lund, Sweden 


¥ 1. A. S. HousEHOLDER, Principles of Numerical Analysis, McGraw-Hill Book Co., Inc., New 
ork, 1953. 
2. E. G. KoGBet.iantz, “Solution of linear equations by diagonalization of coefficients 
— Appl. Math., v. 13, 1955, p. 123-132. 
. zs J. GrREENstTapDT, “A method for finding roots of arbitrary matrices,’”’” MTAC, v. 9, 1955, 
47-52. 


.. - Lorkrn, “Characteristic values of arbitrary matrices,” Quart. Appl. Math., v. 14, 1956, 
p. 267-275. 
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86[C].—NBS Applied Mathematics Series, No. 53, Table of Natural Logarithms 
for Arguments Between Five and Ten to Sixteen Decimal Places, U. S. Govern- 
ment Printing Office, Washington, D. C., 1958, xiii + 506 p., 26 cm. Price 
$4.00. 


This book tabulates log, x, x = 5(.0001)10, 16D, log. x, x = 2(1)10, 40D, and 
log. (1 +x) and —log, (1 — x), x = 10-"(10-*)10-""", _ » = 1(1)13, 25D. In 
addition there are sixteen constants given to 20D. 

With AMS 31 (see MTAC, v. 8, Rev. 1167, 1954, p. 76) these give log, x for 
x = 0(.0001)10, 16D and enough auxiliary material to make the computation of 
any natural logarithm to accuracy compatible with the tables reasonably con- 
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venient. (Actually, either volume suffices without intolerable inconveniences.) To 
this end there is a careful discussion of interpolation methods in an introduction 
written by Arnold N. Lowan. There are no tabulated differences or other aids to 
interpolation ; linear interpolation yields 9D accuracy. 

The volume is a reprint of one volume of the NYWPA “Table of Natural 
Logarithms,” which appeared in four volumes. Two of the volumes (containing 
logarithms of the integers between 1 and 1 00000) evidently will not be reprinted, 
and the third is AMS 31 mentioned above. 

Three trivial misprints and one worse misprint in the earlier printing are noted. 
The trivial misprints all involve arguments. The important misprint corrects the 
value of log, 9.6061 = 2.26239 83133 487638; the originally printed value was 
erroneous in the seventh through the twelfth significant digits. 

The use of these valuable tables since their first printing in 1941 confirms their 
excellence and accuracy. The format is the same as that of AMS 31, which in 
this reviewer's eyes is perfectly adequate. The printing is easily legible, and the 
book is generally up to the high standards of this valuable series published by 
the National Bureau of Standards in scholarship and utility. 


Cc. B. T. 


87(D, E, L, S, X].—Puitre M. Morse & HERMAN FEsHBacu, Methods of Theo- 
retical Physics, Parts I & Il, McGraw-Hill Book Company, Inc., New York, 
1953. 23 cm., xxii X 1978 p. Price $30.00. 


This remarkable two-volume work is already well known to most readers of 
MTAC, but a review seems to be in order as a means of indexing much of its 
content which relates to computation and numerical analysis. 

The volumes are written by physicists for physicists, and they apply analysis 
to problems of physics without attempting to study the complete mathematical 
formulation or background of the methods used. However, the material presented 
certainly is as good a collection as is easily available of partial differential equa- 
tions of mathematical physics. The treatment is largely by separation of variables 
and by transform methods. 

There is no quarrel with the mathematical attainments of these authors, nor 
with their motivation (which, after all, is toward getting the answers to physical 
problems and hence entirely in keeping with the objectives of MTAC). The re- 
viewer, however, has never fully accepted the development of the usual equation 
of the vibrating string through the use of infinitesimals, although it must be 
equivalent in the final analysis to a derivation using tension and inertial forces 
(the last integrated) over all portions of the string. Morse and Feshbach use the 
infinitesimal approach, whereas the integral approach is found, for example, 
in [1]. The reviewer also tends to view with some alarm the use of Dirac delta 
functions, introduced in’ Morse and Feshbach (p. 123) with a warning that 
“differentiation should be attempted only with considerable caution.” This is 
still a minor criticism ; the function gives valuable hints when properly used, the 
authors seem not to misuse it, and the argument here may be almost pedantic. 

A more constructive approach here would be to turn to parts of the book 
which are likely to be of help to readers of MTAC in connection with numerical 
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analysis and computation. These include methods of solution, discussion of some 
of the more important special functions, and some tables. 

Chapter 4 contains sections on gamma and elliptic functions, asymptotic 
series and the method of steepest descent, conformal mapping including the 
Schwarts-Christoffel transformation, and Fourier transforms (along with some 
mention of Laplace and Mellin transforms). There is a table of properties of and 
relations between gamma functions, elliptic functions and theta functions. 

Chapter 5 contains a table of separable coordinates in three dimensions and a 
general table of second order differential equations and their solutions. 

Chapter 6 contains a table of useful eigenfunctions and their properties 
(polynomials of Gegenbauer, Laguerre and Hermite). 

Chapter 7 contains a short table of Green’s functions. 

Chapter 8 contains a section on applications of transforms to the solution of 
integral equations and a short table of integral equations and their solutions. 

Chapter 9 is devoted to approximate methods—perturbation and variational 
methods. It includes a tabulation of methods. 

Chapter 10, devoted to solutions of Laplace’s and Poisson’s equations, lists 
information about and expansions of trigonometric and hyperbolic functions, 
Bessel functions and Legendre functions. 

Chapter 11, on wave equations, takes up cylindrical Bessel functions, Weber 
functions, Mathieu functions, spherical Bessel functions. It also includes a three 
page table of Laplace transforms. 

Chapter 12, devoted to diffusion and wave mechanics, gets to Jacobi poly- 
nomials and semi-cylindrical functions. 

Chapter 13 includes spherical vector harmonics. 

Finally there is a set of tables. These include the following: for x = 0(.2)8, 
sin x, cos x, tanh x, and e~* all 4D, and tan x, sinh x, cosh x, and eé? all 5S; for 
x = 0(.05)2, sin rx, cos rx, tanh rx, and e~**, all 4D, and tan rx, sinh rx, cosh rx, 
and e*? all 5S; the functions a, 6, x and |¢|, 4D, and ¢ (in degrees) 2D, where 
tanh [x(a — 18) ] = 0 — ix =|¢|e~*, tanh ra = 0(.05).95, 8 = 0(.05).45(.025).5; 
the inverse of this function, giving aand 8, 4D, for @ = 0(.2)4and x = 0(.2)1.6(.4)2; 
for x = 0(.1)4(.2)12, in x, sinh x and cosh x (where real), 4D; spherical har- 
monics, P,*(x),a = 0(1)3, 8 = O(1)a, x = 0(.05)1, 4D or 5S; Legendre functions, 
P.P(x) and Q.P(x), a = 0(1)2, 6 = O(1)a, and (a,8) = (+4, 0), (44, 1), (—4, 2), 
(3, 0), x = 1(.2)3(.5)8, 4D or 5S, P.® (ix) and Q,° (ix), x = 0(.2)2(.5)7.5, 5D or 4 
to 5S; Bessel functions J.(x), Na(x), I(x) = i-*Ja(ix), ja(x) = V4/2x Ja+4(x) 
and a(x) = Vx/2x Nos4(x), x = O(.1).2(.2)8; (Here Na(x) denotes a Bessel 
function of the second kind and is equivalent to Y.(x) of Fletcher, Miller and 
Rosenhead [2]); Legendre functions P, (cos v), a = —1(1)9, 6 = 0(5)90°, 4D; 
auxiliary functions C,(x), Ca’(x),5a(x) and 6,’(x) such that J..(z) = C.(z) sin [6.(z) ], 
dJ4(z)/dz = — C,'(z) sin [6.’(z) ], Na(z) = —Ca(z) cos [62(2) ], dNa(z)/de=C,' (z) 
cos [6.’(z)], 6 and 8’ in degrees, a = 0, 1 and x = 0(.1)1(.2)5, a = 2, 3 and 
x = .1(.1).2(.2)5,a = 4,5andx = .4(.2)5,a = 6,7 andx = 1(.2)5,a = 8,9 and 
x = 1.6(.2)5, 4D or 5S and angles 2D; auxiliary functions D,(x), Da’ (x), 5.(x) 
and 6,’(x) related analogously to j.(x) and mq(x), but this 6,, 6,’ is different 
from that of the preceding table, a = 0, 1 and x = .1(.1)1(.2)5, a = 2, 3 and 
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x = .1(.1).2(.2)5, a = 4, 5 and x = .6(.2)5, a = 6, 7 and x = 1(.2)5, a = 8, 9 
and x = 2(.2)5, amplitudes 4D or 5S andangles (in degrees) 2D; periodic Mathieu 
functions, even and odd, Se. (h, cos x), a = 0(1)4, and So, (h, cos x),a = 1(1)4, 
h? = 0(1)9, x = 0(10)90°, 4D; quantities associated with periodic Mathieu 
functions, M,*(h) = Jo** [Sea (h, cos x) Pdx, Jéa, Ca*, 52° such that Je, + iNe, 
= —C,*ie*=" and analogous functions for odd functions, M,°, Joa’, C.°, 52° with 
Joa’ + Noa’ = iC.°e*@’, h? = 0(1)9, a = 1(1)4 and a = 0 for even functions, 4 in 
degrees 2D, other quantities mostly 4D. 

On the whole, this is an impressive collection of information, both in volume 
and in range of topics covered, and it represents an enormous effort on the part 
of its highly competent authors. References are adequate, although no attempt 
was made to include exhaustive bibliographies. Misprints detected by the rather 
cursory reading of the reviewer were trivial. Numerous and well chosen problems 
are included with each chapter. Reference is facilitated by summaries at the ends 
of the chapters, but the broad scope of the book precludes complete reference 
details in the fields covered. While the work is intended as a text book for a course 
the authors have presented at the Massachusetts Institute of Technology for 
several years, many people will find that it is more useful as a guide to be used 
in their courses rather than a text book to be followed faithfully; in this role of 
guide the book gives a valuable list of equations, solutions and significance. The 
nature of the information available in the book is described in the words of the 
authors as follows: ‘‘The present treatise, therefore, is primarily concerned with 
an exposition of the mathematical tools which have proved most useful in the 
study of the many field constructs in physics, together with a series of examples, 
showing how the tools are used to solve various physical problems. Only enough 
of the underlying physics is given to make the examples understandable. 

“This is not to say that the work is a text on mathematics, however. The 
physicist, using mathematics as a tool, can also use his physical knowledge to 
supplement equations in a way in which pure mathematicians dare not (and 
should not) proceed. He can freely use the construct of the point charge, for 
example; the mathematician must struggle to clarify the analytic vagaries of the 
Dirac delta function. The physicist often starts with the solution of the partial 
differential equation already described and measured; the mathematician often 
must develop a very specialized network of theorems and lemmas to show exactly 
when a given equation has a unique solution. The derivations given in the present 
work will, we hope, be understandable and satisfactory to physicists and en- 
gineers, for whom the work is written; they will not often seem rigorous to the 
mathematician.” 

The Table of Contents follows: 


Vol. Chap. 
I 1. Types of Fields 
1.1 Scalar Fields 
1.2 Vector Fields 
1.3 Curvilinear Coordinates 
1.4 The Differential Operator 
1.5 Vector and Tensor Formalism 








II 
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1.6 Dyadics and Other Vector Operators 
1.7. The Lorentz Transformation, Four-Vectors, Spinors 
2.1 The Flexible String 
2.2 Waves in an Elastic Medium 
2.3 Motion of Fluids 
2.4 Diffusion and Other Percolative Fluid Motion 
2.5 The Electromagnetic Field 
2.6 Quantum Mechanics 
Fieids and the Variational Principle 
3.1 The Variational Integral and the Euler Equations 
3.2 Hamilton’s Principle and Classical Dynamics 
3.3 Scalar Fields 
3.4 Vector Fields 
Functions of a Complex Variable 
4.1 Complex Numbers and Variables 
4.2 Analytic Functions 
4.3 Derivatives of Analytic Functions, Taylor and Laurent Series 
44 Miultivalues Functions 
4.5 Calculus of Residues; Gamma and Elliptic Functions 
4.6 Asymptotic Series: Method of Steepest Descent 
4.7 Conformal Mapping 
4.8 Fourier Integrals 
Ordinary Differential Equations 
5.1 Deparable Coordinates 
5.2 General Properties, Series Solutions 
5.3 Integral Representations 
Boundary Conditions and Eigenfunctions 
6.1 Types of Equations and of Boundary Conditions 
6.2 Difference Equations and Boundary Conditions 
6.3 Ejigenfunctions and Their Use 
Green’s Functions 
7.1 Source Points and Boundary Points 
7.2 Green’s Functions for Steady Waves 
7.3 Green’s Function for the Scalar Wave Equation 
7.4 Green’s Function for Diffusion 
7.5 Green’s Function in Abstract Vector Form 
Integral Equations 
8.1 Integral Equations of Physics, Their Classification 
8.2 General Properties of Integral Equations 
8.3 Solution of Fredholm Equations of the First Kind 
8.4 Solution of Integral Equations of the Second Kind 
8.5 Fourier Transforms and Integral Equations 
Approximate Methods 
9.1 Perturbation Methods 
9.2 Boundary Perturbations 
9.3 Perturbation Methods for Scattering and Diffraction 
9.4 Variational Methods 


ead 
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10. Solutions of Laplace’s and Poisson’s Equations 
10.1 Solutions in Two Dimensions 
10.2 Complex Variables and the Two-dimensional Laplace Equa- 
tion 
10.3 Solution for Three Dimensions 
11. The Wave Equation 
11.1 Wave Motion on One Space Dimension 
11.2 Waves in Two Dimensions 
11.3 Waves in Three Space Dimensions 
11.4 Integral and Variational Techniques 
12. Diffusion, Wave Mechanics 
12.1 Solutions of the Diffusion Equation 
12.2 Distribution Functions for Diffusion Problems 
12.3 Solutions of Schroedinger’s Equation 
13. Vector Fields 
13.1 Vector Boundary Conditions, Eigenfunctions and Green’s 
Functions 
13.2 Static and Steady-state Solutions 
13.3 Vector Wave Solutions 
SS: 3: 


- 1. S. L. Sonatev, Uravneniia matematicheskoi fiziki (Equations of Mathematical Physics), 
os) 7 q 

2. A. FLeTcHER, J. C. P. Mrc_er, & L. ROSENHEAD, An Index of Mathematical Tables, Scientific 
Computes ie 23 Bedford Square, London, W. C. 1, 1946. [See MTAC, v. 2, 1946-47, RMT 
233, p. 13-18. 


88[F ].—R. J. Porter, Irregular Negative Determinants of Exponent 3n with their 
critical classes. Part III, from — D = 1 00000 to 1 50000. 215 typewritten 
pages, 25.5 X 10 cm., deposited in UMT file. 


This table of numbers occurring in the theory of quadratic forms completes 
the listing of “critical classes’? to — D = 1 50000. The list of determinants is in 
three parts [UMT 155, MTAC, v. 7, p. 34; UMT 185, MTAC, v. 8, p. 96; 
Review 62, MTAC, v. 9, p. 126]. The determinants are serially numbered and 
listed together with their critical classes in the present table [Review 84, MTAC, 
v. 9, p. 126; Review 113, MTAC, v. 11, p. 275]. The lists were compiled with the 
aid of a UMT list of groups and series [Review 3, MTAC, v. 9, p. 26] by the 
same author. The 58 values of — D < 1 50000 which have exponent of irregu- 
larity 9 are listed and discussed in [1]. The 14 values of — D < 1 50000 which 
have exponent of irregularity 6 are 


17561 70244 91299 1 23539 1 28739 
55555 80755 1 12723 1 25443 1 29355 
67899 91083 1 19195 1 26755 


The other 5764 values of — D < 1 50000 have exponent of irregularity 3. 


J. L. SELFRIDGE 
University of California 
Los Angeles, California 


1. R. J. Porter, “On irregular negative determinant of exponent 9n,” MTAC, v. 10, 1956, 
p. 22-25. 
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89[F ].—C. L. Baker & F. J. GRUENBERGER, Primes in the Thousandth Million, 
1958, 3 p. 83” K 11” + 3917” X 22” blue line ozalid prints deposited in the 
UMT file. 


A list of 47,957 primes found between, and including, 999,000,011 and 
999,999,937 together with the first eighteen 10-digit primes. 
This table came into being as a by-product of another table of primes: the 
5,761,456 eight-digit primes, which will eventually be published. 
A description of the method of construction of the tables on the IBM model 
704 computer is included. 
Ge & Ft. 


90[F ].—O. P. Gupta, ‘‘Partitions into exactly k distinct primes,”” Panjab Univ., 
Res. Bull., No. 107, p. 283-290, May 1957, Hoshiarpur, India. 


The function tabulated is R,(m) the number of partitions of into precisely 
k distinct primes > 1, for all possible values of k and for m = 1(1)300, together 
with the total number R(m) of partitions of into distinct primes. The function 
R,(n) was determined recursively by means of the auxiliary function R;(m, p;) 
defined as the number of partitions of m into k distinct primes of which the 
smallest is p;, so that 

Ri(n, Pp) “ a Ri-i(n ahs Pi, p) 
P>P1 
and 
Ri(n) = L R.(n, p). 
p22 


Sample values are 
R.2(300) = 21, Rs(300) = 13504, R13(300) = 87, R(300) = 53040. 


D. H. LEHMER 
Department of Mathematics 
University of California 
Berkeley, California 


01[F ].—M. S. Cueema & H. Gupta, Tables of Partitions of Gaussian Integers. 
Natl. Inst. of Sciences, India, Mathematical Tables, v. 1, xii + 67 p., 1956. 
30 X 24 cm. Price 15 Rps. Paper Bound. 


The number B(n, m) of partitions of the Gaussian integer m + im into non- 
zero Gaussian integers of the first quadrant is tabulated for m, m = 0(1)50 in the 
first ten pages of the volume, even though B(n,m) = B(m,n). In other words 
B(n, m) is the number of partitions of the bipartite number (, m) into bipartite 
summands (r,s) where 0 < 7, 0 < s,r +5 > 0. On p. 65 the number of unre- 
stricted partitions 

b(n) = B(n, 0) 


is tabulated along with its sum function 


are x p(k) 
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for m = 0(1)50. The rest of the volume is devoted to the auxiliary function 
B,(n, m) generated by 


Ill 
— 


( XL Bn, m)xry") II If (@ — xy) IT (1 — 2*) 
n,m=0 r=0 t(>k s=1 

which was used recursively to compute B(n,m) = B,(n,m). The printing is 
beautifully done. 


D. H. LEHMER 
Department of Mathematics 
University of California 
Berkeley, California 


02[F ].—Harvey Coun, “Some bi-quadratic class numbers,” [MTAC, this 
issue, p. 213-217.] 


This is a tabulation of 446 sets (a, b, R, H) as is described in the attached note. 


93[F, K].—Ricuarp T. Burcu, Approximate Values of Stirling Numbers of the 
Second Kind for the First Hundred Degrees, Department of Defense, 134 p., 
multigraphed, 27 cm. Deposited in the UMT File. 


The Stirling Numbers, u(N, A), of the second kind are defined by the poly- 
nomial identity 


N A-1 
x% = > u(N,A),x™, where x = JJ (x — &). 
a=i ko 


The recursion u(N + 1, A) = AU(N, A) + u(N,A — 1) is used to calculate 
the numbers to four significant figures for N (the degree) :1(1)100 and A:1(1)N. 


A sum U(N) = > U(N,k) and a sum check U’ recursively computed from 
k= 


N 
U(N +1) =1+> 7) U(j) are included. Values are expressed in ‘floating’ 
j=1 


form: 1344 + 12 standing for 1344-10". The number u(JN, A) is the number of 
(ordered) selections of NV items from A categories, at least one item being selected 
from each category. 


J. D. Swirt 
University of California 
Los Angeles, California 


94[G].—K. Yamamoto, “Structure Polynomial of Latin Rectangles and Its 
Application to a Combinatorial Problem,’’ Memoirs of the Faculty of Science, 
Kyushu Univ., S. A., v. 10, No. 1, 13 p., 1956. Fukuoka, Japan. 


The problem mentioned in the title is due to J. Touchard and asks for the 
number NV, of permutations on 


1,2,3,---,m 
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which are discordant with the first k cyclic permutations 


(1, 2,3, ---, m) 
(2, 3, 4, a... 1) 


ee &- ats 


This generalizes the famous probléme des rencontres (k = 1) and the probléme 
des ménages (k = 2). There is given (p. 13) a table of N, in the case k = 3 for 
nm = 3(1)20. Two auxiliary functions 


N,* = Va + 2, N,! = {Nasi + (—1)"N*,41} /n 
where 
Vero +h (0+5=1,0b = —1) 


in Lucas’ companion to the Fibonacci series, are given for m = 3(1)20. 


D. H. LEHMER 
Department of Mathematics 
University of California 
Berkeley, California 


95[G, P ].—Lorenzo LuNELLI, ‘“‘Numerazione delle maglie in una rete completa,” 
Instituto Lombardo di Scienze e Lettere, Rendiconti, v. 91, 1957, p. 903-911. 


[G, P].—EMANUELE Brionp1, ‘“‘Numerazione delle maglie in una rata 
qualsiasi,’’ Instituto Lombardo di Scienze e Lettere, Rendiconti, v. 91, 1957, 
p. 912-926. : 


[X, Z].—Lorenzo LuneEtut, ‘“Determinazione delle maglie in une rete 
mediante una calcolatrice elettronica,” Instituto Lombardo di Scienze e 
Lettere, Rendiconti, v. 91, 1957, p. 927-935. 


These papers are concerned with the calculation of the number of polygonal 
subnets contained in networks of lines connecting nodes (as in electrical networks). 
In particular, the first paper derives a formula for the number of polygonal subnets 
contained in a ‘‘complete net” (wherein each pair of nodes is connected by one 
and only one line). The formula, for ” nodes, is 


My(n) = EA V'("). 


ae 


The derivation is straightforward. 

The second paper derives in a somewhat different fashion a formula for the 
number of polygonal subnets contained in a general net. The procedure described 
considers the ‘‘complementary net” (the additional lines required to produce 
from the given net a complete net with the same number of nodes). Each set of 
¢ lines contained in the complementary net is then classified into one of three 
types: 1) they form a closed polygon of ¢ lines. (In this case there is only one 
polygonal subnet which passes through all ¢ lines and the notation is used 
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M,°(n) = 1); 2) they form a number, y, of separate connected sequences of lines 
(excluding those of the first case) such that at each node at most two lines meet. 
In these situations, the number of polygonal subnets passing through all ¢ 
lines is given by 


Morin) = © 28(*7 77%) @—o-a), 14423; 
q=7+¢ gore 

3) all other situations (involving more than two lines meeting at some node). 

These must have zero polygonal subnets passing through all ¢ lines. 

Let the complementary net of the given net consist of L lines. Let it have, for 
each ¢, m,’ sets of ¢ lines forming y separate connected sequences. Finally, let 
M,(n) be the number of polygonal subnets in the complete net of nodes. Then, 
the number of polygonal subnets in the given net, as derived in the paper, is 
given by the formula 


L i 
M = M,(n) + & (-1)' LD mfMf(n). 


j=l i=—0 


The paper then derives the formula for the number of polygonal subnets in a 
complete net by induction on n. Specifically, by the above formula, 


M,(n — 1) = M.(n) — (m — 1)M,%(n) + & ‘ *) asc. 
It is then shown that M,(n) = M,.(m — 1) + 3(m — 1)M,°(m). The formula for 
M,(n), as given in the first paper, is then derived by induction on n. 

The third paper is concerned with the use of the CRC-102A at the Polytechnic 
Institute of Milan for determination of the number of polygonal subnets in a 
network. The procedure uses a representation of the set of lines in a given subnet 
as a vector with ones and zeros for its components—a one for a given component 
if the corresponding line in the original net is part of the subnet and a zero if 
the corresponding line is not part of the subnet. The original net itself is then 
represented by a matrix (the Poincaré matrix) with a row for each node showing 
by such a vector the subnet of lines meeting at the node. The columns of this 
matrix then indicate the endpoint nodes for each line. If the Poincaré matrix is 
multiplied onto the vector representing a given subnet, a product vector is 
obtained. The given subnet then is a polygonal subnet if and only if the product 
vector consists entirely of twos and zeros. The computer program examines in 
this manner each of the possible subnets of the given net and tabulates the po- 
lygonal subnets. The program includes an output conversion routine to represent 
the make-up of each polygonal subnet in convenient form for the human operator. 
The program is deliberately unsophisticated, and therefore uses no criteria other 
than the make-up of the product vector for determining the polygonal subnets. 


R. M. Hayes 


Business Systems Department 
Magnavox Research Laboratory 
Los Angeles, California 
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96[I, X, Z].—NatTiIonaL PuysicaL Lasporatory, Modern Computing Methods, 
Notes on Applied Science No. 16, Published by National Physical Laboratory, 
Obtainable from Her Majesty’s Stationery Office, London, England, 1957, 
vi + 128 p., 24 cm. Price 10s 6d. net. 


This is a pamphlet which recounts the experience of professional computers 
attacking extensive problems on an automatic digital computer of considerable 
power. It is scholarly in that full use was made of literature known to these 
competent people, but it is not pedantic and it does not contribute particularly 
to the advanced research literature on numerical analysis. 

However, the usefulness of this book to people who are meeting extensive 
computations for the first time is great, particularly because the larger and more 
detailed books in numerical analysis seem to be only mildly influenced by the 
development of large machines. The reviewer is consistent in recommending it to 
each person starting to work in his own group, and any careful reader of the 
pamphlet is likely to profit. So far as the reviewer knows the material in the book 
is not conveniently available in any other place, and considering the size of the 
pamphlet the choice of material presented is excellent. 

Despite this approbation, it should be noted that this material is dated—many 
contributions to extensive computation have been made since the preparation of 
this material. Thus it should be augmented if it is to be used as a guide to a course 
in modern computation. Problems for students would also be required if it were 
to be used as a text. Numerous examples are worked out in the text, however, 
and these illustrate the methods and make reading easier. Occasionally the 
authors or compositors have been careless in displaying the results of computation. 
For example the publisher has noted that in Table 2 on page 49 a proper conclu- 
sion is drawn from incorrectly copied inputs to the table—a type of computing 
error which must be unusual among the distinguished group contributing to this 
book. On page 99 the right column of numbers in an example illustrating Euler’s 
transformation is spaced artistically but not conveniently for following the 
example (which is easy to follow anyhow). Those of us who contend that time 
and effort spent on care and neatness in computing (and, from the reviewer’s 
point of view, even more importantly in coding for automatic machines) are 
particularly good investments hope that these slight faults will be eliminated in 
later editions. 

The table of contents follow: 


Chapters 


Linear Equations and Matrices (1) 

Linear Equations and Matrices (2) 

Roots of Polynomial Equations 

Latent Roots of Matrices 

Finite-difference Methods 

Ordinary Differential Equations (1) 

Ordinary Differential Equations (2) 

Hyperbolic Partial Differential Equations 

Elliptic and Parabolic Partial Differential Equations 


SPs? FPP PS 
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10. Relaxation Methods 
11. Tabulation of Mathematical Functions 
12. Computation of Mathematical Functions 

ie | 


97[K ].—W. J. Dixon & F. J. Massey, JR., Introduction to Statistical Analysis, 
McGraw-Hill Book Co., New York, 2nd ed., 1957, xii + 488 p., 24 cm. 
Price $6.00. 


The first edition of this excellent textbook contained an unusually extensive 
and useful set of tables. An important feature of the second edition is the addition 
of 21 more tables. Besides this, 12 of the tables of the first edition have been ex- 
tended in the second edition ; five have been modified, one was reduced and one 
omitted. The present review will list the new tables and indicate the other changes. 

The new tables are: 


A-4. Cumulative areas under the normal frequency curve to 4D for the 
standardized deviate = —3.25(.05)3.25. 


A-8b (1) gives the .1, .5, 1, 2.5, 5, 10, 90, 95, 97.5, 99, 99.5, 99.9 percentiles 
of the standardized range to 2D for samples of N = 2(1)20 from normal. 


A-8b (4). The variance and the efficiency to 3D or 3S of four estimates of the 
mean, the median, the midrange, the average of the best two ordered observa- 
tions, and the average of all but two extreme observations in samples of 
N = 2(1)20 from normal. 


A-8b (5). Standardized expected values of the order statistics to 3D in 
samples of N = 2(1)20 from normal. 


A-8b (6). The best linear estimates of o and their efficiencies to 3D for samples 
N = 2(1)10 from normal. 


A-10b. Percentile values covering the range .005 to .125 to 3D for the number 
of occurrence in samples of NV from the binomial with » = 4 for NW = 3(1)100. 


A-12a. A chart showing the power curves for one- and two-sided tests for 
means and differences of means from normal with ¢ known for the significance 
levels .01 and .05 for two-sided tests and .005 and .025 for one-sided tests and 
for values of the test criterion in the range 0 to 4.5. 


A-12b. For the same sampling situation as in A-12a, but with o known or 
unknown, values of the test criterion to 2D or more required to make the 
power of the test .1(.1).9, .95, .975, .99, .995 for degrees of freedom 4(1)10, 
12, 16, 24, 36, © and for significance levels .005, .0125, .025, .05 for one-sided 
tests and .01, .025, .05, .1 for two-sided tests. 


A-12c. For the same sampling situation as in A-12b, the sample sizes needed 
in the case of one mean and the equal sample sizes needed in the case of two 
means to attain the power .5(.1).9, .95, .99 for test criterion = .1, .2(.2)2, 3 
for the significance levels, .005, .125, .025, .05 for one-sided tests. 


A-19. Distribution of the signed rank statistic 7. For samples of N = 1(1)20, 
the 100, and 100 (1 — a) percentile values of T for a’s in the range .005—.125. 
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A-20. Distribution of the rank 7’. For all pairs of sample sizes from (1, 1) to 
(10, 10), the cumulative probability distribution of JT’ is given to 3D. 


A-28. Values of p(1 — p), Vp(1 — »), 1 — p?, 1 — (1 — p)?, 2 arcsin Vp and 
2 arcsin ¥1 — » to 4 and 5D for p = .01(.01)1. 


A-29a. Values of (*) a — p)*-* to 4D for N = 2(1)10 and p = .01, 
.05(.05).5 and 4. 


A-29b. All values of co for N = 1(1)20. 


A-29c. Values of N! to 5S for N = 1(1)100. 


A-30a. 95, 97.5, 99, 99.5, and 99.95 percentiles of the distribution of the coeffi- 
cient of correlation to 3D in samples of N = 5(1)20(2)30(10)60, 80, 100, 250, 
500, 1000, © from an uncorrelated normal bivariate universe. 

1+r 


i-r 





A-30b. Values of $ In to 5D for r = —.99(.01).99. 


A-30c. Selected values of =d7, used in the rank correlation coefficient, that 
will be equalled or exceeded with probabilities < .175 in samples of N = 4(1)10. 


A-30d. The probabilities to 3D that the quadrant sum S, for the corner test 
for association, will be equalled or exceeded in samples of 2N, for N = 3(1)5, 
7, © and S = 8(1)23. 


A-31. Reciprocals of N to 5D for N = 1(.01)1.1(.1)9.9. 
A-33. 4D values of VN for N = 1(.01)(10)(.1)99.9. _ 
Tables in the new edition which have been significantly extended beyond their 


counterparts in the first edition are: 


A-5. Percentiles of the ¢ distribution. 60, 70, 80 and 90 percentile values have 
been added and 99.75 percentile values deleted. 


A-6b. Percentiles of the x*/df distributions. .05, .1, 20(10)80, 99.9 and 99.95 
percentile values have been added and the table now gives values for 
1(1)20(2)30(5)60(10)100(20)200(50)500, 750, 1000, 5000 and ~ df. 


A-7c. This is a much extended table of percentiles of the F distribution. For 
the degrees of freedom, »; = 1(1)12, 15, 20, 24, 30(10)60, 100, 120, 200, 500, 
and ve = 1(1)12, 15, 20, 24, 30, 40, 60, 120, ©, the values of F to 3S or 3 or 
more D are given for which the cumulative probabilities are .0005, .001, .005, 
.01, .025, .05, .1, .25, .5, .75, .9, .95, .975, .99, .995, .999, and .9995. 


A-8b(1). This now gives for samples of 2(1)20 from normal, unbiased esti- 
mates of ¢ as a multiple of the sample range to 3D, the variance of the esti- 
mate to 3S and its efficiency to 3D. 


A-8b(3). Modified linear estimates of ¢ in samples of NV from normal are now 
given for N = 2(1)20 with the coefficients to 4D together with their variances 
to 3S and the efficiencies to 3D. 
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to A-8c(1). Upper (and lower) 95, 97.5, 99, 99.5, 99.9 and 99.95 percentiles of 
the ratio of the deviation of the sample mean from the population to the range 
d in samples of N from normal for N = 2(1)20 are given to 3D for the first four 


percentiles and to 2D for the last two. 


A-8c(2). In the new tables the same set of percentiles is given as in A-8c(1) 
a and the same D for the difference of two sample means divided by the average 
of the ranges for sample of sizes N; = Nz = 2(1)20. 


A-8d. For the ratio of the ranges in samples of N; and N; from normal 
(substitute F ratio), .5, 1, 2.5, 5, 95, 97.5, 99 and 99.5 percentiles are now given 
to 2 or 3S again for N;, Nz = 2(1)10. 


A-8e. Criteria for testing for extreme mean. Table 8e of the first edition is 
i- extended to give critical values at the .5%, 2%, 10%, 20% and 30% levels as 
0, well as the 1% and 5% levels. 


A-11. This differs in form from Table 11 of the first edition. It now gives 
the complete distribution to 3D of the total number of runs in samples of 
(Ni, N2) for Ni, Ne = 2(1)10. For N; = Nz = 11(1)20(5)100, the table gives 


t the number of runs that occur with probabilities £.005, .01, .025, .05, .95, 

).. .975, .99, .995, the mean number of runs to the nearest integer, and the vari- 

‘ ance and the standard deviation to 2D. 

; 

5, A-13. Table 13 of the first edition is replaced by charts for the power of the 
analysis of variance test. If for K samples of m each from N(;, o*); i = 1, 
--+, KX; we designate by ¢g* the average squared deviation of the u,’s from 
their mean divided by o?/n, the charts enable one to read the power of the 
test with significance levels .01 and .05 for the appropriate range of values of ¢ 

c for the degrees of freedom, m = 1(1)8 and mz = 6(1)10, 12, 15, 20, 30, 60, @. 
A-14. The values of In (1 — 8)/a to 3D are now given for a, 8 = .001, .005, 
01, (.01).05(.05).25. 

e 
Table 15. The values of the cumulative Poisson distribution to 3D are now 

5 given for Np = .05(.05)1(.1)2(.2)4. 


r A-16. Tolerance factors for normal distributions. This is the same as table 16 
of the first edition with the exception that the tabulation is for sample sizes, 
N = 2(1)27, 30(5)100(10)200, 250, 30C(100)1000, instead of for N = 2(1)- 
102(2)180(5)300(10)400(25)750(50)1000, 2. 


r A-18. Percentiles of the distribution of the ratio of the range to the standard 
, deviation (as defined in this text) in samples of K from normal. Upper 95 and 
99 percentiles are given for y = 1(1)9 degrees of freedom and upper and lower 
values for vy = 10(1)20, 24, 30, 40, 60, 120, ©, all to 3S for K = 2(1)20. 


- Table 19. Determination of the second sample size, of the first edition is 
omitted from the second edition. 


v This extensive and useful set of tables in itself very considerably enhances 
s the value of the second edition as compared to the first. 


Co <a. 
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98[K].—B. L. vAN DER WAERDEN, Mathematische Statistik Springer-Verlag, 
Berlin, Germany, 1957, ix + 360 p., 23 cm. Price DM 49.60. 


This is an example of a book written in the careful style of van der Waerden. 
It includes at the beginning a chart of the type used in his famous Algebra indicate 
earlier chapters which are essential to an understanding of later ones. At the end 
it includes an index of examples classified according to the field in which they 
were generated, and there is a short English-German glossary. 

The purpose of this review is to announce tables in the bodk, and not par- 
ticularly to discuss material covered. However, it might be noted that the author 
describes the central question of statistics as that of how far quantities estimated 
from random samples may vary from the corresponding ideal values (translation 
approximate only). To attack this question he studies the axiomatic development 
of probability and its applications to statistics, including studies of confidence 
limits and confidence bands. He then turns to the theory of estimation, and then 
to the testing of hypotheses. 

Tables include the normal distribution (area) @(#), ¢ = —3.9(.1)—3(.01)- 
3(.1)3.9, 4D, the inverse function for the normal distribution V(x), x = 0(.001).99, 
2D; g and g?, the (one sided) .05, .1, .5, 1, 2.5 and 5% points for the normal 
distribution and the x? distribution with one degree of freedom, 2D; exact and 
asymptotic 1 and 5% points for Smirnoff’s test, m = 5, 8, 10(10)50, 4D and the 
ratio of the values 3D; exact and asymptotic 1 and 5% points for Kolmogorov’s 
test for m = 5(5)30(10)80 and the 5% values for m = 90, 100, 4D with ratios 3D; 
.1, 1, and 5% points of x? with f degrees of freedom, f = 1(1)100, 3S with suitable 
formula for larger f; (one sided) .05, .5, 1, and 2.5% points for Student’s test 
with f degrees of freedom, f = 1(1)30(10)60(20)100, 200, 500, », 4S; 5% points 
for F-test with f; and f. degrees of freedom in numerator and denominator re- 
spectively, f; = 1(1)20(2)30(10)60(20)100, fz = 1(1)30(2)50(10)100(25)150, 200, 
300, 500, 1000, 3D; the 1% points with f, as above and fz = 2(1)30(2)50(5)60- 
(10)100(25)150, 200, 300, 500, 1000, and the .1% points for f; = 1(1)10(5)20, 30, 
50, 100 and f2 = 2(1)20(2)30(10)60(20)100, 200, 500, «, 3S; (one sided) .5, 1 
and 2.5% points for the sign test, m = 5(1)100, integral; tables for Wilcoxon's 
test of the number of inversions in a series with g marks x and h marks y with 
g = 2and hk = 5(1)10, and g = 3(1)10, & = g(1)10, the probabilities not exceed- 
ing 5% are listed 3S; tables for the X-test, (one sided) .5, 1, and 2.5% points 

r 
for X = x (- +1 
g — h = 0(1)5, 3S, where g and h are as in Wilcoxon’s test above, n = g + h, 
r is the set of rank numbers of x in the sequence of x and y, and y is the inverse 





). and corresponding values of X for m = 6(1)50 and 





1 n 
normal distribution function; also an auxiliary table of Q = rv ( . x -) , 
r=] 


n = 1(1)150, 3D; (one sided) .1, 1, 2, and 5% points for the correlation coefficient, 
f =m — 2 = 1(1)20(5)50(10)100, 3D. Many of these are reprinted from other 
sources, which are acknowledged. 


.. a &. 
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99[K ].—H. A. Davin, “The ranking of variances in normal populations,”’ Amer. 
Stat. Assn., Jn., v. 51, 1956, p. 621-626. 


This paper considers two procedures for ranking, or grouping, the variances 
of? (t = 1, ---,%) of k normal populations by use of mean square estimates s7 
each based on » degrees of freedom. 

Let X; < X2 <---< X; be an ordered sample of k independent observations 
from a unit normal distribution. Table I gives, to 3S and for k = 2(1)12, values G 
such that the expected number of differences, X14; — X:, which exceed G is 
a = .05 and a = .01. The title of Table I refers to these as “Critical 100a% 
levels . . .’’ and the text material points out that they are not percentage points. 
If X, = log s? Table I can be used to rank the variances o/. 

Table II gives, to 3S, for k = 2(1)12, for » = 2, 4(1)10, 12, 15, 20, 30, 60, @, 
and for a = .05, a = .01, values of Ra(v) such that the expected number of 
ratios s*,,;/s? (when the s? have been ordered by size) which exceed Ra(v) is a. 
Again the title refers to ‘Critical 100a% levels. . . .’’ The title has a misprint, 
the division bar being printed as m. 

Table III (titled “Significant Maximum F-Ratios in a 1% Level Multiple 
Fmax-Ratio Test. (Normal Variation Assumed)’’) gives, to 3S, for a = .01, for 
k = 2(1)12, and for vy = 2(1)10, 12, 15, 20, 30, 60, ©, values Ca(v,k) which 
are the upper 100{1 — (1 — a)*'}% points of the maximum ratio of the s/?. 


F. J. MAssEy, Jr. 
University of California 
Los Angeles 24, California 


100[K ].—D. G. CuarMan, “Estimating the parameters of a truncated gamma 
distribution,”’ Ann. Math. Stat., v. 27, 1956, p. 498-506. 


To facilitate computation of maximum likelihood estimates of parameters of 
the gamma or Type III distribution with density function 


S (x) = i, e7etz—e) (x on c)*; x>c 
r'(d) rer et 


= 0; a4" 
in the case c = 0, a table of the function 
b = y"[In  — Inx], 


where £ designates the sample mean and In x designates the mean of natural 
logarithms of the sample observations, is given for [In # — In x] = .01(.001).99, 
.1(.005).495, .50(.01)1. For values of the argument from .01 —.049, b is given to 4D. 
Elsewhere it is given to 3S. The author states that a complete tabulation of (5) 
with its first and second differences, correct to one figure in the fifth decimal, 
for 6 = .01(.01)2, 2(.02)5, 5(.1)20, 20(1)100, is available in mimeographed form 
from the Laboratory of Statistical Research, University of Washington. 

A new procedure bearing certain resemblances to a minimum chi square 
method is suggested for estimating parameters of a truncated gamma distribution. 
The asymptotic variance-covariance matrix of these estimates is determined. 
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The author points out that his method is applicable to a number of other truncated 
distributions whether truncation is in the tails or the center of the distribution. 


A. C. COHEN, JR. 
The University of Georgia 
Athens, Georgia 


101[K].—R. J. BuEHLER, ‘“‘Confidence intervals for the product of two binomial 
parameters,’ Amer. Stat. Assn., Jz., v. 52, 1957, p. 482-493. 


The author considers the problem in which two binomial populations have 
been sampled and that &; failures are observed in a sample of m; from the first 
population and k; failures are observed in a sample of m2 from the second and one 
is interested in finding the number C,,n,(%:, 2; a) having the property that 


Prob {0 < PiP2 < Cain, (Ri, 23 a)} >a forall O< Pi, P2 <1 


where P; and P; are the probabilities associated with the respective statistically 
independent populations. That is, one is interested in obtaining an upper con- 
fidence interval in the sense of Neyman for the product of two probabilities of 
statistically independent binomial populations. 

The author proposes that the above problem be solved by taking 


Can, (hi, Re; a) = Ca, (her, Ver) -Cx, (hee, Vex) 
where Cy,(ki, Va) and C,,(k2, Va) are such that 
Prob {0 < Pi < Ca,(Ri; Va)} > Va for O< P; < 1,4 = 1,2. 


Approximate solution to the above inequalities are obtained under the assumption 
that m is large and C is small by using Poisson terms for the binomial terms. 
These approximations are tabulated. 

Table 1 gives values of 2C,,(k; a) to 3S for k = 0(1)6, 8, 10, 20, 30 anda = .90, 
.95 and .99. 

To obtain the confidence limit C,,,n,(21, 2; a) which are tabulated in Table 3, 
the author resolves the lack of uniqueness problem associated with the solution 
of the inequality by imposing a generalized condition of (1) Regularity: If k; < ke 
then C,(R1;a) < Cn(R2,a) and (2) Shortness: C,(k;a@) should be as small as 
possible. The 3S values for 2%2Cn,n,(R1, R2; a) are tabulated for k; = 0(1)5 and 
ke = 0(1)29 for a = .90 and .95. The symmetry allows the interchange of k; 
and ke. 

The final Table 5 gives unsymmetrical values for the confidence limits 
NyN2Cn,n, (Ri, 22; 0.9) to 1D for ki, ke = 0(1)4 which are shortened by use of an 
auxiliary random element in a fashion similar to that suggested by Stevens [1] 
and Eudey [2]. 


C. F. Kossack 
Purdue University 
Lafayette, Indiana 


1. W. L. StEvEns, ‘‘Fiducial limits of the parameter of a discontinuous distribution,” Bio- 
metrika, v. 37, 1950, p. 117-129. 

2. M. W. Eupey, “On the treatment of discontinuous random variables,”’ Stat. Lab. Technical 
Report No. 13, Univ. of Calif., Berkeley, 1949. 
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102[K ].—Josrern Berkson, ‘Tables for the maximum likelihood estimation of 
the logistic function,’ Biometrics, v. 13, 1957, p. 28-34. 


Table 2 in this paper gives to 5D the function p = 1/[1 + e~*'] and (directly 
beneath the value of ») the function w = p(1 — p), for argument / = 0(.01)4.99. 
The text states that the tables are correct to +1 in the fifth place and that linear 
interpolation in the table will also give values correct to +1 in the fifth place. 

Suppose P(x) = 1/[1 + e~“+*»)]. Then P(x) is said to be a logistic function 
of x. If, for each of a set of values x1, ---, xz, there are performed n; independent 
binomial trials with probability of success P; = P(x,), the results may be used 
to estimate the parameters a and @ in a number of ways. The table of values 
of the functions p and w facilitates the calculation of the Maximum Likelihood 
estimates of a and 8. The text explains the method and gives a numerical example 
(table 1). 

PauL MEIER 
University of Chicago 
Chicago, Illinois 


103[K ].—A. J. Duncan, “Charts of the 10 per cent and 50 per cent points of 
the operating characteristic curves for fixed effects analysis of variance F 
tests, a = 0.01 and 0.05,”” Amer. Stat. Assn., Jn., v. 52, 1957, p. 345-349. 


The symbols and expressions a, 8 (ordinate of OC curve), f1, f2, OC (1-power) 
have meanings usual in testing hypotheses and (fixed effect) variance analysis. 
Duncan gives four (very small) charts for ¢ (defined below) 


(1) a@ = .05,6 = .10, fi = 1(1)8, fe = 6(1)10, 12, 15, 20, 30, 40, 60, @ 


(2) a@= .05, 8 = .50, fi as above, fz as above 
(3) a= .01, 8 = .10, f: as above, fs as above 
(4) a= .01, 8 = .50, fi as above, f2 as above 


where ® is defined for (a) one-way analysis, (b) two-way analysis, (c) Latin 
Squares. For example, for a Latin Square, 


# = L 6,?/0? 
k 
for testing treatment effects, the basic linear model being 


Xiyg = wt ti + 0; + Oe + espe. 


These charts (described by the author as “essentially a special condensation of 
the Pearson and Hartley charts’’) should be examined along with similar (power) 
tables and charts by Fox, [1], Lehmer [2], Pearson and Hartley [3], and the 
pathbreaking paper by Tang [4]. 

HAROLD FREEMAN 


Massachusetts Institute of Technology 
Cambridge, Massachusetts 








238 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


Martin Fox, “Charts of the power of the F-test,” Ann. Math. Stat., v. 27, 1956, p. 484-497. 
CMPAC. v. 11, Rev. 28. 7] 

2. EMMA LEHMER, “Inverse tables of probabilities of errors of the second kind,” Ann. Math. 
Stat., v. 15, 1944, p. 388-398. 

3. E. S. Pearson & H. O. HARTLEY, “Charts of the power function for analysis of variance 
tests derived from the noncentral F-distribution,” Biometrika, v. 38, 1951, p. 112-130. [MTAC, 
v. 8, 1954, RMT 1183, p. 84.] 

4. P.C. TANG, “The power function of the analysis of variance test with tables and illustrations 
of their use,”’ Stat. Res. Memoirs, v. 2, 1938, p. 126-149 + tables. 


104[K ].—Bruno DE FinettT1, ‘Una legge riguardante l’estinzione nei processi 
di eliminazione,” Giornale dell’ Instituto Italiano degli Attuari, v. 15, 1953, 
p. 94-99, 


A population, initially of m elements, is subjected to successive trials with 
elimination probability g (for each individual, each trial) until complete extinc- 
tion. The author gives a recurrence formula and an exact series formula for 
P,,” = the probability that the number of survivors eliminated at the last step 
(before extinction) is m, m = 1,2, ---,m. An approximate formula is derived, 
Pa= g"/m|log (1 — g)|, which, [1], is not the limit of P,," as » becomes infinite. 
As is shown in [1], P,," approaches a periodic function of log m for whose values 
P,, is an average for large n. 

A table of P,,” is given for the cases through m = 10 for g = 4. Exact expres- 
sions are indicated and values (quotients) are given up to maximum 6S. There 
is an error of +1 in the least significant digit for P.®. Round-off is uniform over n 
for fixed m but (starting at » = 3) works to the disadvantage of >, P.” = 1. It 
is of interest to note that for the case of g = 4: P} = P3; P5° < Py’ < P}® < Pit. 


W. H. Martow 
The George Washington University 
Logistics Research Project - 
Washington, D. C. 


1. Grusepre Orttaviant, “A proposito della legge di estinzione nei processi di eliminazione,”’ 
Giornale dell’ Instituto Italiano degli Attuari, v. 15, 1953, p. 100-114. 


105[K, W, Z].—A.tson Dore, ‘A bibliography on the theory of queues,” Bio- 
metrika, v. 44, 1957, p. 490-514. 


An extensive and valuable bibliography on the theory of queues including 
related topics and papers devoted to computational aspects of the subject. 


oe 


106[L].—K. A. Karpov & S. N. Razumovsk1, Tabliisy integral’nogo logarifma 
(Tables of the logarithmic integral). Akad. Nauk SSSR, Vychislitel’ny Tsentr, 
Matematicheskie Tabliisy (Acad. Sci. USSR, Computational Center, Mathe- 
matical Tables). Izdatel’stro Akad. Nauk SSSR (Press of the Academy of 
Sciences of the USSR). Moscow, 1956. 320 p., 26.7 cm. Price 33.30 rubles. 


The main table (pages 11-311) lists values of the logarithmic integral 


liz = [Panos (x > 0), 
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understood as a principal value when x > 1, to 7S for 
x = 0(.0001)2.5(.001)20(.01)200(.1)500(1)10,000(10)24,990. 


There are no differences except in the first five pages (x < .5), where first and 
second differences are provided. 

Pages 312-315 contain an auxiliary table of lix — In|1 — x| to 6D without 
differences for x = .95(.0001)1.0499. Page 316 tabulates 4#(1 — 7) to 4D without 
differences for ¢ = 0(.001)1. 

The main table was computed on the high-speed electronic calculating machine 
(BESM) of the Academy of Sciences. It was found expedient to use, chiefly, 
numerical quadrature by Simpson’s rule, reserving the formula lix = Ei (In x) 
for the computation of pivotal and checking values. The error is stated not to 
exceed 0.6 final units. A correction slip pasted in gives the following correction 
on p. 186: 


li 25.82, for 0.1175643-10? read 0.1176643-10*. 


Although the range 0 < x < 24,990 does not include the high values (up to 
at least x = 10*) of interest in number theory, and although extra figures may be 
found elsewhere for a few arguments, this is nevertheless by far the largest table 
of the logarithmic integral available. It will be noticed that throughout the 
interval .0100 < x < 24,990 there is always at least a three-figure argument, and 
that the main table runs to 300 pages. None of its predecessors comprises more 
than a page or two. In these circumstances one cannot seriously dissent from the 
claim that this is the first table of the logarithmic integral which is anything like 
complete. It can hardly fail to be found of use. 

On pages 318-319 are given corrections to four other tables published in the 
same series in 1954 and 1955. 

A. F. 


107[L].—E. N. Dexanosinze, Tabliisy isilindricheskikh funkisit ot doukh pere- 
mennykh (Tables of cylinder functions of two variables). Akad. Nauk SSSR, 
Vychislitel’ny Tsentr, Matematicheskie Tabliisy (Acad. Sci. USSR, Computa- 
tional Center, Mathematical Tables). Izdatel’stvo Akad. Nauk SSSR (Press of 
the Academy of Sciences of the USSR), Moscow, 1956. 494 p. + 2 leaves of 
diagrams, 26.7 cms. Price 50 rubles. 


These tables were computed at the Institute of Precise Mechanics and Com- 
putational Technique on the initiative of Professor P. I. Kuznetsov. The functions 
tabulated are Lommel functions of two variables, U,(w, z), V,(w, 2), which are 
related to Bessel functions and were first introduced in connection with the study 
of diffraction. Apart from the use of » and w instead of m and y, the functions are 
as described in MTAC, v. 1, 1944, p. 258 and FMR Index, p. 301. 

The tables (pages 11-492) are for y = 1 and 2, the four functions U;, U2, 
Vi, V2 being listed in parallel columns. All values are to 6D and there are no 
differences. There is a separate table for each of the following values of w: 


w = .5(.02)1.2(.05)4(.1)10. 
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For w < 6.2, the second variable takes the values of z = w(.01)W, where W is 
always fairly close to 4¥w; for w > 6.3, on the other hand, z = w(.01)10. The 
resulting region of tabulation in the (w, z) plane is illustrated in Figure 3. 

Various formulas are given for extending the scope of the tables. For example, 
computation for suffixes vy which are integers greater than 2 is possible if the 
present tables are supplemented by tables of Bessel functions J of integral order; 
negative values of w and z are covered by symmetries; while a further valuable 
resource is that, for integral », the U and V functions are interchanged if w is 
replaced by 2*/w. Asymptotic expansions for large w and z are given. Some 
numerical examples are worked out. 

The computation of the tables is briefly sketched in mathematical terms, but 
nothing is said about the machine or machines used. The errors are stated not to 
exceed unity in the sixth decinal place. The mere four pages of the Introduction 
contain so many formulas that they will be largely intelligible to mathematicians 
without a knowledge of Russian. 


A. F. 


108[L, S].—Kurt ALpER & AaGE WINTHER, Tables of the Classical Orbital 
Integral in Coulomb Excitation, Mat. Fys. Medd. Dan. Vid. Selsk., Bind 31, 
No. 1, 1956, 74 p., 9 cm. Price 3.00 kr. 


This paper contains a tabulation of certain functions which arise in the 
classical description of Coulomb excitation. The functions tabulated here are the 
orbital integrals J,,(v, ¢), defined by 


[cosh w + ¢ + iVe — 1 sinh w}*dw 
[e cosh w + 1 }+# : 





IK(v, ) = f 6 { eile sinh w+) } 


the differential cross-section functions 





dfen(v, $) dfun(v,f) _. 
a” and q Biven by 
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-|Denu(v, £) |?-cot? »/2 -sin-* »/2; 


and the total cross-section functions fm(v, ¢) and fin(», ¢) given by 


fm(§) = 24 [eee sin ydy 


fur(t) = 2e f PD sin ‘ile. 


The tables of the orbital integrals were obtained by numerical integration using 
the high-speed electronic computer BESK at Stockholm. 
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The orbital integrals were calculated for the following values of the parameters: 


vy = 10°(10°)180° 
and 
¢ = 0.0(0.1)1.0(0.2)2.0,4.0 for = 1 and 2 
¢ = 0.2(0.2)1.0(0.5)2.0,4.0 for A = 3 and 4, 
w= —A, —A +2, ---X. 


The tables are said to be accurate to the number of decimals given, which is six, 
except for small »(10°). 

The differential and total cross section functions are tabulated for the following 
parameter values: 


ae g) and fm(f) for \ = ® ; 3, and 4 
“ae 8) g) and fin(f) for \ = 1 and 2. 


In the introduction to the tables properties of the orbital integrals are de- 
scribed. A review of the Coulomb excitation process wherein these integrals are 
discussed is given in the following publication: K. Alder, A. Bohr, T. Huus, B. 
Mottelson and A. Winther, Rev. of Mod. Phys. v. 28, 1956, p. 432. 

Copies of these tables may be obtained from 


Secretary 

Det Kongelige Danske Videnskabernes Selskab 
Dantes plads 5 

Copenhagen V, Denmark, 


at a cost of 10.00 kr. 


Lioyp D. Fospick 
University of Illinois 
Urbana, Illinois 


109[L, V].—Jack N. Nre.sen, Tables of Characteristic Functions for Solving 
Boundary- Values Problems of the Wave Equation with Application to Supersonic 
Interference, Report NACA TN 3873, 1957, 245 p., diagrams and Tables, 
26 cm. 


The present Technical Note deals with the problem of calculation of the 
pressure distribution on quasi-cylindrical bodies of revolution including the 
different possible cases of interference. The principles of the referring theory had 
been outlined in detail in several previous reports of Jack N. Nielsen and other 
authors (see References). This theory contains the so-called W,,-functions, which 
had been tabulated for different values of the quantities x (streamwise distance) 
and 1 (radial distance of unity corresponding to the body surface), Wo occurring 
in the formula for pressure coefficient at the body surface, the higher order 
W-function being needed in connection with a calculation of interference pres- 
sure fields. 
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Nielsen stated some disadvantages of these W,,(r, 1)-functions, if the external 
pressure field (and not only the pressure at the body surface) was to be calculated. 
In this case singularities occurred and the results obtained as differences of two 
large numbers became inaccurate. Therefore a new method was developed based 
on a more general set of functions W,,(x, r), r being the radial distance from body 
center line. Nielsen found that with this new function it is as easy to compute 
the pressures away from the body as those on the body. 

The new W,,(x, r) functions are defined as the inverse Laplace transformation 
of a function containing modified Bessel functions of the second kind: 


= J-1 rl K.,(sr) =| 
Wa (x,r) = L [-« ete ; 





They are the solution of the partial differential equation 
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with the boundary condition 


_ (1/8) — (m*/r) 4 (3/8) + (m?/r) 


r r3 


Wn 





n=x-r+1=0 


OWn 
or 





=0 at r=1. 


The author discusses the different methods of evaluating this formulae par- 
ticularly adapted to automatic computation. He found that the following 3 
methods if used in a combined manner are to be considered as the most suit- 
able one: 


1. Power series in x with coefficients depending on m and r 
2. Integral-equation methods 
3. Characteristics using finite difference equations 


Thus the calculated tables of W,,(x, r) cover a range of the order m from 1 
to 10. The range of x is less for the higher order functions than the lower order 
functions because they become smaller faster. The actual ranges of the parameter 
x and r are 


r (all values for each value 





m x m and r) 

0 1 0 to 10 L, Bhs. dee be 
: a 0 to2 2.0, 3.0, 4.0 

4 5678910 0to5 6.0, 8.0, 10.0 


On the whole 69,000 values are contained in the prescribed tables. 
The author has intensively investigated the question of the accuracy and 
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came to the following results; 
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As to the physical significance of the W,,(x, r) function, the author explains that 
these functions are to be considered as cylindrical pressure waves associated with 
a step in the radius of a streamwise quasi-cylindrical body. 

The author finally interprets the application of the W,,(x, 7) functions to some 
aerodynamic examples which partly have hitherto not appeared in the literature, as 


1) Calculation of the pressure field due to an indented body 

2) Flow field of an infinite cylinder impulsively moved in axially direction in 
a viscous incompressible field 

3) Wave drag of corrugated bodies 

4) Arbitrary quasi-cylindrical body alone 

5) Body-body or shock-body interference 

6) Wing-body interference 

7) Wortex-panel interference 

8) Non-aerodynamic applications 


The reviewed Technical Note represents without doubt a very remarkable 
contribution for a general solution of boundary value problems of the second kind 
involving the wave equation in three dimensions with approximately circular 
cylindrical boundaries or involving the unsteady heat-conduction equation in two 
space dimensions with nearly circular boundaries. 

Reference is made on some publications which deal with the same or with 
similar problems; listed below are some of the most important. 

ALFRED WALz 
National Bureau of Standards 
Washington, D. C. 


1. Jack N. NIELSEN, “Quasi-cylindrical theory of wing-body interference at supersonic speeds 
and comparison with experiment,”” NACA Report 1252, 1955. : 

2. Jack N. Nretsen, “General theory of wave drag reduction for combinations ame 
oe cm bodies with an application to swept-wing and body combinations,” NACA T 
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3. M. J. LicuTumtt, “Supersonic flow past bodies of revolution,” R. & M., No. 2003, British 
A. R. C., January 3, 1945. 

4. D. G. RANDALL, “Supersonic flow past | bodies of almost circular cross- 
British R. A. E. TN Aero 2404, November, 1 
_5. G. K. BATCHELOR, ‘The skin friction on infinite ~ moving parallel to their length,” 
. R. C. 15,105, Fluid Motion Subcommittee F. M. 1772, Aug. 6, 1952. 
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110[P, Z].—W. W. Dosrowotsk1, Theorie der Mechanismen zur Konstruktion 
ebener Kurven, Akademie-Verlag, Berlin, 1957, viii + 134 p., 24 cm. Price 
DM 19, bound. 


This work was part of an extensive symposium on mechanisms and machines 
and appeared in Trudy Seminara Teorii po Mashin i Mekhanizmov 9, no. 36 
(1950), Akademiya Nauk SSSR. 

The constructions of projective geometry are implemented in mechanical 
devices to draw curves. A point constrained to a line is replaced by a pin con- 
strained in a straight slot. A line passing through a fixed point is replaced by a 
slotted bar pivoted about a fixed pin. The transfer of angular rotation from one 
line to another is accomplished by locking the opposite sides of a parallelogram 
linkage to the rotating lines. The demonstrations are quite detailed and elemen- 
tary with very little use of algebra. Some of the mechanisms have more theoretical 
than practical value. 

A conic section can be drawn by a mechanism which simulates Pascal’s 
theorem. If five points A, B, C, D, E of the conic are given, the variable sixth 
point M traces the conic. Fixed slots through BE and EC are used to guide pins 
which move, also, in three rotating slots which are hinged at H (the intersection 
of AB and CD) and points A and D. If the pin L lies on the triplet of lines (BE, 
HL, DL) and the pin K on (EC, HL, AK), then the pin M will be the intersection of 
AK and DL. Variations of the mechanism are shown corresponding to special cases 
of the theorem, namely, when any five elements are given (points and tangents). 

Third and fourth degree curves are generated by compounding the simpler 
quadratic linkages, or by the use of four-bar and crosshead mechanisms. Mecha- 
nisms for quadratic transformations and Cremona transformations are described. 
Some higher degree curves, like the Lissajous figures, can be made by a combina- 
tion of two functions of the same parameter. An inversor mechanism can convert 
a curve into another of higher degree. 


MICHAEL GOLDBERG 
Bureau of Ordnance 


Navy Department 
Washington, D. C. 


111[S, X, Z].—S. FLucce, Editor, Handbuch der Physik (Encyclopedia of Physics), 
Band 1, Mathematische Methoden I (Mathematical Methods 1), Springer- 
Verlag, Berlin, 1956, vii + 364 p., 24.5 cm. Price DM 72. 


[S, X, Z].—S. FLucce, Editor, Handbuch der Physik (Encyclopedia of Physics), 
Band 2, Mathematische Methoden II, Springer-Verlag, Berlin, 1955, vii + 520 
p., 24.5 cm. Price DM 88. 


The two volumes contain as complete a compendium of useful mathematical 
methods and theorems as one is ever likely to compile at any one time as part of 
an encyclopedia of physics. It would be impracticable to give more than a brief 
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indication of the contents and depth of the separate sections prepared by various 
authors. 

The first three sections of Volume I were written by Josef Lense and cover (a) 
calculus, Lebesgue integration, complex variable theory, ordinary differential 
equations—89 pp. ; (b) partial differential equations culminating in the classifica- 
tion of the semi-linear second order equation in two independent variables as 
elliptic, parabolic or hyperbolic and the Riemann method for solving the linear 
hyperbolic equation—30 pp. ; and (c) elliptic functions and integrals—27 pp. 

The next section on the special functions of mathematical physics was written 
by J. Meixner and gives a condensed but thorough development of the properties 
of all of the commonly used functions and relates them to the study of partial 
differential equations by the method of separation of variables—71 pp. 

The final section of Volume I, by F. Schlégl, gives a detailed account of the 
techniques for studying the boundary value and eigenvalue problems for second 
order partial differential equations and includes the related aspects of the theory 
of integral equations and the calculus of variations (use is made of the Dirac delta 
function in studying the Green’s function). Most of the material relates to elliptic 
equations and less space is devoted to the parabolic and hyperbolic equations— 
135 pp. 

Volume II begins with a thorough treatment of selected topics in modern 
algebra by G. Falk which contains a discussion of linear algebra, the theory of 
group representations and the representation theory of algebras, with some brief 
indication of the application to mechanics—116 pp. 

The second section is by H. Tietz and covers vectors, analytic geometry, 
projective geometry, differential geometry of curves and surfaces, Riemann 
geometry, tensor analysis, theory of spinors, contact and canonical transforma- 
tions for a Hamiltonian system of equations—81 pp. 

The next section by I. N. Sneddon is the only one written in English and is 
devoted to a detailed study of many integral transforms, their properties and 
their application to typical physical problems. It is supplemented by a discussion 
of Lebesgue integration, a development of the theory of linear operators and 
functionals in Hilbert and Banach spaces, and a short description of the theory 
of distributions (of L. Schwartz)—S51 pp. 

The section on numerical and graphical methods is written by L. Collatz. He 
discusses the operation of the slide rule, nomography, smoothing of data, approxi- 
mations of functions, location of roots of polynomials and more general functions, 
solution of linear and non-linear systems of algebraic equations by direct and 
iterative methods, determination of eigenvalues of matrices; finite differences, 
interpolation and numerical integration ; numerical methods for ordinary differ- 
ential equations, initial value problems and mixed initial-boundary value prob- 
lems for partial differentia] equations, methods of Ritz and Trefftz; integral and 
functional equations. The descriptions of methods are accompanied by adequately 
worked out examples—122 pp. 

The last section of 28 pages is by H. Biickner and describes the operation and 
logical structure of analogue and digital computers and the coding of problems. 

The above studies are of high quality and should serve as a valuable guide to 
the branches of mathematics that are most usable today in classical and modern 
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Physics. The treatments are generally self contained, but references to the 
mathematical literature are given to make it possible to keep the size of these 
volumes manageable. Mention is made of the fact that the subject of Probability 
Theory is not contained in these books. Otherwise the two books are a well- 
rounded and remarkably well coordinated set of chapters in what might be called 
Applied Mathematics in spite of its high purity. 

E. I. 


112[V].—J. J. Stokes, Water Waves, Interscience Publishers, New York, 1957, 
xxviii + 567 p., 23 cm. Price $12.00. 


This book gives a lucid, well written, and connected account of the theory 
of wave motion in incompressible fluids with a free surface and subjected to 
gravitational and other forces. The author pays special attention to some of the 
more recent work in the field to which he has made significant contributions during 
and since World War II. He has taken considerable pains to supply background 
material in hydrodynamics, some of the mathematical toois needed and an 
extensive bibliography. 

There are four main parts of the book: Part I, comprising Chapters 1 and 2, 
presents the basic hydrodynamic theory for non-viscous incompressible fluids. 
It also contains the derivation of the equations of the two principle approximate 
theories which form the basis for the discussion in the remainder of the book. 
The first of these is a consequence of the assumption that the wave amplitudes 
are small. The second comes from the assumption that the depth of the liquid is 
small. In both cases the approximate theories are derived from formal expansions 
with respect to a small dimensionless parameter and methods are given for de- 
termining the higher order terms. The derivation of the shallow water theory 
follows that first given by K. O. Friedrichs. 

Part II, which consists of Chapters 3 to 9 inclusive, describes the linear theory 
obtained when it is assumed that the wave amplitudes are small. The main mathe- 
matical tool used is potential theory. This part is divided into three subsections: 
(A) dealing with wave motions that are simple harmonic oscillations in time; 
(B) dealing with unsteady or transient motions that arise from initial disturbances 
starting from rest; and (C) dealing with waves on a running stream. 

At the end of subdivision A, the author outlines other methods for dealing 
with boundary problems which differ from those he discussed at some length 
earlier. In particular, the Wiener-Hopf method is briefly discussed, and the author 
“hazards the opinion that problems solvable by the Wiener-Hopf technique will, 
in general, prove to be solvable more easily by other methods. . . .” This re- 
viewer is among those who do not share this opinion. 

Part III, consisting of Chapters 10 and 11, deals with shallow water theory, 
the approximation resulting from the assumption that the depth of the liquid is 
small. This assumption leads to a system of non-linear partial differential analo- 
gous to the equations describing the motion of compressible gases. 

The first part of Chapter 10 contains a detailed discussion of the theory of 
characteristics based on the discussion given by Courant and Friedrichs in their 
book, Supersonic Flow and Shock Waves. Simple waves and propagating discon- 
tinuities are discussed. Chapter 10 contains a detailed discussion of a variety of 
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concrete problems including conditions for the occurrence of a bore and a hy- 
draulic jump, the motion resulting from the breaking of a dam, and the motion 
of frontal discontinuities in the atmosphere. It concludes with a discussion of 
some applications of the linearized version of the shallow water theory. 

Chapter 11, entitled ‘‘Mathematical Hydraulics,” deals with flows and wave 
motions in rivers and other open channels with rough sides. The fundamental 
equations derived in Chapter 10 are modified to include the force of resistance 
due to the rough sides which is taken to be proportional to the square of the 
velocity. This chapter contains an illustration of the length to which the author 
has gone to present portions of the book in such a way that they can be read, to 
a large extent, independently of the rest of the book. It contains a brief discussion 
of the method of characteristics which was given in Chapter 10; Figure 11.4.1, 
page 473, is identical with Figure 10.2.2, page 298. 

A numerical method for calculating solutions of the differential equations is 
presented briefly. The method is based on the method of characteristics, but is 
applied in the x, ¢ plane. The reader is referred to the paper of Courant, Isaacson 
and Rees for a discussion of the convergence of the method. No mention of 
stability questions is made. 

The result of calculations for the prediction of an actual flood in the Ohio and 
at its junction with the Mississippi are reported on. A comparison of the predicted 
with the observed flood is also made. The author states ‘‘there is no doubt that 
this method of dealing with flood waves in rivers is entirely feasible since it gives 
accurate results without the necessity for unduly large amounts of expensive 
computing time on a machine such as the UNIVAC.” The author also argues that 
such computations are cheaper and at least as good and as feasible as those made 
with the “analogue computers’’ previously used, namely scale models of long 
rivers or river systems. 

Part IV, consisting of Chapter 12, discusses a few analytical solutions of the 
equations describing exact non-linear theory. This chapter deals with problems 
whose solutions can be expressed in power series. The book ends with an exposition 
of Levi-Civita’s theory of progressive waves of finite amplitude in water of infinite 
depth which satisfy exactly the non-linear free surface conditions. 

This stimulating book which contains a selection of material chosen because 
it interests the author provides an excellent treatment of many topics in the 
theory of incompressible flow. It will be useful to both students and research 
workers in the field. 

A. H. T. 


113[W, Z].—Anprew D. Bootu, L. BRanpwoop, & J. P. CLEAVE, Mechanical 
Resolution of Linguistic Problems, Academic Press, Inc., New York, 1958, 
vii + 306 p., 21.5 cm. Price $9.80. 


Modern electronic computers have been operating for more than ten years. 
Despite extensive experience with them and a host of innovations, little is known 
about their limitations. 

It has been shown that some computers can play checkers, and some have 
played chess after a fashion, and some have even been forced reluctantly to learn 
a little. Each of these demonstrations throws some light on the ability of the 
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computer, and each of them up to now has resolved the question in favor of 
ability rather than of limitation. The exploration of these limits is an outstanding 
scientific problem. 

Another demonstration which in addition promises to be productive of useful 
results is that of translation from one language to another. The first of several 
demonstrations was made in 1953 by International Business Machines and Pro- 
fessor Dostert of Georgetown University. Others have followed. But it still 
remains to be shown that the machines can turn out useful translations. 

Professor Booth has been thinking of this problem since 1947, and he and his 
colleagues have put into this somewhat uneven book a resume of their work 
and ideas. 

There is a chapter on APEXC, the machine they have at the University of 
London. This is a modest machine, inadequate for the large task of translating 
languages, and the fact that they have successfully used it to demonstrate many 
of the steps of translating is a tribute to the ingenuity of the authors. 

The basic processes of making concordances and using dictionaries are de- 
scribed. Catalogues of sentence structures can be made in the same way. 

There is a chapter on the analysis of the style of Plato, a problem of interest 
to linguists, historians, and philosophers. Plato’s style changed as he grew older, 
with one striking discontinuity at about age sixty. This change in style enables 
scholars to divide the works into two groups, those earlier and those later than 
the change. Refinements in these studies of style lead to extensive burdensome 
drudgery at which a computer should excel. This is an excellently written chapter 
on a highly technical subject. 

Getting back to translating words rather than counting them, there is a dis- 
cussion of the mechanics of breaking words into stems and endings in order to 
have short dictionaries. Compound words, such as abound in technical German, 
can also be broken. If they are not broken, the size of the dictionary gets out of 
all bounds. Words with multiple translations, such as English ‘“‘but’’ into German 
“aber” or ‘‘sondern,’’ or “‘know’’ into ‘‘kennen’’ or ‘‘wissen,’’ must be handled 
with regard to the context. The two words on either side are frequently enough for 
this purpose. The comparison of the expenses of word splitting and full dic- 
tionaries is confused. 

The methods by which the computer could do the work of stem recognition, 
identification of parts of speech, and acting on the dictionary information about 
syntax are described in terms of analogous operations with punched cards, which 
are apparently thought to be more picturesque than the arcane operations of a 
computer. The most time-consuming operation will be that of consulting the 
dictionary. Considerable thought has been given to this, with the conclusion that 
the best way will be a modification of the bracketing method, in which the middle 
entry of the dictionary is consulted first, thereby establishing which half contains 
the word, then this operation is repeated on the proper half, and so on. In the 
introduction Booth says that this method was first given a practical form in 1955, 
but the reviewer is sure that versions were in use before that time. Some estimates 
of numbers of steps needed are made based on Zipf’s law (attributed to Estoup). 
Some methods of compressing words to fit in small boxes are discussed, but 
nothing about the origins of these ideas is said, nor is any evaluation of their 
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usefulness made. The reviewer is under the impression that these are Russian 
ideas, and that they are probably not very useful. Ways in which the machine 
could resolve ambiguities by determining the subject matter are given. For 
instance a word which could be either “resistance” or “‘opposition’’ would be 
referred back to an earlier frequency count. This count of key words would de- 
termine that the subject was electricity, in which case the first would be used, or 
politics, whence the choice might be the latter. Very polished translations might 
be made by having a large collection of literary quotations, from which phrases 
of just the right flavor could be selected; to succeed this would require large 
storage and sophisticated programming, and a fast machine. 

Braille II offers a fascinating opportunity to use the principles of mechanical 
translation against a simpler problem. Braille II uses a system of shorthand with 
a moderately complex set of rules, such as “do not use a contraction to bridge 
the components of compound words,” or, some pairs of words can be treated as 
one word ‘“‘where the sense permits.’’ These rules can be reasonably well met with 
a dictionary of a thousand entries. Unfortunately—from the viewpoint of a 
computer enthusiast—this transcription is now being done by volunteers gratis, 
leaving no opportunity for the expensive machine. 

Translation from French to English has been planned in some detail for 
APEXC. The small size of the memory permits only 250 French words in the 
vocabulary besides rules enabling the machine to recognize masculine from 
feminine, singular from plural, and most verb forms. Words are split into stems 
and endings. Not all the letters of a word need to be used ; for example in recevrai 
only rec is needed to identify the stem and only rai the ending. If the vocabulary 
were more extensive more complication would be needed. There are plans for 
46 verb endings. The authors assert that grammar and ambiguous endings are 
the only sources of real difficulty. No mention is made of how strange words will 
be handled. The procedures seem well thought-out, and sound as though they will 
work. We can look forward to seeing some examples of its output. 

The chapter on translating German into English is quite in contrast to that 
on French. It is 162 pages long, more than half the book. It is not based on a 
machine experiment but only on an armchair study. Again in contrast with the 
chapter on French it looks toward the difficulties of the future rather than the 
opportunities of the present. The study is based on Oswald and Fletcher's ‘‘Pro- 
posals for the Mechanized Resolution of German Syntax Patterns’”’ in Modern 
Language Forum XXXVI, 1951. Most of the examples of German used to illustrate 
the study are from a critique of Thomas Mann’s novels and tend to emphasize the 
difficulties to be encountered rather than the infrequency of these unusual forms. 
One concludes that this untutored machine is to translate Mann at his deepest 
immediately or it will be considered a failure. Oswald and Fletcher have listed 
82 syntactical units which can be used by a machine to analyze sentences. A 
second part of the study is the examination of a particular piece of text (7000 
words) dealing with the specific subject of electron optics. Counts were made of 
the vocabulary and the syntactical rules which would need to be invoked. The 
number of words needed was 960. Each case of multimeaning, where in different 
contexts different English words are required, is taken up. Most of them are pre- 
positions or adverbs. The variations in meaning are illustrated by auf dem Tisch, 
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auf dem Tanz, and auf dem Lande, in which auf is respectively ‘‘on,”’ ‘‘at,’’ and “‘in.”’ 

There is a resume of an article by I. S. Mukhin, Proc I. E. E. 103 (B), 463-472, 
1956, telling of a machine program for translating English into Russian. Trans- 
lation in this direction is more difficult than in the opposite, for the machine must 
synthesize the endings of the Russian words from the word order of the English. 
Each Russian word in the dictionary has 17 associated grammatical facts. After 
all these facts have been marshalled the Russian sentence is constructed. It may 
even happen that a word with its 17 grammatical properties is suppressed. If an 
English word is capable of more than one Russian equivalent, depending on con- 
text, this fact is recorded in the dictionary from which the machine (which is 
BESM) then has recourse to a second and more elaborate dictionary to resolve 
the ambiguity. Idioms are handled the same way. 

A very short chapter discusses the proposal of a meta-language or intermediate 
language understood only by the machine, the object being to simplify translating 
among many languages. With such a meta-language translation from any one of 
N languages to any other could be done by 2N programs. Without it N(N — 1) 
might be necessary. The authors point out that if instead of a meta-language one 
of the WN languages is selected to play a central role, then the number of programs 
is reduced to 2(N — 1), smaller yet. 

Last there is some discussion of a possible computer designed especially for 
translating languages. Armchair daydreaming about computers with peculiar 
properties is not without its rewards. It commonly brings out the fact that the 
modern computer is nearly optimal, and that abandoning any of its main prop- 
erties for other advantages is rarely attractive. In this case the authors suppose 
that since the dictionary will be changed only rarely the ability to selectively alter 
the memory can be abandoned (perhaps only in part, although this point is not 
clear) in favor of a very large store. Using glass discs very large dictionaries 
(several million words) with small access times are economically feasible. It 
appears later in the discussion that the authors would like to be able to compile 
concordances and word counts; this would require a selectively alterable store. 
With a serial or serial-parallel machine words of variable length are feasible and 
at the same time certain economies become possible. The authors also recommend 
replacing multiplication and division by other instructions. One such replacement 
would be a transfer of control conditioned on a comparison of two words. Another 
would be to locate and copy in the accumulator an entry from the dictionary. An 
interesting suggestion would enable the machine to insert an additional word into 
a list of words which was already monotonic non-decreasing. This would have the 
new word in the register ; it would then iterate the operation of copying the word 
at x into x + 1 and then comparing it with the word in the register. If the word 
in the register is smaller continue with address x decreased by one; if it is large or 
equal copy from the register into x + 1 and go to the next operation. 

Despite its uneven quality this book is an interesting exposition. It has a 
lively approach and promises exciting advances for the future. With enough effort 
machines can not only translate but can translate adequately, or even better. 


H. H. CAMPAIGNE 


National Security Agency 
Research and Development 
Fort George G. Meade, Maryland 
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114[X].—R. A. Buckincnam, Numerical Methods, G. P. Putnam’s Sons, New 
York, 1957, xii + 597 p., 21.5 cm. Price $15.00. 


This book, very comprehensive and very well written, is a combination of a 
reference book and a text book on standard numerical methods. The emphasis is 
on techniques and procedures, the proofs are as a rule omitted. The author gives 
a very detailed discussion of most of the methods available with abundant refer- 
ences to the literature and each important technique is followed by an example 
with all the necessary comments. Every chapter contains a collection of problems. 
Here is the list of topics: An introduction to computation ; Lagrangian methods of 
interpolation, differentiation and integration; Newton’s methods of divided dif- 
ferences ; Differences at constant interval ; Symbolic methods ; Solution of ordinary 
differential equations; Polynomial and other algebraic equations; Fitting of data 
by the method of least squates; Direct solutions of simulaneous linear equations; 
Matrices and determinantal equations; Indirect methods for linear equations; 
Ordinary differential and Fredholm integral equations; Functions of two vari- 
ables; and Partial differential equations. The book ends with useful appendices 
on: Relations between powers and factorials; Summary of difference formulas 
with remainder ; Lagrangian formulas for differentations and integration ; Orthog- 
onal polynomials for curve fitting; Index on useful formulas in text. 

The reviewer noticed the following omissions: a very incomplete discussion on 
ill conditioned matrices and an absence of warning that normal systems of equa- 
tions have a tendency to be ill-conditioned and are almost certainly ill-conditioned 
for systems of order six or greater. Other omissions, like the Runge-Kutta method, 
are mentioned by the author in the preface. 

This book, in the reviewer’s opinion, is probably the best on the subject and 
is highly recommended not only for mathematicians but for all scientists with 
any interest in computing. 

T. LEsER 


Ballistic Research Laboratories 
Aberdeen Proving Ground, Maryland 


115[X ].—N. Aronszajn, A. Douc.is, &C. B. MorrEy, jr., Editors, Transactions 
of the Symposium on Partial Differential Equations, Interscience Publishers, 
New York, 1956, vi + 334 p., 26 cm. Price $6.50. 


The Program of the Berkeley Symposium on Partial Differential Equations, 
June 20—July 1, 1955, had two parts. One consisted of twenty hours of expository 
lectures on linear elliptic equations, linear hyperbolic equations, non-linear partial 
differential equations, and approximations to solutions of differential equations. 
The other was composed of thirty-two 45-minute lectures on a wide variety of 
topics. Twenty-eight of these papers, first published in Communications on Pure 
and Applied Mathematics, v. 1X, No. 3, 1956, have been reissued in book form. A 
program of such length, with participants drawn from the leading centers of 
research in partial differential equations in the United States, not to mention 
contributions by a number of visiting scholars from abroad, would be expected 
to give a comprehensive picture of recent American acitivity in this field of mathe- 
matics. The book contains material to suit all interests, from casual beginners’ to 
specialized experts’, from proofs to conjectures. The subjects range from specific, 








252 REVIEWS AND DESCRIPTIONS OF TABLES AND BOOKS 


concrete, but delightful results obtained by simple, subtle devices from very 
special properties of the particular problems considered, to broad general conclu- 
sions derived by very abstract techniques. 

To discuss all twenty-eight papers adequately in a short review would be 
impossible. Readers interested in numerical analysis and applied mathematics 
would find material of fairly immediate usefulness and value in at least half of 
the papers. The present discussion will arbitrarily be restricted to these, and even 
they will be described essentially by titles. Three papers are explicitly devoted to 
numerical analysis. John Todd has given elegant systematic stability analyses 
for several methods for approximating parabolic and hyperbolic differential equa- 
tions. H. F. Weinberger and G. E. Forsythe have determined bounds for the 
eigenvalues of vibrating membrane problems in terms of eigenvalues for corre- 
sponding difference equations. L. E. Payne’s new isoperimetric inequalities repre- 
sent another approach to the basic problem of the two preceding papers. Many 
other papers deal with techniques and results useful for solving problems in 
applied mathematics. Noteworthy among these are A. E. Heins’ excellent account 
of the scope and limitations of the method of Wiener and Hopf and J. J. Stoker’s 
discussion of the determination of appropriate radiation conditions for steady 
state problems by passage for the limit in related nonsteady wave problems. 
Special representations of harmonic functions by S. Bergman, polyharmonic 
functions by A. Huber, and of singular and regular solutions of Cauchy problems 
for Euler-Poisson-Darboux equations by J. B. Diaz typify results that find fre- 
quent application. Information such as that about the nature of discontinuities 
of solutions of quasi-linear equations discussed by Y. W. Chen, about singularities 
of generalized axisymmetric potentials considered by Erdelyi, and several papers 
of estimates, growth, and local behavior of solutions has obvious importance for 
both pure and applied mathematics. In this connection, L. Bers’ very lucid sur- 
vey of local properties of solutions of elliptic partial differential equations de- 
serves especial mention. 


Joun H. GIEsE 
Ballistic Research Laboratories 
Aberdeen Proving Ground, Maryland 


116[X].—F. B. HrLpEBRAND, Methods of Applied Mathematics, Prentice-Hall, 
Inc., New York, 1956, xi + 523 p., 21 cm. Price $8.50. 


“The principal aim of this volume is to place at the disposal of the engineer 
or physicist the basis of an intelligent working knowledge of a number of facts 
and techniques relevant to four fields of mathematics which usually are not 
treated in courses of the ‘Advanced Calculus’ type, but which are useful in varied 
fields of application. The text includes the result of a series of revisions of material 
originally prepared in mimeographed form for use at the Massachusetts Institute 
of Technology.”’ The four fields are: Matrices, determinants and linear equations, 
Calculus of variations, Difference equations, and Integral equation; these are 
covered in chapters of about equal length. Each chapter is followed by an impres- 
sive collection of examples, about 100 to a chapter, and a list of references. 

Readers who are mainly concerned with numerical analysis will find much 
useful material in this book. Some is not where it might be expected, e.g., the 
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extremal properties of characteristic numbers of a matrix is discussed in the 
examples. 
The appearance of the book suffers from a very narrow margin. 
Joun Topp 


California Institute of Technology 
Pasadena, California 


117[Z ].—A. Gret, ABACS or Nomograms, Philosophical Library, New York, 
1956, ix + 225 p., 22 cm. Price $12.00. 


This is a translation of a book recently published in France and still available 
in a French edition [1]. Without sophistication it tells how to construct the sim- 
plest nomograms—Cartesian nomograms (essentially contours sketched on a 
plane with cartesian or other coordinates), alignment charts, and a iew others. 
There is no real mathematical treatment, and the work can be followed by anyone 
with knowledge of graphs similar to that obtainable in our secondary schools. 

The presentation is mainly by working out examples, each of which is stated 
as a formula pertaining to some shop or engineering problem. 

The utilitarian nature of the book is illustrated by an appendix consisting of 
several sections. One concerns the advantages of Cartesian nomograms, and is 
largely concerned with their retention of accuracy even when the material on 
which they are drawn is stretched, deformed, or torn; this permits photography 
without worry about shrinking of paper during development of the print, for 
example. A second section of the appendix concerns advantages of alignment 
charts—accuracy, ease of reading, ease of interpolation, and so on. Another sec- 
tion lists several common forms of relations (f/f: + fe + fs: = fifefs is one type 
listed, where the functions f; are to be taken into account by suitable choice of 
scales) with recommended types of charts. 

This is a useful text for those who want to construct nomograms, which may 
be the cheapest, the most convenient and the most widely used computational 
aids. Hence it has some considerable importance. 

As noted above, it has no theory or bibliographic references. Readers requiring 
more theory would probably be happier with any of several books listed in the 
bibliography below, say [9], [10], [7] and possible [1]. Some bibliographical 
references are found in all these and in [5]. A few helpful tables are included in 

[3]. Otherwise, the present book is similar in content to those listed below, and 
if it is chosen by a reader this choice must be based on arrangment of material, 
wealth of examples (and many of the others listed below use the method of exposi- 
tion by worked examples), or simplicity of exposition. 

While the reviewer does not object to the utilitarian nature of the book, he 
would have been more favorably impressed if a few examples of the greatest 
attainments in the art of nomography has been cited ; the common General Radio 
reactance charts relating inductance capacity, reactance and frequency for sinu- 
soidal voltages [11] and similar charts have saved electrical engineers enough 
time to pay for all the world’s research on nomographs, and the nomograms by 
Rybner [12] are examples of a professional nomographer of a high level of com- 
petence attacking a problem of seeming forbidding extent. It seems unfortunate 
that such examples are not made available to the student reading Giet’s book. 
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An abstract of the table of contents follows; there is no index. 


I. Relations Between Two Variables 


1. Diagrams 
2. Scales 


II. Cartesian ABACS 


3. Relations Involving Three Variables 

4. Superimposed Cartesian ABACS for Three Variables 

5. Superimposed Cartesian ABACS for Four Variables 

6. Superimposed Cartesian ABACS for Relations Involving Variables 
(n > 4) 


III. Alignment Charts 


7. Parallel Coordinates 
8. Graphical Representation of a Relation of the First Degree Between 
Three Variables 
9. Standard Forms of Alignment Chart for Relations Between Three 
Variables 
10. Charts for Four Variables 
11. Charts for Three Variables Requiring the Use of Auxiliary Variables 
12. Combined Charts for Several Equations Involving Common Variables 


IV. Alignment Charts Not Based on Parallel Coordinates 


13. Charts with Concurrent Scales 

14. Circular ABACS 

15. Double-Alignment Charts for a Relation Between Four Variables of the 
Form fife = fafa 


V. Relations Between » Variables (n > 4) 


16. Relations of the Form f: + fe+fs+---+f.=0 
17. Relations of the Form f; - fo- fs: + + fa-1 = fa 


Appendix 
Choice of Methods: Choice of Type of ABAC 


<= 
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12. J. RyBNER, Nomograms of Com ‘yperbolic Functions, Gjellerups, Copenhagen, 2nd 
edition, 1955. [MTAC, v. 11, 1957, Rev. 80, p. 207.] 
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118[Z].—T. E. Iva, Editor, Electronic Computers, Philosophical Library, Inc., 

New York, viii + 167 p., 22 cm. Price $10.00. 

This book is a non-mathematical introduction to the principles and applica- 
tions of computers employing tubes and other electronic devices. The treatment 
has been made very general in order to provide a broad background of the entire 
field of computing. 

The first chapter of the book discusses some of the outstanding contributions 
made to the computing field over the past three hundred years. 

The second chapter discusses the general principles of computing. Comparisons 
are made between pencil and paper, desk calculator, and electronic computer 
methods, in doing a particular type of mathematical operation. 

The next three chapters are devoted to analogue computers. Typical circuits 
of analogue computers are explained in general terms. Descriptions are given on 
how the circuits operate, how they are constructed, and combined to fulfill the 
requirements for a particular application. The types of inputs required, scaling, 
types of outputs and some of the applications for which an analogue computer 
are best suited are described. 

The remaining six chapters are devoted to digital computers. The basic arith- 
metic operations are discussed along with some of the different types of circuits 
required to perform these operations. Storage systems such as magnetic drums 
and discs, magnetic tape, magnetic core arrays, and cathode ray tubes are dis- 
cussed and comparisons are made between the relative speed of the different 
systems. 

The applications of digital computers to data processing, machine tool control, 
production planning are discussed in general terms. 

The book gives the reader a painless introduction to the principles of operation 
and application of high-speed electronic computers and computer components. A 
background in mathematics or a previous knowledge of computer circuitry is not 
required for a thorough understanding of the contents. In clear and concise terms, 
the author covers the subject in eleven well-written chapters. The book is written 
primarily for electronic technicians and students. Engineers in other fields who 
wish to become acquainted with the subject of computers will find the text both 
interesting and informative. The book lays a firm foundation for subsequent 
detailed reading in the digital computer field. Persons in the sales, managerial, 
administrative, and executive fields would find it quite helpful. It is not intended, 
however, that the book be of use to engineers who are well informed in computing 
techniques and equipment. The excellent illustrations of components and cir- 
cuitry serve well as an aid to comprehension. 


CuHEsTER H. WALLIN 
Computing Laboratory, BRL 
Aberdeen Proving Ground, Maryland 


119[Z].—R. K. Livestey, Digital Computers, Cambridge University Press, 
London, 1957, viii + 53 p., 21 cm. Price $1.75. 
This slim volume of 53 pages is published in the Cambridge engineering series. 
It amounts to a brief introduction to digital computers for engineers who are 
prospective users of machine. The author starts by analyzing a simple calculation 
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(the solution of three simultaneous linear algebraic equations in three unknowns) 
to show what arithmetical and logical steps are involved. He then describes a 
simple computer which could carry out this calculation. Next a brief and sketchy 
treatment is given of the physical equipment needed for memory and input- 
output units ; no attempt is made to show how control and information processing 
functions are actually carried out. There follows a chapter on programming, 
stressing the utility of subroutines and showing how they can be organized into 
complete programs. The book is concluded by a chapter giving a brief sketch 
of typical engineering problems requiring extensive computations. 

The writing is clear, unpretentious, and easy to follow. The book could be 
understood, with a little additional explanation, by an interested reader of limited 
scientific background. 

eT ten 


120[Z ].—G. A. MontTGoMERIE, Digital Calculating Machines, Blackie and Son 
Limited, Glasgow, Scotland, 1956, vii + 262 p., 22 cm. Price 30s. net. 


This is an elementary practical guide without pretentions of scholarship and 
without efforts to remanticize computers. The material covered includes all com- 
mon types of computers, but much is sparse, particularly with regard to larger 
computers. Illustrations of common machines are included, but nothing like com- 
plete coverage of machines available (particularly those easily available outside 
but not within the United Kingdom) is attempted. 

Additional references are given at the end of each chapter; the author’s view 
that reference material for most topics is not useful or is hard to obtain is not 
entirely shared by the reviewer, but, on the other hand, it seems to the reviewer 
that by the time a reader has advanced to a point where-he needs such references 
they may become known through other contacts. In any case the bibliographic 
references are not complete. 

In short, what is put together is an elementary work giving the first principles 
of operation of the various types of machines discussed. It is a practical work, 
although it does not group and exhibit recipes for particular calculations. Square 
roots are taken by repeated subtraction on page 95, for example, with a suggestion 
that one of the methods described later in the book be used for most work. On 
page 117, square roots are again mentioned and logarithms are suggested. Finally, 
on page 141 the common iteration, x, = $(x,-1 + @/x,_1) is introduced as the 
best way to find a!, but never is there given a rule for continuing by this method 
an extraction started by successive subtraction, which is easily brought about on 
some machines. Still, such methods are usually (and properly) brought to users’ 
attention by the manufacturers of the various machines, and description in a 
small book such as the present with no reference to a laboratory instrument would 
entail difficult exposition. 

There is an index of processes treated, on pages 143-144, and this has some of 
the uses of a set of recipes. 

Analysis of strong and weak points of machines and their applications like 
that in [1] is missing almost entirely. The Mechanical principles are treated, but 
not deeply, or to the extent some are described in [2 ]. 

Thus, what is put together is a sound but not advanced text, neither deep nor 
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encyclopedic in nature, but containing material not easily available elsewhere in 
English. It is mainly useful in connection with manually operated machines, here 
mainly presenting a description for those who have not yet used such machines 
to any great extent. It describes many but not all types of machines. There is 
undoubtedly a set of readers who will find such a text more useful than [2], [3], 
[4], [5] or [6], for example. 

The table of contents follows: 


Chapters 


I. Introduction 
II. Key-responsive Adding Machines (Group I, Type 1) 
III. Key-set Adding Machines (Group I, Type 2) 
IV. Four-rule Machines, Lever-set, Hand-operated (Group I, Type 3) 
V. Four-rule Machines, Key-set and Electrically-operated (Group I, 
Type 3, cont.) 
VI. Four-rule Automatic Machines (Group I, Type 4) 
VII. Computing with Four-rule Machines 
VIII. Punched-card Machines (Group II) 
IX. Special-purpose Machines (Group III) 
X. Universal Sequence-controlled Machines (Group IV) 
XI. Electronic Sequence-controlled Machines (Group IV cont.) 
XII. Programming of Sequence-controlled Machines 
ee TF. 
L. J. Comrie, “The cn of commercial calculating machines to scientific computing,” 
= “AC v. 2, 1946, p. 149-159. 
F. J. Murray, The ‘Theory of Mathematical Machines, King’s Crown, New York, 1947 and 
(revised) 1948. 
. F. A. Witters, Mathematische Maschines and Instrumente, Akademie, Berlin, 1951. 
r F. A. WIters, Mathematisch Instrumente, Gruyter, Berlin, 1926. 


5. W. MEYER zuR CAPELLEN, Mathematische Instrumente, Akademie, Leipzig, 1949. 
6. E. HErsBurGH, Modern Instruments and Methods of Calculation, Bell, London, undated. 


121[Z].—Epmunp C. BERKELEY & LAWRENCE WAINWRIGHT, Computers, Their 
Operation and Applications, Reinhold Publishing Corporation, New York, 
1956, ix + 366 p., 22 cm. Price $8.00. 


The authors admit in their preface that it is a hopeless task to gather all the 
information about computing machinery. Their approach is to “do the best one 
can to report information already learned—incomplete as it is—indicating areas 
where other information is needed.”” The appeal is to the general reader, and it is 
truly unfortunate that the book is so disorganized as to be confusing to the novice 
in the field. 

One assumes from the first chapters that the approach is for the non-expert, 
yet the section on analog equipment (written by the second author) does not 
hesitate to use derivatives and other mathematical operators with the assumption 
of previous knowledge on the reader’s part. 

The peculiar mixture of the general and the specific which is the hallmark of 
the first author leads him to make some misleading and some inaccurate state- 
ments. A few quotations will perhaps suffice to show what is meant. 

On the operation of overflow he states on page 41: ‘“‘Another useful operation 
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is overflow, in a counter; this is the production of a number in the counter which 
is beyond its capacity, giving rise to a carry signal that can be employed for useful 
purposes. This would happen for example in a six-place counter holding 999784 
if 4000 were added.” 

Unfortunately, he says nothing more about the usefulness of overflow, nor 
what purposes the carry signal could be employed for. In addition, the quotation 
serves to mar the distinction between adders and counters, which he does not 
clarify elsewhere. 

On accuracy he states on page 58: “‘An automatic digital computer is sur- 
prisingly accurate. One of the earliest machines, the Harvard IBM Mark I, repre- 
sented numbers of 23 decimal digits. . ” This, by the author’s own defi- 
nitions in the glossary appended to the book, is not accuracy, but precision. A 
few sentences further on he repeats the error: ‘“‘Accuracy to such a high degree 
is not often useful.’’ Again it is precision which is meant. Finally to conclude this 
exhibition, and make “‘confusion worse confounded,’”’ we have: ‘‘One machine, 
the IBM Electronic Selective Sequence Calculator, . . . calculated the motion 
of the moon with new and remarkable accuracy, locating the moon precisely in the 
heavens at six-hour intervals for four centuries.’’ (Italics are the reviewer's.) 

In Section IV he deals with types of automatic computing machines that are 
not digital or analog computers. In drawing this artificial distinction he only 
succeeds in completely confusing the reader. He lists, among other non-computers 
on page 162, an astronomical telescope aiming equipment to adjust the direction 
of a telescope to follow a particular star field. This is as clear an example of an 
analog computer as one would wish, and the author himself says so on page 293: 
“‘And automatic analog computers are responsible for keeping a telescope pointed 
at a heavenly body while the astronomer watches it through the eyepiece.” 

It would have been far better to avoid such artificial distinctions if he was not 
able to be consistent in his definitions. Such loose thinking could be demonstrated 
by many more references to the text. 

Section III, written by Lawrence Wainwright, is, by contrast, ‘‘a good deed 
in a naughty world.” It is a clear elementary exposition of the functions and ele- 
ments of analog computers, with examples of various types. The book would be 
worthwhile if it were all on a par with this section. 

Section VII, on applications, spends a good deal of space on arguments as to 
the usefulness of computers in fields where they have long made their mark. The 
actual descriptions of applications in some fields are very terse. In mathematics, 
the whole field of numerical analysis is ignored. Computers, according to this 
author, have determined whether numbers are prime, and calculated the values 
of x and e to 2000 decimal places. The basic importance of being able to ask the 
computer the right questions so as to get the right answers is nowhere stressed. 
Even the fact that computers have opened a whole new realm of unsolved prob- 
lems is glossed over. 

The book concludes with a section of references, a roster of organizations and 
computing services, and a glossary. The index reveals no mention of the name of 
von Neumann, an omission which the author should find it hard to justify. 

I have rarely read a book which could make me angry by the mere style of 
the author; this one does. I have tried to approach it with the objectivity of a 
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person new to the subject, as a general reader might. Even this did not help; the 
fuzzy thinking obscured all the difficult parts, and the “‘cute’”’ analogies which 
compare a computer to a large railroad yard with switching stations are of no 
avail in the serious business of understanding what a computer really is. 


F. HOLLANDER 
University of California 
Los Angeles, California 


122[Z].—Georce F. Fores, Digital Differential Analyzers, An Applications 
Manual for Digital and Bush Type Differential Analyzers, Third Edition, 
Published by the author, Pacoima, California, 1956, xii + 154 + 19 p. of 
Appendices and 9 p. of index sheets, bound, soft cover. Price $7.50. 


[Z ].—GEorGE F. Forses, Digital Differential Analyzers, An Applications 
Manual for Digital and Bush Type Differential Analyzers, Fourth Edition, 
Published by the author, Pacoima, California, 1957, iii + iv + 154 + 22p. 
of Appendices and 8 p. Index sheets. Price $5.00. 


This book is strictly as described in the title, an applications manual. A very 
brief introduction to the principles of digital differential analyzers is given in 
Chapters I and II, not sufficient to explain their operation fully to one 
unacquainted with devices of this type. However, this is not a serious defect, 
since the detailed information on the use of these machines contained in the 
remaining twenty-seven chapters will only be of interest to those who really 
understand them anyway. One’s general impression is that the author presents 
this material carefully and precisely, and that it must be very useful to the 
audience to whom it is addressed. 

The Third and Fourth Editions are essentially identical. The author states 
that some few corrections have been made in the Fourth, while the only additions 
are those of a few extra Exercises and a few extra items in the Bibliography. 


Cc. Vv. L. S. 


123[Z].—O. J. Dau, Autocoding for the Ferranti Mercury Computer (MAC), 
NDRE Report No. 24, Oslo, Norway, 1957, iii + 92 + 11 p. of appendices. 


This is the description of the external features of an algebraic compiling sys- 
tem similar to FORTRAN written for the Ferranti Mercury Computer. In addi- 
tion to the floating-point representation of numbers which is built into the 
machine, the system allows arithmetic operations on complex numbers, double- 
length floating-point numbers, and short integers. Programs are written as 
sequential lists of directives, or instructions, which may themselves be elaborate 
algebraic formulas composed of pre-defined variables and constants. For example, 
the directive to compute 


x/2(1 + x*)§ ++ $1n (x + (1 + x*)4) 
is programmed as 
(X/2)Q(1 + X*2) + (/2)IEXP(X + Q(1 + X*2))-Y 
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where X denotes a single variable previously defined, Q the square root, IEXP 
the inverse exponential, * the power, + addition, and / division. Functions avail- 
able directly in the symbolism are exponential, sin, cos, tan, square root and their 
inverses as well as more complicated functions such as integral and derivative. 
Other functions can be incorporated as subroutines. 

The built-in B-registers are not referred to directly; a postfix method of 
indexing sets of variables is used. For example, the matrix element a;; is written 
as An12 where index nl is used to count rows and n2 to count.columns. For all 
programs the extents of fields of variables must be defined in a set of field founda- 
tions preceding the instructions. 

The Mercury high-speed store consists of 16 pages of 32 floating-point num- 
bers each. This is backed up by a magnetic drum from which a page at a time 
may be transferred. Large matrices, for example, must be stored on the drum and 
transferred as needed to the high-speed store. Matrices are frequently filled out 
to multiples of 32 by adding zero elements. A field in the H. S. store is referred to 
by a symbol such as FA, on the drum as FMA. Operations on an entire field in 
the H. S. store are possible without counting. For example, the directive 0-FA 
sets the whole field FA to zero. Facilities are also available for constructing com- 
pound logical statements. In the report, programs are shown for various matrix 
and vector operations, the solution of linear equations by the Jordan-elimination 
method, and there is a discussion of indexing for non-rectangular fields. 

The normal method of looping, using counters and conditional jump instruc- 
tions, can be used but a more sophisticated method, using operators, is also 
described. These operators cause the sequence of directives within their range of 
influence to be executed systematically for all possible values of an index stated 
in the operator. r 

The system is said to prepare programs comparable to manually coded ones. 
Storage allocation of program and data is still the responsibility of the programmer 
if these cannot all be contained in the H. S. store. No indication of experience with 
the system is given nor is there very much information about the mechanism of 
compilation. It appears to be a comprehensive scheme admirably suited to the 
needs of a seasoned programmer. The more elementary sections of the report 
would undoubtedly permit even the beginner to use the system. 


J. N. P. Hume 
University of Toronto 
Toronto, Canada 


124[Z].—O. J. Dani, Multiple Index Countings on the Ferranti Mercury Com- 
puter, NDRE Report No. 23, Norwegian Defence Research Est., Oslo, 
Norway, 1957, iv + 43 + 10 p. of diagrams. 


This report describes a systematic method of handling multiple countings in 
programs for problems such as linear algebra or tensor calculus. The form of 
counting depends on how a set, or field, of data is stored in the machine. The 
dimensionality of a field stored completely in the high-speed store is the number 
of indices describing an element of the field, e.g., 1 for a vector, 2 for a matrix. 
The range of an index is known as an extent of the field. In general an extent may 
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be a function of other indices; fields whose extents are constant are rectangular 
fields. The H. S. store of a computer is essentially a one-dimensioned field and 
the problem as stated is ‘‘to store an n-dimensioned field in a linear one.’’ For 
example, matrices can be stored as sets of linear subfields, end-to-end, e.g., row 
by. row. Use of the secondary store, the magnetic drum, can increase the number 
of indices. For example, if a single column of the matrix occupies several sectors 
of the drum, countings must proceed both on sectors and elements within sectors. 
Various methods of partitioning matrices are discussed. Flow diagrams and 
Mercury Computer Programs are given for a variety of partitioning methods. In 
fitting a field into the machine storage, when it cannot be divided into an integral 
number of similar subfields, zeros may be added to complete subfields. It is found 
best to complete at least all H.S. store subfields with zeros allowing for incom- 
pleteness on the drum only. Countings may avoid operations on the zeros but the 
gain in time for large matrices appears to be offset by the higher degree of com- 
plexity in the innermost countings. 

J. N. P. Hume 
University of Toronto 
Toronto, Canada 


125[Z ].—LorENzo LUNELLI, “Un comando di computo di cicli per una calco- 
latrice elettronica aritmetica a tre indirizzi,’’ La Ricerca Scientifica, v. 27, 
1957, p. 3381-3394. 


This paper describes a ‘‘cycle computing instruction’’ added to the CRC-102 A 
digital computer located at the Centro di Calcoli Numerici of the Milan Poly- 
technical School. Essentially, the command combines counter modification (by 
subtraction of binary one from specified sections of the contents of a counter cell— 
any desired cell of the memory) with branch point on zero results in any section 
so specified. The feature of interest in the command is the implicit capability of 
combining several separate counts in a single cell with branching possible on each 
such count. The command is mechanized as follows: Let the command be repre- 
sented by Cmymzm;. Then, m, is the counter cell; (m2) specifies the division of 
(m,) into separate counters by locating contiguous binary ones in the correspond- 
ing bit positions of (mz) (individual counters are therefore separated by at least 
one bit) ; the command subtracts one from the least significant bit of each of those 
sections of (m,) for which (m2) has such a contiguous sequence of binary ones and 
replaces (m,) by the result ; m; specifies the location of the next command if none 
of the sections so specified becomes all zero; otherwise the next command is taken 
from the address following that of the command being executed. The use of this 
command is particularly effective for address modification in a repetitive sub- 
routine. The paper includes illustrations of this type of operation. 

An appendix to the paper describes a proposed command with somewhat 
similar purpose but different mechanization. The proposed command is similar 
in function to the familiar ‘‘B-box’’ except that it includes the modification of the 
B-box with the actual augmentation of the command, rather than with the 
branching decision. The proposed mechanization is as follows: Let the command 
be Cmymm;; then m, is the location of the command to be executed ; mz is the 
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““B-box,”’ and its contents modify (mm) ; (m3) is the modifier of (mz) ; the command 
results in the execution of the command stored in m, as modified by (mz) and in 
the replacement of the original contents of m2 by the result of (m2) + (ms). 

No time estimates are given for either of the two commands. However, the 
first one appears to be fairly cheap, requiring perhaps one or two word times 
more than a normal subtract command. On the other hand, the second, proposed, 
command requires four memory look-ups in addition to those of the command 
being modified and executed as part of the operation. Under many circumstances 
this would be very costly in time. These commands are quite representative of 
some of the complex built-in sub-routines which can fairly easily be added to a 
micro-programmed machine such as the CRC-102A. Such built-in sub-routines 
are much faster than the corresponding programmed sub-routine. On the other 
hand, it is characteristic of the micro-programmed approach that the resulting 
built-in sub-routines are slower than they would be if designed otherwise; the 
advantage of micro-programming is the relatively efficient use of equipment. 


R. M. Hayes 
Business Systems Department 
Magnavox Research Laboratory 
Los Angeles, California 


TABLE ERRATA 


265. Epwin P. Apams, Smithsonian Mathematical Formulae and Tables of Elliptic 
Functions, Second reprint, The Smithsonian Institute, Washington, D. C., 
1947. See also MTAC, v. 1, p. 191, 325; v. 2, p. 46, 352-3; v. 3, p. 314, 423; 
v. 6, p. 236. . 


P. 122, formula 6.42, no. 4, the sign before the >> should be —. 

P. 127, formula 6.475, no. 2, the formula is correctly printed in this reprint. See 
MTAC, v. 2, p. 46. 

P. 139, formula 6.821, in the numerator of the fraction in the second summation 
(involving sines), for 1, read n. 

P. 139, formula 6.822, the formula should read 
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E. G. H. Comrort 
Ripon College 
Ripon, Wisconsin 


266. F. J. MAssEY, jr., ‘Distribution table for the deviation between two sample 
cumulatives,’’ Annals of Math. Stat., v. 23, 1952, p. 435-441. 


The entire table was recomputed using the IBM 701 electronic computer at 
the University of California’s Computer Center, Berkeley, California. In addition 
to the following 46 errors there are approximately twenty other entries in which 
the last digit has been rounded incorrectly. The 16 corrections marked with an 
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asterisk were previously reported by J. L. Hodges, Jr., in Arkiv for Matematik, 


1958, in print. 
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University of California 
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-00833 
-99126 
-99467 
-18943 
65703 
-71795 
-77762 
-93629 
83566 
01518 
68511 
-98416 
-98900 
99393 
-67045 
-90123 
92317 
35336 
-96564 
99281 
-70465 
-99448 
99637 
.99731 
99799 
99864 
82837 
-72675 
-76664 
57501 
-76578 
53765 
.99423 
99624 
73908 
.00020 
97431 
-98093 
.98554 
-98915 
.99216 
-96528 
98044 
98586 
.94026 
12214 


Read 


a 


01677 
.99184* 
.99534* 
-20140 
.68687* 
-74779* 
.80746* 
.94406* 
86014 
01598 
.67935 
.98601* 
.99054* 
.99548* 
-70447 
.90489* 
.92957* 
36252 
.96973* 
.99296* 
-73189 
.99519 
.99707 
99801 
-_99870 
-99934 
82715 
-73230 
-17734 
57581 
-75919 
55837 
.99475* 
.99669* 
73811 
00197 
97511 
.98173 
98635 
98996 
-99298 
.96531 
98101 
.98637* 
.94023 
12114 
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NOTES 
Association Francaise de Calcul 


The editors have been informed of the founding of the Association Francaise 
de Calcul. The first President is the well known astronomer Andre Danjon, of the 
French Academy of Sciences. 

The Association intends to develop all aspects of numerical analysis, comput- 
ing and associated fields ; its members will include people inierested in these topics 
from any profession—engineers, economists, statisticians and so on. 

The address of the new Association is: 

Association Francaise de Calcul 
Institut d’Astrophysique 
98 bis, Boulevard Argo 
Paris 14, France 
ae 2 


Milton Abramowitz (1915-1958) 


Dr. Milton Abramowitz, head of the Computation Laboratory of the Applied 
Mathematics Division at the National Bureau of Standards, died suddenly of a 
heart attack at his home on Saturday, July 5. 

Dr. Abramowitz was a valued younger scientist of the Bureau. Starting in 
1938 as an original member of the Technical Planning Staff of the New York 
Mathematical Tables Project, the scientific program of which was under the 
direction of the Bureau, he moved upward to become in 1954 the Chief of the 
Bureau’s Computation Laboratory. In addition to the supervision of the Labora- 
tory’s computation facilities, he directed the preparation of mathematical tables 
and performed research in numerical analysis and mathematical physics. Among 
his many accomplishments, he is known for his significant contribution to the 
planning and preparation of the well-known volumes of basic mathematical tables 
published in the Mathematical Tables Series of the original New York group and 
in the Applied Mathematics Series of the Bureau, and most recently for his 
editorship of a comprehensive Handbook of Mathematical Functions, which the 
Bureau is preparing under the sponsorship of the National Science Foundation. 

Born in Brooklyn, New York in 1915, Dr. Abramowitz received his B.A. in 
mathematics in 1937, and his M.A. in 1940, both from Brooklyn College. He 
earned his Ph.D. from New York University in 1948. 

Dr. Abramowitz received the Gold Medal in mathematics from Brooklyn 
College in 1933, and was awarded the Meritorious Civilian Service Award from 
the U. S. Navy in 1945. He was a member of the American Mathematics Society, 
and the Washington Academy of Sciences. 


Horace Scudder Uhler 
1872-1956 
By J. W. Wrench 


The art of highly accurate computation lost one of its most gifted exponents 
and practitioners when Professor Horace S. Uhler died in Meriden, Connecticut 
on December 6, 1956, following a brief illness. 
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Professor Uhler was born in Baltimore, Maryland on August 5, 1872, and, 
following his education as a physicist, he spent all his adult life in the academic 
world. He received a baccalaureate degree from Johns Hopkins University in 
1894, following which he was an instructor in the Academy of Northwestern 
University, laboratory director at the University School in Baltimore, and gradu- 
ate student and fellow at Johns Hopkins, where he received his doctorate in 1905. 

The following year he served as Carnegie assistant in physical chemistry at 
Johns Hopkins University, whence he was appointed instructor in physics at 
Yale College, Yale University. This was the beginning of an association that was 
to extend throughout the rest of his life, except for an interlude, 1925-26, when 
he was the director of the department of physics at Gettysburg College. In 1937 
he received an honorary master’s degree from Yale University, and concurrently 
he was elevated to the rank of professor of physics. 

The title of Professor Emeritus of Physics was officially conferred on Professor 
Uhler on the occasion of his retirement in 1941. Enlistment of many of the younger 
members of the faculty during the second World War led to his service as a special 
lecturer in physics for a period of two years. After returning to retirement in 1943, 
he devoted all his energies to his avocation, mathematical computation, which to 
a progressive extent had occupied his attention during the preceding half-dozen 
years. 

He was pre-eminent in the field of elaborate computation because of his scru- 
pulous attention to each detail and because of his extraordinary efforts to insure 
the impeccable quality of his published results. Professor D. H. Lehmer once re- 
ferred to Uhler’s ‘indefatigable perspicacity” in describing his practice of subject- 
ing all his publications to the most careful scrutiny. This habit led to the occa- 
sional discomfort of editors who did not always meet his typographical standards 
of perfection! 

Uhler’s numerical research focused principally on factorials and their recipro- 
cals, natural logarithms, and Mersenne numbers, the computations being per- 
formed by means of electrically operated desk calculators. 

His exact values of factorials, which appeared in a succession of publications, 
and his elaborate table of reciprocals of factorials to 475D remain unrivalled to 
the present time. A definitive compilation of radix logarithms is his 137D table of 
such data, privately published in 1942. Supplementary material is available in a 
series of papers giving approximations (ranging from 155D to 330D) to the 
natural logarithms of small primes, as well as to the modulus of the common 
system of logarithms and to z. 

At the suggestion of Professor R. C. Archibald he began in 1944 the pioneer 
investigation of the primality of the six remaining Mersenne numbers 2? — 1, 
corresponding to p = 157, 167, 193, 199, 227, and 229, for which no identification 
was then available. His numerical proofs that these integers are all composite 
were confirmed by subsequent calculations performed on the National Bureau of 
Standards Western Automatic Computer (SWAC). 

A sequel to this investigation was the computation of the first seventeen 
perfect numbers, following the identification by SWAC of all Mersenne primes 
for exponents p not exceeding 2303. 

A chronological list of Professor Uhler’s significant mathematical publications 
is appended. Copies of certain unpublished manuscripts have been deposited in 
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the library of Brown University. These include a table of factorials of all integers 
from 201 to 300, inclusive, and a table of factorials decomposed into prime factors 
for consecutive integer arguments ranging from 1001 to 3400, inclusive. The latter 
constitutes an elaborate extension of a table of M. Alliaume published in 1928. 

In addition to his numerous papers in mathematics, Professor Uhler contrib- 
uted a number of papers in the fields of geometrical optics, spectroscopy, and 
x-rays. 

He was a charter member of the Mathematical Association of America, a 
member of the American Mathematical Society, the London Mathematical 
Society, the Optical Society of America, a life member of the Société Frangaise 
de Physique, and a fellow of the American Association for the Advancement of 
Science, and the American Physical Society. He was a former associate editor of 
the American Mathematical Monthly, the Physical Review, the American Journal 
of Science, and the Journal of the Optical Society of America. 

A portrait of Professor Uhler, accompanied by biographical and bibliographical 
information, appears in the book Mathematical Table Makers, by R. C. Arch- 
ibald, published by Scripta Mathematica, New York, New York, 1948. 

The indelible impression left on the memory of the writer, who enjoyed 
Professor Uhler’s friendship for the past two decades, was that of a genuine scholar 
of broad interests and attainments, who set for himself the highest standards of 
clarity and accuracy in his research. 


1. “On the numerical value of i*,” American Mathematical Monthly, v. 28, 1921, p. 114-116. 

2. “A new table of reciprocals of factorials and some derived numbers,” Conn. Acad. of Arts 
and Science, Trans., v. 32, 1937, p. 381-434. 

3. “Log w and other basic constants,” Nat. Acad. Sci., Proc., v. 24, 1938, p. 23-30. [See also 
MTAC, v. 1, 1943, RMT 95, p. 55.] 

4. “Recalculation and extension of the modulus and of the logarithms of 2, 3, 5, 7, and 17,” 
National Acad. Sci., Proc., v. 26, 1940, p. 205-212. [See MTAC, v. 1, 1943, RMT 96, p. 56.] 

5. “The coefficients of Stirling’s series for log T(x),’”” Nat. Acad. Sci., Proc., v. 28, 1942, p. 
59-62. [See MTAC, v. 1, 1943, RMT 97, p. 56. 

6. Original Tables to 137 Decimal Places of Natural Logarithms for Factors of the Form 
1 +n-10-?, Enhanced by Auxiliary Tables of Logarithms of Small Integers, New Haven, Conn., 
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